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High  throughput  (combinatorial)  materials  science  methodology  is  a  relatively  new  research 
paradigm  that  offers  the  promise  of  rapid  and  efficient  materials  screening,  optimization,  and 
discovery.  The  paradigm  started  in  the  pharmaceutical  industry  but  was  rapidly  adopted  to 
accelerate  materials  research  in  a  wide  variety  of  areas.  High  throughput  experiments  are 
characterized  by  synthesis  of  a  “library”  sample  that  contains  the  materials  variation  of  interest 
(typically  composition),  and  rapid  and  localized  measurement  schemes  that  result  in  massive  data 
sets.  Because  the  data  are  collected  at  the  same  time  on  the  same  “library”  sample,  they  can  be 
highly  uniform  with  respect  to  fixed  processing  parameters.  This  article  critically  reviews  the 
literature  pertaining  to  applications  of  combinatorial  materials  science  for  electronic,  magnetic, 
optical,  and  energy-related  materials.  It  is  expected  that  high  throughput  methodologies  will 
facilitate  commercialization  of  novel  materials  for  these  critically  important  applications.  Despite 
the  overwhelming  evidence  presented  in  this  paper  that  high  throughput  studies  can  effectively 
inform  commercial  practice,  in  our  perception,  it  remains  an  underutilized  research  and 
development  tool.  Part  of  this  perception  may  be  due  to  the  inaccessibility  of  proprietary  industrial 
research  and  development  practices,  but  clearly  the  initial  cost  and  availability  of  high  throughput 
laboratory  equipment  plays  a  role.  Combinatorial  materials  science  has  traditionally  been  focused 
on  materials  discovery,  screening,  and  optimization  to  combat  the  extremely  high  cost  and  long 
development  times  for  new  materials  and  their  introduction  into  commerce.  Going  forward, 
combinatorial  materials  science  will  also  be  driven  by  other  needs  such  as  materials  substitution 
and  experimental  verification  of  materials  properties  predicted  by  modeling  and  simulation,  which 
have  recently  received  much  attention  with  the  advent  of  the  Materials  Genome  Initiative.  Thus, 
the  challenge  for  combinatorial  methodology  will  be  the  effective  coupling  of  synthesis, 
characterization  and  theory,  and  the  ability  to  rapidly  manage  large  amounts  of  data  in  a  variety  of 
formats.  ©  2013  AIP  Publishing  LLC.  [http://dx.doi.org/10.1063/L4803530] 
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I.  INTRODUCTION  AND  HISTORICAL  PERSPECTIVE 

High  throughput  (combinatorial)  materials  science  is  a 
relatively  new  paradigm  for  materials  discovery  and  optimi¬ 
zation  that  arose  from  the  pharmaceutical  industry1-4  and  rap¬ 
idly  spread  to  other  fields,  several  of  which  have  matured  and 
been  the  subject  of  review  articles:  organic  coatings,5  sensor 
materials,6  inorganic  materials,7  and  functional  materials.8-10 
The  idea  of  utilizing  deliberately  created  composition  gra¬ 
dients  in  co-deposited  thin  films  to  obtain  significant  portions 
of  phase  diagrams  can  be  traced  back  to  about  fifty  years  ago 
for  various  ternary  systems.11-14  In  this  review,  we  concen¬ 
trate  on  applications  of  combinatorial  materials  science  for 
electronic,  magnetic,  optical,  and  energy-related  materials. 
As  the  emphasis  of  this  review  is  on  applications  of  these 
functional  materials,15’16  library  synthesis  techniques,  and 
properties  measurement  and  analysis  methods,  are  only  cov¬ 
ered  as  they  apply  to  the  materials  of  interest.  Further,  some 
high  throughput  fields  of  great  importance,  such  as  discovery 
and  optimization  of  catalysts  (well  accepted  and  integrated 
into  the  field),  are  largely  omitted  from  this  review.  In  con¬ 
trast  to  many  of  the  topics  covered  in  this  review,  high- 
throughput  catalyst  development  is  perhaps  more  prevalent  in 
industry  than  in  academia;17-21  a  number  of  critical  reviews 
and  perspectives  have  been  published,2’22-27  which  we  rec¬ 
ommend  to  the  reader.  Another  area  where  combinatorial 
strategies  have  played  a  significant  role  in  advancing  the  field 
is  in  polymeric  and  soft  materials,  and  recent  reviews  of  this 
field  are  available28’29 

The  combinatorial  materials  science  approach24’30-33  is 
characterized  by  creation  of  a  “library”  sample  that  contains 
variations  of  the  materials  parameter  of  interest  (typically 
composition,  although  this  paper  will  discuss  many  others), 
and  rapid,  local,  and  automated  measurements  to  interrogate 
the  library  for  the  property  of  interest,  resulting  in  the  gener¬ 
ation  of  large  data  sets.  These  data  sets  are  analyzed  and  gen¬ 
erally  displayed  as  a  response  surface  (i.e.,  property  as  a 
function  of  variable  of  interest). 

Combinatorial  methodologies  are  generally  faster  and 
less  expensive  (likely  tme  even  taking  initial  equipment 
investments  into  account)  than  traditional  (e.g.,  “one  composi¬ 
tion  at  a  time”)  approaches.  Further,  since  all  experiments  are 
carried  out  on  the  same  library  sample  with  the  same  measure¬ 
ment  tool  over  a  short  time  period,  most  systematic  errors  are 
eliminated,  thus  generating  a  comprehensive  and  reliable  data 
set.  Combinatorial  materials  science  excels  in  the  observation 
of  trends  and  in  screening,  i.e.,  finding  the  “sweet  spot”  of  the 
response  surface,  which  enables  rapid  materials  optimization. 

In  spite  of  the  large  number  of  papers  published  in  the 
field  of  combinatorial  materials  science,  industry  has  not 


adopted  the  new  paradigm  on  a  large  scale,  and  there  has 
been  some  skepticism.  Some  of  this  is  due  to  the  nature  of 
the  library  films,  which  may  contain  local  variations  such  as 
microstructural  changes  and  stress  variations  that  result 
directly  from  compositional  variations;  thus,  the  libraries  are 
perceived  as  flawed,  in  that  although  only  one  variable  is 
intentionally  varied,  other  properties  may  change  uninten¬ 
tionally,  resulting  in  a  perceived  “unscientific”  approach. 
Other  reasons  for  skepticism  are  based  on  practicality,  such 
as  the  high  cost  of  a  tool  set  for  synthesizing  library  samples, 
the  paucity  of  commercially  available  high  throughput  mea¬ 
surement  instruments,  and  the  highly  skilled  labor  needed  to 
run  a  combinatorial  laboratory.  While  it  is  difficult  to  fully 
defend  the  notion  that  the  composition  of  a  library  film  can 
be  varied  while  maintaining  all  other  properties  constant,  the 
use  of  a  greater  number  of  analysis  techniques  should  allow 
these  complicated  interactions  to  be  sorted  out.31 

Another  source  of  criticism  is  the  frequent  necessity  to 
process  the  combinatorial  library  samples  using  different 
growth  and  post-growth  processing  techniques  than  will  be 
employed  in  the  production  of  the  material  for  its  intended 
application.  This  criticism  can  be  mitigated  to  some  extent 
by  realizing  that  in  combinatorial  analysis,  the  “slope”  of  the 
property  with  regard  to  the  parameter  variation  may  be  more 
important  than  the  actual  value  of  the  property  at  any  given 
value  of  the  parameter.  Also,  combinatorial  methodologies 
may  not  be  useful  for  determination  of  properties  that  are 
highly  structure  sensitive.  In  any  case,  in  an  ideal  combinato¬ 
rial  experiment,  the  results  should  be  benchmarked  to  some 
key  data  that  have  been  determined  by  more  traditional 
means;  in  this  review,  we  will  endeavor  to  supply  such  infor¬ 
mation  whenever  possible. 

Combinatorial  materials  science  has  traditionally  been 
focused  on  materials  discovery,  screening,  and  optimization 
to  combat  the  extremely  high  cost  and  long  development 
time  of  new  materials  and  their  introduction  into  commerce. 
Going  forward  this  new  paradigm  may  be  driven  by  other 
needs  such  as  materials  substitution  and  experimental  verifi¬ 
cation  of  materials  properties  predicted  by  modeling  and 
simulation. 

II.  ELECTRONIC  MATERIALS 
A.  Silicon  microelectronics 
1.  Advanced  gate  stack  materials 

Continuous  scaling  of  complementary-metal-oxide- 
semiconductor  (CMOS)  transistors  to  meet  the  evolving 
requirements  of  the  International  Technology  Roadmap  for 
Semiconductors  (ITRS)  has  made  the  traditional  gate  stack, 
Si02  gate  dielectric  and  polycrystalline  Si  (poly-Si)  gate 
electrode,  unsuitable  for  future  integrated  circuit  devices.34 
Further,  as  is  shown  in  Fig.  1,  only  about  a  dozen  elements 
were  involved  in  integrated  circuit  manufacturing  in  1990,  as 
opposed  to  over  fifty  by  the  year  2000.  The  introduction  of 
so  many  new  elements  and  materials  in  the  form  of  their 
alloys  and  compounds  should  have  provided  an  excellent 
opportunity  for  the  widespread  use  of  combinatorial  method¬ 
ologies  in  the  microelectronics  industry.  However,  the 
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FIG.  1 .  Increasing  utilization  of  different 
elements  in  the  silicon  microelectronics 
industry,  as  a  function  of  time. 
Reproduced  with  permission  from  IBM 
Corporation. 


industry  was  slow  in  adopting  such  high  throughput  techni¬ 
ques,  perhaps  due  to  the  high  cost  of  the  suite  of  deposition 
and  testing  equipment  required.35  Nonetheless,  some  excel¬ 
lent  high  throughput  materials  research  has  been  focused  on 
Si  microelectronics. 

a.  High-K  gate  dielectrics.  Since  the  advent  of  the  inte¬ 
grated  circuit,  Fig.  2,  Si02  has  been  the  gate  dielectric  mate¬ 
rial  of  choice  for  silicon  microelectronics.36  However,  since 
high  gate  capacitance  is  necessary  to  produce  the  required 
drive  currents  for  integrated  circuit  devices,  and  further, 
capacitance  is  inversely  proportional  to  gate  dielectric  thick¬ 
ness,  Si02  gate  dielectrics  have  been  scaled  to  ever  thinner 
dimensions.  Because  Si02  has  a  relatively  low  dielectric 
constant,  k ,  of  3.9,  such  scaling  resulted  in  Si02  layers  so 
thin  (^1.0  nm)  that  excessive  leakage  current  compromised 
their  performance.37  Alternate,  “high-K”  gate  dielectrics 
were  therefore  required.38 

The  first  high  throughput  experiments  in  high  dielectric 
constant  electronic  oxides  (< albeit  for  memory  cells,  not  gate 
dielectrics)  were  by  van  Dover  et  al.39  Using  multiple  sput¬ 
tering  guns  to  create  continuous  composition  spread  (CCS) 
libraries,  they  found  compositions  in  the  Zr-Sn-Ti-0  system 
that  displayed  high  figures  of  merit  (based  on  the  product  of 


FIG.  2.  Cross  section  of  an  integrated  circuit  transistor  device,  showing  the 
gate  stack  (gate  metal  electrode  plus  gate  dielectric),  the  source  and  drain 
regions,  and  the  channel  (in  between  the  source  and  drain). 


the  capacitance  per  unit  area  and  the  operating  voltage), 
shown  in  Fig.  3.  They  studied  about  thirty  different  ternary 
oxide  systems,  sometimes  making  four  thousand  measure¬ 
ments  per  day  on  a  single  library. 

High-/c  gate  dielectric  materials  have  been  the  subject  of 
several  combinatorial  studies,  starting  with  an  investigation 
of  the  Hf02-Y203-Al203,40^2  Fig.  4,  and  Zr02-Y203-Al203, 
Fig.  5,  systems.40  Pulsed  laser  deposition  (PLD)  or  ion  beam 
sputtering  was  used  to  deposit  these  CCS  library  films, 
Fig.  6,  and  a  scanning  microwave  microscope,43,44  Fig.  5, 
was  used  to  measure  the  dielectric  properties  of  the  films 
without  having  to  fabricate  capacitors.  In  addition,  the  use  of 
X-ray  diffraction  (XRD)  analysis  enabled  Chikyow  et  al 40 
to  identify  an  amorphous  Hf02-Y203-Al203  composition, 
having  the  ratio  6:1:3,  as  the  most  promising  gate  dielectric 
candidate.  Since  the  dielectric  constant  of  Hf02  (^20)  is  the 
highest  of  the  three  components  in  their  system,  it  is  not  sur¬ 
prising  that  the  amorphous  dielectric  of  choice  is  rich  in  that 
oxide.  XRD  evidence41  indicated  the  presence  of  cubic, 
monoclinic,  and  amorphous  phases  in  the  Hf02-Y203-Al203 
system.  Non-combinatorial  studies  of  Hf02  (Ref.  45)  and  the 
Hf02-Al203  (Ref.  46)  system  are  in  agreement  with  the  basic 
finding  of  cubic  and  monoclinic  phases,  and  the  fact  that 
A1203  addition  stabilizes  the  amorphous  phase.  Whereas  in 
the  combinatorial  study  a  Hf02  to  (Y203  +  A1203)  ratio  of 
6:4  was  sufficient  to  stabilize  the  amorphous  phase  after  a 
700  °C  anneal  for  30  min,  single  composition  Hf02-Al203 
films  with  a  Hf02  to  A1203  ratio  of  1.5:1  was  stable  in  the 
amorphous  phase  up  to  900  °C,  and  films  with  a  Hf02  to 
A1203  ratio  of  1:1.5  remained  amorphous  after  a  1050  °C 
rapid  thermal  anneal.46  This  result  highlights  the  importance 
of  benchmarking  combinatorial  experiments  with  conven¬ 
tional  experimental  approaches. 

Another  recent  high  throughput  study  looked  at  the 
RE203-Al203-Hf02  systems,  where  RE  is  the  rare  earth  ele¬ 
ments,  in  this  case  La,  Ce,  and  Y.47  The  goal  of  the  experi¬ 
ment  was  to  determine  which  compositions  in  these  systems 
might  grow  epitaxially  on  Si  (100)  surfaces,  and  thus  perhaps 
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FIG.  3.  Figure  of  merit  (FOM)  of  ZrxSnyTiz02_^  dielectric  thin  films:  (a)  FOM  as  a  function  of  position  on  the  6.2  x  6.5  cm  substrate,  obtained  by  multiplying 
the  capacitance  and  voltage  measured  at  each  point,  and  normalizing  to  the  capacitor  area.  The  location  of  the  Zr,  Ti,  and  Sn  guns  is  indicated  schematically 
with  respect  to  the  position  of  the  substrate;  (b)  same  data  mapped  onto  a  conventional  ternary  composition  diagram.  Reprinted  by  permission  from  R.  B.  van 
Dover  et  al.,  Nature  392,  162  (1998).  Copyright  1998  Macmillan  Publishers  Ltd. 


FIG.  4.  Ternary  composition  spread  film  grown  on  a  Si(100)  substrate.  The 
equilateral  triangle  in  the  center  corresponds  to  the  ternary  phase  diagram 
and  the  three  isosceles  regions  around  it  correspond  to  three  binary  phase 
diagrams.  Reproduced  with  permission  from  K.  Hasegawa  et  al.,  Appl.  Surf. 
Sci.  223,  229  (2004).  Copyright  2004  Elsevier  BY. 


mimic  the  low  defect  density,  continuous  random  network 
structure  of  the  native  Si02/Si  interface.  They  discovered  a 
new  ionic  structure  that  accommodated  local  atomic  altera¬ 
tions  as  easily  as  native  Si02  and  could  serve  as  a  model  for 
advanced  high-K  gate  dielectrics.  One  desirable  composition, 
although  exhibiting  high  leakage  current  most  likely  due  to 
its  small  bad  gap  offset,  was  a  22  nm  film  of  the  composition 
Ce6(A103)5.1(Si03.5)o.9.  In  the  course  of  performing  this 
work,  the  authors  introduced  a  method  for  high  throughput 
structural  characterization  of  CCS  libraries  using  high 
energy  (synchrotron)  X-ray  beams  in  transmission  geome¬ 
try.48  Using  lOOnm  thick  films  in  the  Ce203-Al203-Hf02 
system,  they  were  able  to  obtain  texture,  epitaxial  relation¬ 
ship,  and  phase  information  in  their  combinatorial  library 
film,  Fig.  7.  They  were  able  to  identify  fifteen  different 
phases  because  the  films  were  relatively  thick  and  they  used 
an  intense  (synchrotron)  beam.  In  general,  laboratory  XRD 
analyses  of  such  combinatorial  library  films  result  in  the 
identification  of  only  the  most  predominant  phases,  such  as 
has  been  observed  for  the  TaN-AIN  (gate  metal  electrode) 

49 

system. 

PLD  has  been  the  film  growth  technique  of  choice  for 
several  gate  dielectric  combinatorial  experiments.  Although 
PLD  is  not  a  commercially  viable  deposition  technique 
because  the  process  often  involves  the  creation  of  myriad 
particles,  it  is  a  very  valuable  laboratory  technique  since  it 


FIG.  5.  Schematic  illustration  of  the 
microwave  microscope  for  dielectric 
property  mapping,  and  a  measured  exam¬ 
ple  for  a  (Zr02)x(Y203)y(Al203)z  combi¬ 
natorial  library  sample.  Reproduced 
with  permission  from  T.  Chikyow  et  al., 
IEEE  Cat.  No.  06EX1294  (2006). 
Copyright  2006  IEEE. 
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FIG.  6.  (a)  Schematic  illustration  of  the 
ternary  and  binary  composition  spread 
method  with  sample  rotation  and  moving 
mask  system,  first  employed  by 
Fukumura  et  al.104  Here,  the  PLD  depo¬ 
sition  method  is  shown  as  an  example, 
but  ion  beam  sputtering  can  also  be 
used.  Reproduced  with  permission  from 
T.  Chikyow  et  al.,  IEEE  Cat.  No. 
06EX1294  (2006).  Copyright  2006 
IEEE,  (b)  SCHEMATIC  showing  the 
synthesis  of  a  combinatorial  “library”  of 
the  binary  system  A-B.  A  thin  wedge 
(wedge  height  ~0.5  nm)  of  A  is  depos¬ 
ited,  and  the  substrate  is  rotated  180°.  A 
complementary  wedge  of  B  is  then  de¬ 
posited.  This  is  repeated  “N”  times,  until 
the  desired  film  thickness  is  achieved.  If 
the  substrate  is  rotated  120°  instead,  a 
ternary  (shown  from  above  in  the  lower 
left),  rather  than  a  binary  alloy  system 
may  be  synthesized.  Reproduced  with 
permission  from  M.  L.  Green  et  al.. 
Microelectron.  Eng.  84,  2209  (2007). 
Copyright  2007  Elsevier  B. 


(b) 


results  in  nearly  perfect  stoichiometry  transfer  between  the 
target(s)  and  the  deposited  film.  Recent  studies  were  aimed 
at  optimizing  the  compositional  coverage  of  PLD  films 
while  maintaining  a  uniform  thickness  across  the  sub¬ 
strate.50,51  In  these  studies,  Al203-Hf02-Y203  combinato¬ 
rial  libraries  were  deposited  by  the  “natural”  CCS 
technique,  i.e.,  films  were  sequentially  deposited  at  various 
(x,  y)  locations  on  the  substrate  wafer,  from  different  tar¬ 
gets;  no  shutters  or  masks  were  used,  unlike  the  case  illus¬ 
trated  in  Fig.  6.  The  composition  and  thickness  of  such  films 
can  be  tuned  by  adjusting  the  laser  spot  size  and  fluence, 
background  gas  pressure,  and  target  geometry.  Using  a  sim¬ 
ple  bimodal  cosine  power  distribution  function  for  the  depo¬ 
sition  rate  from  the  three  individual  targets,  the  composition 
and  thickness  of  ternary  library  films  could  be  predicted. 
The  experimentally  determined  film  thickness  (by  spectro¬ 
scopic  reflectometry)  and  composition  (by  electron  probe 
microanalysis  wavelength  dispersive  spectroscopy  (EPMA- 
WDS)  results  agreed  with  the  predictions  to  within  measure¬ 
ment  resolution,  as  is  shown  in  Figs.  8  and  9.  This  modeling 
approach  results  in  a  high  degree  of  confidence  that  one  can 
predict  thickness  and  compositional  coverage,  making  the 
combinatorial  PLD  tool  very  useful  for  routine  library 
design,  while  avoiding  costly  and  time  consuming  experi¬ 
mental  verification  of  library  composition. 


In  another  PLD  study  of  gate  dielectrics,52  combinato¬ 
rial  library  films  of  the  Hf02-Ti02-  Y203  system  were 
grown.  Ti02  has  a  relatively  high  dielectric  constant,  k  ~  80, 
but  is  not  thermodynamically  stable  on  Si.53  Thus,  the  goal 
of  this  experiment  was  to  screen  for  a  gate  dielectric  compo¬ 
sition  that  had  a  higher  dielectric  constant  than  that  of  Hf02 
(about  20),  was  thermally  stable  on  Si,  and  also  be  stable  as 
an  amorphous  phase.  Interestingly,  it  was  found  that  films  in 
the  Hf02-Ti02-Y203  library  exhibited  a  visible  boundary 
separating  amorphous  and  crystalline  phases,  Fig.  10,  and 
that  this  boundary  could  be  manipulated  by  varying  substrate 
temperature  or  Ti02  distribution.  Further,  the  authors  were 
able  to  eliminate  the  visual  boundary,  and  therefore  realize 
two  microstructural  extremes;  the  entire  Hf02-Ti02-Y203 
ternary  film  could  be  amorphous  after  deposition  at  300  °C, 
or  crystalline  in  the  Hf02-Y203  binary  system.  Electrical 
measurements  demonstrated,  Fig.  11,  that  the  dielectric  con¬ 
stant  of  the  library  films  could  be  varied  between  about  10 
and  120.  Surprisingly,  some  of  the  highest  dielectric  con¬ 
stants  are  observed  near  the  Ti02-Y203  binary  system,54 
containing  only  small  fractions  of  Hf02. 

Another  recent  study  focused  on  the  Y203-A1203  sys¬ 
tem,  with  the  goal  of  forming  high  dielectric  constant,  ther¬ 
mally  stable,  yttrium  aluminate  gate  dielectrics.55  The 
library  films  were  grown  by  PLD  using  discrete  masks,  and 
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FIG.  7.  Bubble  plot  reconstruction  of  the  crystal  phase  diagram  for  the 
Ce02-Al203-Hf02  CCS  thin  him  on  a  Si(001)substrate.  The  various  phases 
are  numbered.  The  lower  figure  is  an  HEXRD  diffraction  image  taken  from 
the  polycrystalline  Ce02  part  of  the  composition  spread.  Reprinted  with  per¬ 
mission  from  D.  A.  Kukuruznyak  et  al.,  Appl.  Phys.  Lett.  91,  071916 
(2007).  Copyright  2007  American  Institute  of  Physics. 


the  effects  of  Al  doping  and  annealing  on  the  thermal  stabil¬ 
ity  of  the  Y203/Si  gate  stack  were  studied  by  X-ray 
photoemission  spectroscopy  (XPS)  and  X-ray  absorption 
spectroscopy.  Although  Al  doping  was  effective  in  increas¬ 
ing  the  band  gap  of  Y203,  Fig.  12,  it  did  not  suppress  yttrium 
silicate  formation  via  diffusion  of  Si02  from  the  interfacial 
layer  into  the  yttrium  aluminate;  annealing  promoted  the  for¬ 
mation  of  silicate  compounds  in  the  film  via  mixing  of  the 
silicate  layer  and  the  Si  substrate.  The  observed  band  gap 
bowing  in  the  Y203-A1203  system  was  attributed  to  para¬ 
bolic  changes  in  the  conduction  band  offset  as  a  function  of 
Al  concentration. 

The  development  of  high-throughput  chemical  vapor 
deposition  (CVD)  methodologies  for  high-/c  gate 
dielectrics56-58  is  an  important  effort,  since  CVD  and  a 
related  technology,  atomic  layer  deposition  (ALD)59  may  be 
the  manufacturing  methods  of  choice  for  these  critical  films. 
Thus,  it  is  important  to  explore  the  feasibility  of  performing 
combinatorial  experimentation  in  such  a  platform.  While  it  is 


non-trivial  to  generate  large  composition  variation  in  a  CVD 
chamber,  they  can  be  achieved  by  controlling  the  flow  of 
precursor  gases.60,61  Equally  important  to  achieving  compo¬ 
sitional  variations  on  the  wafers  is  the  need  to  achieve  pro¬ 
cess  parameter  (i.e.,  pressure  and  flow)  variations. 

b.  Metal  gate  electrodes.  Substitutions  for  the  poly-Si 
gate  electrode,  the  second  component  of  the  CMOS  gate 
stack  of  Fig.  2,  using  true  metals  or  metalloids  that  would 
provide  higher  carrier  densities  and  therefore  lower  resistiv¬ 
ities  as  well  as  less  carrier  depletion,  were  not  developed  until 
after  the  high-K  gate  dielectrics.  Therefore,  the  selection  of 
metal  gate  substitutes  for  poly-Si,  which  must  be  compatible 
with  the  new  high-K  gate  dielectrics,  is  not  as  advanced.  For 
high  performance  CMOS  applications,  the  metal  gate  electro¬ 
des  must  have  work  functions  (Om)  that  are  aligned  with  the 
conduction  (4.05  ±  0.2  eV)  and  valence  (5.17  ±  0.2 eV) 
bands  of  Si  for  n-channel  metal-oxide-semiconductor 
(NMOS)  and  p-channel  metal-oxide-semiconductor  (PMOS) 
applications,  respectively.34 

In  one  of  the  first  high  throughput  studies  of  metal  gate 
electrodes,  Ni-Ti-Pt  ternary  libraries  were  deposited  by 
magnetron  sputtering  through  moving  masks  on  Hf02  dielec¬ 
trics,  and  studied  systematically  by  WDS,  micro-XRD,  and 
capacitance-voltage  (C-V)  and  current-voltage  (I-V)  analy¬ 
ses.62  WDS  results  showed  that  compositional  coverages  up 
to  about  90%  for  Ni  and  Ti,  and  75%  for  Pt  were  achieved  in 
the  library  film,  as  can  be  seen  in  Fig.  13.  A  more  negative 
flatband  voltage  (Vfb)  was  observed  close  to  the  Ti-rich  cor¬ 
ner  than  to  the  Ni-  and  Pt-rich  comers,  implying  a  smaller 
®m  near  the  Ti-rich  corner,  and  higher  ®m  near  the  Ni-  and 
Pt-rich  corners.  Not  only  is  that  observed,  as  can  be  seen  in 
Fig.  13(b),  but  also  the  ®m  values  measured  for  the  nearly 
pure  metals  (i.e.,  near  the  comers)  in  this  combinatorial 
experiment  are  in  very  good  agreement  with  bulk  literature 
values  for  Ni,  Ti,  and  Pt.  Thus,  combinatorial  methodologies 
proved  to  be  accurate  and  useful  in  surveying  the  large  com¬ 
positional  space  of  ternary  alloy  metal  gate  electrode 
systems. 

When  using  a  moving  shutter  for  library  film 
deposition63-65  to  deposit  a  wedge  of  material  at  a  time,  as 
shown  in  Fig.  6,  the  shape  of  the  wedge,  i.e.,  amount  depos¬ 
ited  as  a  function  of  distance  in  the  plane  of  the  wafer,  must 
be  calibrated  very  carefully  if  one  is  to  have  a  sufficiently 
accurate  idea  of  the  compositional  spread.  The  wedge  height 
is  very  often  small  (a  few  tenths  of  a  nm)  to  enable,  even  in 
the  as-deposited  condition,  atomic  mixing  in  the  direction 
perpendicular  to  him  growth.  Figure  14  illustrates  the  use  of 
Rutherford  backscattering  (RBS)  as  a  calibration  procedure 
for  a  wedge  of  TaN  deposited  by  sputtering  with  the  use  of  a 
moving  shutter.49,65  RBS  is  of  course  very  helpful  for  cali¬ 
brating  deposition  profiles  using  any  deposition  method. 

High  throughput  experiments  were  also  performed  on 
the  Mo-Ru66  and  Pt-W67  metal  gate  electrode  systems  depos¬ 
ited  on  high-/c  dielectrics.  As  can  be  seen  in  Fig.  15,  the  addi¬ 
tion  of  10%-30%  Mo  in  Ru  electrodes,  deposited  on  Hf-Si- 
O-N,  enables  tuning  of  the  values,  likely  due  to  Mo  seg¬ 
regation  at  the  metal  gate/high- k  interface.  Further,  carbon 
incorporation  was  found  to  be  effective  in  refining  the  metal 


Downloaded  20  Aug  2013  to  131.215.127.215.  This  article  is  copyrighted  as  indicated  in  the  abstract.  Reuse  of  AIP  content  is  subject  to  the  terms  at:  http://jap.aip.org/about/rights_and_permissions 


231 1 01  -7  Green,  Takeuchi,  and  Hattrick-Simpers 


J.  Appl.  Phys.  113,  231101  (2013) 


FIG.  8.  Comparison  of  simulated  ((9a) 
and  (9b))  and  experimentally  determined 
((9c)  and  (9d))  composition  maps  for 
PLD  libraries  in  the  Hf-Al-YOx  system. 
The  library  in  (9a)  was  deposited  at 
P02=  1.33  Pa,  and  the  library  in  (9b)  at 
Po2=  13.3  Pa.  Reproduced  with  permis¬ 
sion  from  P.  K.  Schenck  et  al.,  Appl. 
Surf.  Sci.  254,  781  (2007).  Copyright 
2004  Elsevier  BY. 


FIG.  9.  (a)  Photograph  and  (b)  predicted 
thickness  map  for  a  PLD  Hf-Al-YOx 
library  him  deposited  with  P02  =  13.3  Pa. 
Reprinted  with  permission  from  N.  D. 
Bassim  et  al.,  Rev.  Sci.  Instrum.  78, 
072203  (2007).  Copyright  2008  American 
Institute  of  Physics. 
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FIG.  10.  Photographs  of  four  Hf02- 
Ti02-Y203  library  films,  each  deposited 
at  a  different  temperature.  Those  library 
films  deposited  at  T>400°C  show  a 
clear  boundary  separating  the  Ti02-rich 
region  from  the  Hf02-Y203  region. 
Reprinted  with  permission  from  J.  L. 
Klamo  et  al.,  J.  Appl.  Phys.  107,  054101 
(2010).  Copyright  2008  American 
Institute  of  Physics. 
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FIG.  11.  Ternary  diagram  showing  the  predicted  compositional  coverage 
(small  points)  of  the  500  °C  Hf02-Ti02-Y203  library  film,  and  the  meas¬ 
ured  relative  dielectric  constant  (filled  colored  points)  extracted  from  the 
C-V  characteristics  of  the  as-deposited  film.  The  overlaid  line  represents  the 
approximate  location  of  the  crystalline-amorphous  boundary.  Reprinted  with 
permission  from  J.  L.  Klamo  et  al.,  J.  Appl.  Phys.  107,  054101  (2010). 
Copyright  2010  American  Institute  of  Physics. 

gate  electrode  grain  size,  thereby  adding  another  Vrb  control 
mechanism.  Similarly,  a  wide  control  range  (~1.2  V)  of  Vfb 
was  observed  for  (Pt-W)/La203  gate  stacks.  Again,  the  effect 
was  attributed  to  charge  transfer  at  the  metal  gate/La203 
interface,  which  induced  an  increase  in  the  effective  work 
function,  especially  at  high  Pt  composition  ratios,  as  can  be 
seen  in  Fig.  16.  The  control  effect  was  much  stronger  for 
La203  than  for  Hf02  gate  dielectrics. 

Metalloids  such  as  metal  nitrides  and  carbides  have 
attracted  attention  as  potential  gate  metal  electrodes  due  to 
their  superior  thermal  and  chemical  resistance  compared  to 
metals.  Among  the  many  metal  gate  candidates,  TaN  has 
shown  very  useful  ®m  tunability  through  alloying  with  vari¬ 
ous  elements.68-71  In  particular,  Ta!_xAlxNy  alloys  are  easily 
deposited,  with  very  good  electrical  and  chemical  proper¬ 
ties.70  The  Ta!_xAlxNy  system  is  an  ideal  combinatorial 
materials  science  problem  because  systematic  measurement 
of  ®m  across  the  wide  composition  range,  x,  is  not  trivial, 
since  capacitor  fabrication  and  characterization  based  on  a 
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Ti:  4.37  eV  [17] 

FIG.  13.  (a)  Map  of  the  flatband  voltage  shift  extracted  from  the  C-V  char¬ 
acteristics  at  1  kHz  across  the  Ni-Ti-Pt  library,  and  (b)  map  of  the  extracted 
work  functions  across  the  library.  Round  symbols  in  the  three  corners  denote 
literature  values.  Reprinted  with  permission  from  K.  S.  Chang  et  al.,  Appl. 
Phys.  Lett.  89,  142108  (2006).  Copyright  2008  American  Institute  of 
Physics. 

“one-composition-at-a-time”  approach  are  extremely  time 
consuming;  thus  very  little  data  are  available  for  this  metal 
gate  alloy  system.72  Recently,  a  combinatorial  approach  was 
applied  to  ®m  extraction  for  the  metal  gate  electrode  system 
Ta!_xAlxNy  deposited  on  Hf02.49  Over  two  thousand  capaci¬ 
tors  on  four  identical  metal  gate  libraries  (< albeit  with  differ¬ 
ing  capacitance  due  to  Si02  underlayers  of  different 
thicknesses)  were  automatically  measured,  from  which 
shifts  were  extracted  and  ®m  determined.  Figure  17  shows 
the  array  of  test  capacitors  in  the  library  sample;  although  a 
composition  range  of  0.05  <  x  <  0.85  was  achieved  on  the 
library,  electrical  measurements  were  only  performed  on 
samples  with  0.05  <  x  <  0.50,  because  higher  Al  concentra¬ 
tions  were  not  suitably  conductive.  As  can  be  seen  from  Fig. 
18,  ®m  of  the  Ta!_xAlxNy/Hf02/Si02/Si  gate  stack  structures 
can  be  tuned  as  a  function  of  gate  metal  composition  over  a 
wide  composition  range,  as  well  as  by  annealing.  From  these 
results,  it  was  suggested  that  Ta0.9AlaiNL24  gate  metal  elec¬ 
trodes  may  be  useful  for  PMOS  applications.  Figure  18  is  a 
noteworthy  example  of  the  power  of  the  combinatorial  tech¬ 
nique.  In  this  graph,  ®m  is  shown  as  a  function  of  gate  metal 
composition,  with  about  thirty  compositions  represented,  as 


FIG.  12.  Al  concentration  dependence 
of  (a)  valence-band  offset  and  (b)  band 
gap  by  annealing  samples  at  400,  600, 
and  800  °C.  The  dashed  lines  indicate 
(a)  linear  and  (b)  parabolic  fitting  results. 
Reprinted  with  permission  from  S. 
Toyoda  et  al.,  J.  Vac.  Sci.  Technol.  A 
28(1),  16-19  (2010).  Copyright  2010 
American  Institute  of  Physics. 
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FIG.  14.  Compositional  calibration  of  a  TaN-AIN  library  sample  achieved 
by  performing  RBS  analysis  on  a  wedge  profiled,  moving  shutter-deposited 
TaN  (for  example)  film.  Reproduced  with  permission  from  M.  L.  Green 
et  al.,  Microelectron.  Eng.  84,  2209  (2007).  Copyright  2007  Elsevier  BY. 
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FIG.  15.  At  CRu  =  0%,  30%,  and  100%,  the  values  in  Vfb  are  typical  of 
metal/high-K;  interfaces,  i.e.,  much  narrower  Vfb  controllability  on  high -k 
films  than  work  function  values.  However,  at  CRu  =  30%,  a  wide  controll¬ 
ability  of  0.58  V  has  been  achieved  on  HfiSiON  after  500  °C  due  to  Mo  seg¬ 
regation  at  the  interface.  Modification  of  the  micro  structure  due  to  carbon 
incorporation  plays  a  key  role.  Reproduced  with  permission  from  2007 
IEEE  Electron  Devices  Meeting  -  IEDM  2007,  345.  Copyright  2007  IEEE. 
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FIG.  17.  (a)  Composition  map  of  the  Taf_xAlxNy  composition  spread,  as 
determined  by  WDS.  The  dimension  of  the  library  is  15  mm  x  15  mm.  Each 
dot  represents  a  capacitor,  and  (b)  plot  of  the  composition  variation  across 
the  sample.  A  wide  composition  range,  x  =  0.05  to  0.85,  was  achieved. 
Reproduced  with  permission  from  K.  S.  Chang  et  al.,  IEEE  Trans.  Electron 
Devices  55,  2641  (2008).  Copyright  2008  IEEE. 


opposed  to  isolated  data  for  only  a  few  of  those  composition 
points  found  in  the  literature.  Further,  the  data  instantly 
reveal  two  other  phenomena.  First,  one  can  see  the  consistent 
increase  in  ®m  after  anneals  at  higher  temperatures  than  the 
initial  forming  gas  anneal  (500  °C),  and  second,  the  thermal 
instability  after  the  1000  °C  anneal  of  gate  electrodes  con¬ 
taining  less  than  about  27%  Al.  Consistent  and  comprehen¬ 
sive  sets  of  working  data  such  in  Fig.  18  are  almost  non¬ 
existent  in  the  semiconductor  microelectronics  literature. 

Another  metalloid  system,  Ta-C-N,  has  attracted  atten¬ 
tion  as  a  potential  metal  gate  electrode,  also  due  to  the  possi¬ 
bility  of  ®m  tuning,  high  temperature,  chemical  stability,  and 
good  mechanical  properties.  Combinatorial  Ta-C-N  library 
films  were  deposited  by  magnetron  sputtering  on  Hf02,73 
and  WDS  and  fluorescence  yield  near  edge  spectroscopy 
(FYNES)  were  used  to  quantitatively  determine  the 
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FIG.  16.  Slope  factor  S  (derivative  of  Vfb  as  a  function  of  work  function),  as 
a  function  of  Pt-W  composition  after  a  450  °C  forming  gas  anneal  followed 
by  a  350  °C  oxidizing  gas  anneal.  Note  that  the  S  factor  of  La203  overlaps 
that  of  Hf02  except  for  an  anomalous  increase  in  the  range  0.7<Pt<  1. 
Reproduced  with  permission  from  T.  Chikyow  et  al.,  IEEE  Cat.  No. 
06EX1294  (2006).  Copyright  2006  IEEE. 


FIG.  18.  Plot  of  the  extracted  work  functions  (<Dm)  for  the  Tai_xAlxNy  com¬ 
position  spreads.  <Dm  was  extracted  as  a  function  of  Al  content  after  a  form¬ 
ing  gas  anneal  (500 °C),  and  900  °C  and  1000  °C  rapid  thermal  anneals. 
After  the  1000  °C  anneal,  for  x  <  0.27,  ®m  could  not  be  mapped  due  to  deg¬ 
radation  of  the  capacitors.  Reproduced  with  permission  from  K.  S.  Chang 
et  al.,  IEEE  Trans.  Electron  Devices  55,  2641  (2008).  Copyright  2008  IEEE. 
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composition  across  these  films.  Interestingly,  scanning 
micro-XRD  determined  that  a  solid  solution  of  Ta(C,  N)x 
forms  and  extends  to  compositions  that  were  previously 
unknown  (0.3  Ta  <  0.5  and  0.57  <  Ta  <  0.67),  as  is  shown  in 
Fig.  19.  In  addition,  the  thermal  stability  of  the  Ta-C-N/ 
Hf02  library  film  was  studied  using  high  resolution  transmis¬ 
sion  electron  microscopy  (TEM),  which  showed  that  the 
Ta-C-N/Hf02/Si02/Si  gate  stack  exhibited  good  thermal 
stability  up  to  950  °C.74 

III.  FERROELECTRIC,  PIEZOELECTRIC,  AND 
MULTIFERROIC  MATERIALS 

Ferroelectrics,  which  exhibit  a  spontaneous  electric 
polarization  that  can  be  reversed  by  the  application  of  an 
external  electric  field,  are  an  important  class  of  functional 
materials,  with  a  variety  of  applications.75,76  In  thin  film 
form,  they  are  important  for  non-volatile  memory  devices. 
Recently,  their  applications  for  tunable  microwave  devices 
have  been  explored.  Depending  on  the  intended  application, 
various  properties  such  as  dielectric  constant  and  its  tunabil- 
ity,  dielectric  loss,  spontaneous  polarization,  coercive  field, 
and  breakdown  voltage  must  be  optimized.  Many  ferroelec¬ 
tric  materials  are  also  good  piezoelectrics,  which  are  ubiqui¬ 
tous  in  a  wide  range  of  sensor  and  actuator  applications.77,78 
One  important  research  direction  in  the  piezoelectrics  com¬ 
munity  is  the  identification  of  lead-free,  high  piezoelectric 
coefficient  replacements  for  the  commonly  used  compound 
Pb(Zr,Ti)03  (PZT). 

Ferroelectric  materials  often  have  a  perovskite  or  related 
crystal  structure.  Thus,  the  emphasis  of  some  of  the  first 
combinatorial  studies  was  on  synthesis  of  perovskite  struc¬ 
tures  in  thin  film  libraries  with  a  large  compositional  varia¬ 
tion.  Chang  et  al .79  employed  the  precursor  approach  to 
fabricate  (Ba,Sr)Ti03  thin  films  and  systematically  investi¬ 
gate  the  effect  of  transition  metal  (TM)  dopants  on  their 
dielectric  constants  and  losses.  Different  thicknesses  of 
amorphous  BaF2,  SrF2,  Ti02  layers  were  deposited  to  control 
the  desired  final  A-site  composition  of  the  (Ba,Sr)Ti03  thin 
film  at  each  site  on  the  library.  A  thermal  process  involving  a 
series  of  annealing  steps  to  interdiffuse  the  precursor  layers 


(and  dopants),  remove  F2,  and  crystallize  the  titanate  com¬ 
pounds  was  developed.  Cross-sectional  TEM  revealed  that 
quasi-epitaxial  single  crystalline  thin  films  can  be  fabricated 
using  the  precursor  approach  on  lattice  matched  substrates. 
A  modified  precursor  approach  has  also  been  developed  for 
fabricating  compounds  containing  high  vapor  pressure  ele¬ 
ments  such  as  Pb.80 

Figure  20,  from  Chang  et  al.?9  shows  the  composition 
distribution  of  a  doped  (Ba,Sr)Ti03  thin  film  library.  This 
was  also  the  first  study  where  quantitative  scanning  micro- 
wave  microscopy  (SMM)81  was  used  as  an  effective  tool  to 
scan  and  screen  dielectric  properties  of  thin  film  libraries.  An 
important  finding  from  this  study  was  that  high  valence  TM 
dopants  such  as  W  are  effective  in  reducing  the  dielectric 
loss  by  minimizing  the  generation  of  oxygen  vacancies  dur¬ 
ing  the  crystallization  process.  A  similarly  designed  library 
was  also  reported  in  this  paper  to  study  the  effect  of  dopants 
on  the  properties  of  (Ba,Sr)Ti03  high-K  capacitors  with  elec¬ 
trodes.  It  was  found  that  the  same  high  valence  TM  ions  are 
instrumental  in  reducing  the  leakage  current  and  boosting 
the  breakdown  voltage  of  the  capacitors.  Further  research  on 
single  composition  samples  showed  that  W  doping  of 
(Ba,Sr)Ti03  thin  films  indeed  results  in  reduced  dielectric 
loss.82 

A  recent  surge  of  research  activity  in  tunable  microwave 
devices  for  phased  array  antenna  and  other  communication 
devices  has  brought  much  attention  to  the  microwave  proper¬ 
ties  of  ferroelectric  thin  films  in  general.  The  100  MHz  to 
100  GHz  (approximately)  frequency  band  is  a  crucial  range 
for  these  materials,  where  domain  wall  motion  and  dielectric 
relaxation  associated  with  the  soft  phonon  mode  are  known 
to  occur.83  Thus,  it  is  important  to  understand  the  frequency 
dependence  of  their  dielectric  (real  and  imaginary)  constants. 
In  one  study,  multifrequency  mode  variable  temperature 
SMM  was  used  to  map  the  frequency  (1-10  GHz),  composi¬ 
tion,  and  temperature  dependence  of  Ba!_xSrxTi03  near  the 
Curie  temperature.84  The  epitaxial  Bai_xSrxTi03  composition 
spreads  were  fabricated  using  PFD  and  the  layer-by-layer 
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FIG.  19.  Ta-C-N  library  film  composition  data  plotted  in  ternary  composi¬ 
tion  space  (the  solid  and  dashed  black  line).  Ta(C,N)x  and  Ta2(C,N)x  equi¬ 
librium  phase  regions  are  also  shown.  Reprinted  with  permission  from  K.  S. 
Chang  et  al.,  Appl.  Phys.  Lett.  96,  192114  (2010).  Copyright  2010 
American  Institute  of  Physics. 


FIG.  20.  Dielectric  loss  mapping  of  doped  BaTi03  films  measured  by  scan¬ 
ning  microwave  microscopy.  The  region  defined  by  the  red  box  are  compo¬ 
sitions  doped  with  W.  The  loss  is  low  in  general  in  this  region.  Reprinted 
with  permission  from  H.  Chang  et  al.,  Appl.  Phys.  Lett.  72,  2185  (1998). 
Copyright  1998  American  Institute  of  Physics. 
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gradient  wedge  approach  similar  to  that  depicted  in  Fig.  6. 
Figure  21,  from  Chang  et  a/.,85  shows  the  frequency  disper¬ 
sion  of  the  dielectric  constant,  as  measured  by  SMM.  The 
fact  that  the  broad  peak  is  observed  in  the  range  x  =  0. 2-0.4, 
where  the  Curie  temperature  is  closest  to  room  temperature 
(RT),  indicates  that  there  is  coupling  between  the  microwave 
field  and  with  the  local  soft  mode  mediated  by  the  presence 
of  structural  defects  such  as  threading  dislocations. 

Other  combinatorial  investigations  of  ferroelectric  mate¬ 
rials  focused  on  chemical  synthesis  techniques  for  preparing 
libraries,86’87  reflecting  the  versatility  of  high  throughput 
methodologies. 

In  searching  for  new  materials  with  high  piezoelectric 
coefficients,  a  particular  emphasis  is  placed  on  identifying 
compositions  at  structural  boundaries  called  morphotropic 
phase  boundaries  (MPBs).  At  MPBs,  materials  undergo 
structural  phase  transitions,  and  the  MPB  is  vertical  in  a 
temperature-composition  phase  diagram,  i.e.,  it  is  isocompo- 
sitional.  At  such  compositions,  small  perturbations  in  the 
material  caused  by  electric  field  or  stress  can  lead  to  an 
enhanced  piezoelectric  response.  Thus,  the  most  commonly 
used  piezoelectric  materials  are  at  MPB  compositions  such 
as  PbZro.52Tio.4sO3  and  (PbMg1/3Nb2/303)o.7  -(PbTiO3)0.3.  In 
a  recent  CCS  library  study,  BiFe03,  an  intrinsic  multiferroic 
material  that  displays  rich  and  complex  physical  properties 
was  selected  as  a  starting  compound  to  screen  for  a  novel 
lead  free  piezoelectric  material,  specifically  at  MPB 
compositions.88-90  A  series  of  Bi!_x(RE)xFe03  pseudobi¬ 
nary  CCS  epitaxial  libraries  (200  nm  thick)  were  fabricated 
epitaxially  at  600  °C  on  SrTi03  (STO)  substrates  by  layer- 
by-layer  PLD  (see,  e.g.,  Fig.  6). 

Quantitative  piezoresponse  force  microscopy  (PFM) 
was  performed  to  measure  the  out-of-plane  piezoresponse. 
The  PFM  scan  of  a  Bi!_xSmxFe03  CCS  library  indicated 
substantial  enhancement  of  the  piezoresponse  near  the  com¬ 
position  Bi_0.8-o.85Sm~o.2-o.i5Fe03,  where  the  structure 
was  found  by  XRD  to  undergo  significant  change  in  the  out 
of  plane  lattice  constant.  The  effective  piezoelectric  coeffi¬ 
cient  values,  d33,  were  then  measured  across  the  composition 
spread,89  Fig.  22.  At  x  =  0.13-0.15,  d33  displays  a  rapid 
increase,  peaking  at  x  =  0.14;  beyond  this  composition,  d33 


BaTi03  Ba^Sr^TiOj  SrTiOj 

FIG.  21.  Normalized  dielectric  constant  dispersion  vs.  composition  on  a 
library  sample,  as  measured  by  microwave  microscope.  The  dispersion  is 
defined  as  (dielectric  constant  at  0.95  GHz — dielectric  constant  at  4.95  GHz)/ 
dielectric  constant  at  0.95  GHz.  The  dispersion  is  largest  in  the  range 
x  =  0.2-0.4.  Reprinted  with  permission  from  K.  S.  Chang  et  al.,  Appl.  Phys. 
Lett.  79,  4411  (2001).  Copyright  2001  American  Institute  of  Physics. 


FIG.  22.  Discovery  of  a  lead  free  morphotropic  phase  boundary.  The  piezo¬ 
electric  coefficient,  d33,  is  mapped  on  a  Bi1_xSmxFe03  composition  spread. 
A  strong  peak  is  observed  at  x  =  0.14.  Reprinted  with  permission  from  S. 
Fujino  et  al.,  Appl.  Phys.  Lett.  92,  202904  (2008).  Copyright  2008  American 
Institute  of  Physics. 

rapidly  decreases.  The  measured  remnant  polarization,  and 
high  field  d33  in  this  case  are  similar  to  values  previously 
reported  for  epitaxial  thin  films  of  the  lead  based  MPB  com¬ 
pounds  mentioned  above.  Detailed  microstructural  studies 
have  shown  that  the  structural  transition  for  this  MPB  is  a 
rhombohedral  to  orthorhombic  structural  transition.  Planar 
TEM  of  a  single  composition  sample  at  x  =  0.14  showed  the 
presence  of  nanoscale  twinned  triclinic  domains  20-50  nm  in 
size.  The  occurrence  of  such  nano  twins  and  the  stress 
accommodation  associated  with  them  is  a  signature  of  an 
adaptive  ferroelectric  phase  at  the  MPB,  which  enables  the 
high  piezoelectric  coefficient.91  Subsequent  work  has 
revealed  that  a  chemical  pressure  effect,  due  to  RE  substitu¬ 
tion  into  the  Bi  site  (perovskite  A-site),  induces  the  structural 
transition  and  the  concomitant  enhancement  in  piezoelectric 

.  99 

properties. 

In  another  study,  the  B-site  substitution  effect  on  the 
ferroelectric  properties  in  (Bi,Sm)Fe03  was  investigated  to 
systematically  investigate  the  relative  effects  of  A-site  and 
B-site  doping,  using  a  pseudotemary  composition  spread.  Sc 
was  selected  as  the  substitutional  dopant  for  the  B-site.93  The 
Sc  ion  has  a  3+  valence,  which  is  not  likely  to  cause  an 
increase  in  leakage  current.  To  make  this  library,  the  moving 
shutter  technique  (e.g.,  Fig.  6)  was  implemented  in  an  epitax¬ 
ial  PLD  process.  On  the  triangularly  shaped  library,  Fig. 
23(a),  the  concentration  of  each  dopant  (A-  or  B-site  dopant) 
varies  along  one  side  of  the  triangle.  The  three  vertex  compo¬ 
sitions  are  BiFe03,  Bi(Fe0.8Sc0.2)O3,  and  (Bi0.8Sm0.2)FeO3.  In 
Fig.  23(a),  all  RT  hysteresis  loops  taken  across  the  spread  are 
displayed.  Ferroelectric  properties  exhibit  clear  changes  in 
response  to  the  substitution  induced  structural  evolution.  In 
the  region  near  the  BiFe03  vertex,  where  the  dopant  composi¬ 
tion  is  relatively  low,  unclosed  loops  are  observed,  indicative 
of  high  leakage  current.  As  the  substitution  of  Sm  and  Sc  is 
increased,  the  square  shaped  ferroelectric  hysteresis  loop  with 
a  saturation  polarization  of  about  70  /iC/cm2  begins  to  appear, 
indicating  that  both  substituted  species  lead  to  a  reduction 
of  the  leakage  current;  an  antiferroelectric  phase  may  be 
playing  a  role  in  suppressing  the  leakage  current.  Comparison 
of  Figs.  23(a)  and  23(b)  shows  that  the  occurrence  of  the 
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FIG.  23.  (a)  Ferroelectric  and  structural  properties  of  a  BiFe03- 
Bi(Fe0.8Sco.2)03-(Bio.8Sm0.2)Fe03  pseudoternary  composition  spread. 
Ferroelectric  hysteresis  loops  measured  at  each  spot  are  plotted.  The  red  line 
marks  the  MPB,  where  the  square  hysteresis  transitions  into  double  hystere¬ 
sis  loops.  The  dotted  line  is  where  the  hysteresis  loops  become  square 
shaped,  indicating  that  the  films  are  non-leaky,  and  (b)  ferroelectric  and 
structural  properties  of  a  BiFe03-Bi(Fe0.8Sco.2)03-(Bio.8Smo.2)Fe03  pseudo¬ 
ternary  composition  spread.  Normalized  intensity  of  the  yd  Oil }  superstruc¬ 
ture  spot,  corresponding  to  the  occurrence  of  an  antiferroelectric  secondary 
phase,  is  plotted  as  a  function  of  position.  The  dotted  line  is  the  same  com¬ 
position  line  as  the  dotted  line  in  Fig.  23(a),  indicating  that  the  presence  of 
the  antiferroelectric  phase  plays  a  role  in  controlling  the  leakage  conduction 
in  the  films.  Reproduced  with  permission  from  D.  Kan  et  al.,  Integr. 
Ferroelectr.  Ill,  116  (2009).  Copyright  2009  Taylor  and  Frances,  Inc. 


square-shaped  hysteresis  loops  coincide  with  the  composition 
region  having  the  greatest  superstructure  spot  intensity.  Also, 
a  closer  look  at  the  trend  shows  that  in  the  composition  region 
where  only  Sc  is  substituted,  the  hysteresis  loops  are  charac¬ 
teristic  of  leaky  ferroelectric s,  indicating  that  Sm  is  more 
effective  than  Sc  in  reducing  the  leakage  current.  As  the  Sm 
concentration  reaches  14%,  the  ferroelectric  hysteresis  loops 
exhibit  a  transition  to  double  hysteresis  loops,  clearly  shown 
Fig.  23(a).  The  complexity  of  this  ferroelectric  system  high¬ 
lights  the  efficacy  of  the  combinatorial  approach  for  screening 
the  entire  composition  space,  and  correlating  structural  and 
ferroelectric  effects. 

Bulk  piezoelectric  materials  are  required  for  many  appli¬ 
cations,  and  bulk  composition  spread  techniques  have  been 
developed;94,95  such  techniques  can  be  used  to  investigate 
how  microstructural  properties  change  as  a  function  of  com¬ 
position  in  bulk  diffusion  couple  samples,  where  two  end 


member  compositions  are  hydrostatically  compacted  and  sin¬ 
tered  to  induce  diffusion  and  a  resultant  composition  gradi¬ 
ent.  Fig.  24  shows  the  variation  in  the  electric-field  induced 
transition  in  a  bulk  ceramic  composition  spread,  where  a 
subtle  substitution  effect  of  a  small  amount  of  La  and  Fe  into 
the  A-site  of  Pb(Zr0.6Ti0.4)O3  was  systematically  investi¬ 
gated.94  A  bulk  library  has  also  been  used  to  study  the  effect 
of  processing  temperature  on  the  structural  variation  of  pie¬ 
zoelectric  materials.96 

The  last  decade  has  seen  an  explosion  of  research  on 
multiferroic  materials,  which  are  defined  as  exhibiting  two 
or  more  ferroic  properties  (ferroelectricity,  ferromagnetism, 
and  ferroelasticity)  in  a  single  material.  Such  materials  can 
exhibit  interesting  coupling  between  the  different  ferroic 
properties,  leading  to  intriguing  novel  device  applications 
for  sensors,  actuators,  tunable  high  frequency  communica¬ 
tion  devices,  and  novel  logic  devices.97  Among  different 
types  of  multiferroic  materials,  nanocomposite  materials 
have  been  well  studied,  since  they  have  the  potential  to  real¬ 
ize  RT  device  operation  due  to  cross  coupled  ferroelectric/ 
ferromagnetic  properties.  To  study  the  formation  of  nano¬ 
composite  multiferroic  thin  films,  PLD  based  CCS  libraries 
were  designed  to  directly  probe  continuous  “mixing”  of  a 
ferroelectric  material  and  a  ferromagnetic  material  on  an 
MgO  (100)  substrate,98  Fig.  25.  One  end  of  the  spread  was 
ferroelectric  PbTi03,  and  the  other  end  was  the  ferromag¬ 
netic  insulator  CoFe204.  A  modified  version  of  the  epitaxial 
layer-by-layer  wedge  deposition  method  was  used  to  synthe¬ 
size  this  library.99  SMM  was  used  to  map  the  dielectric  con¬ 
stant,  and  RT  superconducting  quantum  interference  device 
(SQUID)  microscopy  the  magnetization,  as  a  function  of 


FIG.  24.  Bulk  composition  spread.  XRD  patterns  near  the  111  reflections  as 
a  function  of  electric  field  amplitude,  (a)  La-rich  region,  (b)  co-doped  La/Fe 
region,  and  (c)  Fe-rich  region  of  (Lai_xFex)o.o2Pt>0.98Zr0  6Ti0.403.  Vertical 
scale  is  linear  intensity.  Bottom  right  shows  a  photograph  of  the  sample  with 
the  positions  of  the  data  points  indicated.  Reprinted  with  permission  from 
J.  L.  Jones  et  al.,  Appl.  Phys.  Lett.  93,  152904  (2008).  Copyright  2008 
American  Institute  of  Physics. 
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i  CoFe204  PbTi03 

—  Average  Composition  [%] 

FIG.  25.  Multiferroic  properties  of  a  (CoFe204)i_x(PbTi03)x  CCS  library  as 
a  function  of  average  percentage  of  PbTi03.  Remnant  magnetization  was 
measured  by  scanning  SQUID  microscopy  and  high  dielectric  constant 
region  was  mapped  by  microwave  microscopy.  Region  in  the  middle  was 
found  to  display  coexisting  magnetic-ferroelectric  properties  (replotted  from 
data  in  Murakami  et  al9S).  Reprinted  with  permission  from  Murakami  et  al., 
Appl.  Phys.  Lett.  87,  112901  (2005).  Copyright  2005  American  Institute  of 
Physics. 

composition  in  the  CCS  library  (PbTi03)i_x(CoFe204)x. 
Figure  25,  from  Murukami  et  al .,98  shows  that  there  is  a 
region  in  the  middle  of  the  spread,  which  displays  coexisting 
magnetism  and  ferroelectricity  (indicated  by  enhanced 
dielectric  constant).  A  detailed  micro  structural  investigation 
has  revealed  that  different  compositions  displayed  different 
nanocomposite  nanostructures.  In  a  follow-up  study,  SMM 
was  used  to  measure  the  magnetoelectric  coupling  coeffi¬ 
cient.100  The  (PbTi03)i_x(CoFe204)x  experiments  demon¬ 
strated  that  one  can  study,  using  high  throughput  techniques, 
libraries  of  disparate  crystal  structures  (PbTi03  is  a  perov- 
skite,  while  CoFe204  is  a  spinel).  This  body  of  work  also 
illustrates  how  one  can  explore  the  synthesis  of  multifunc¬ 
tional  materials  in  general  where  the  two  end  members  may 
display  a  wide  variety  of  functionalities  beyond  the  three 
ferroic  properties.  Recently,  a  similar  technique  was  used  to 
systematically  investigate  synthesis  of  self-assembled  verti¬ 
cal  nanocomposites  in  the  BiFe03-CoFe204  system.101 
Figure  26  shows  the  scanning  electron  microscopy  (SEM) 
images  taken  at  different  average  composition  across  the 
spread.  Depending  on  the  average  composition,  the  strain  in 


the  film  and  the  magnetization  were  found  to  continuously 
vary.  This  work  showed  that  in  such  structures,  magnetoelas¬ 
ticity  dominates  the  magnetic  properties,  and  illustrates  how 
one  can  tune  and  design  the  microstructure  and  resulting 
properties  of  composites. 

IV.  MAGNETIC  MATERIALS 

Magnetic  phenomena  have  historically  played  profound 
roles  in  the  advancement  of  human  civilization;  the  techno¬ 
logical  foundation  of  modern  society  is  built  upon  develop¬ 
ments  of  magnetic  materials.102  Today,  magnetic  materials 
and  devices  are  ubiquitous  in  everyday  life.  For  example,  the 
steady  rate  at  which  the  maximum  energy  product,  (BHmax), 
of  permanent  magnets  has  increased  over  the  past  decades 
led  directly  to  continuous  improvements  in  electric  motor 
technology;  further,  magnetic  memory,  a  key  technology  for 
information  technology,  recording  and  sensing  devices  is  at 
the  heart  of  a  $100  x  109  industry. 

Demand  for  higher  performance  magnetic  materials  in  a 
wide  range  of  technologies  continues,  and  these  applications 
represent  ideal  platforms  for  implementing  combinatorial  strat¬ 
egies.  At  the  same  time,  fundamental  scientific  investigations 
into  novel  magnetic  phenomena  continue  to  uncover  exotic 
effects  such  as  multiferroicity  and  topological  dielectricity. 
Exploration  of  such  effects  and  their  underlying  physical 
mechanisms  is  directly  tied  to  mapping  of  their  phase  dia¬ 
grams.  To  this  end,  CCS  experiments  have  been  exploited  to 
discover  new  compositions  with  improved  properties,  as  well 
as  to  understand  composition- structure-property  relationships. 

A.  High  throughput  magnetic  characterization 
techniques 

A  number  of  non-destructive  scanning  probe  techniques 
have  been  used  to  map  magnetic  properties  from  combinato¬ 
rial  libraries.103-108  As  is  the  case  for  combinatorial 
measurements  in  general,  a  significant  challenges  lies  in 
extracting  accurate  quantitative  values  of  magnetic  proper¬ 
ties  from  a  particular  location  on  a  library  film,  such  that  the 
value  is  in  agreement  with  that  obtained  from  a  discrete 


FIG.  26.  SEM  images  of  the  top  view  of 
the  him  made  from  the  single  target 
and  the  composition  spread  samples  of 
BiFe03  (BFO)-CoFe204  (CFO)  (1-5). 
As  average  composition  ratio  of  BFO 
and  CFO  is  varied,  surface  morphology 
is  seen  to  change  continuously.  Reprinted 
with  permission  from  N.  M.  Aimon 
et  al.,  Appl.  Phys.  Lett.  100,  092901 
(2012).  Copyright  2012  American 
Institute  of  Physics. 
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sample  of  the  same  composition,  using  a  standard  “bulk” 
measurement  technique  such  as  a  SQUID  magnetometer. 
Significant  variations  in  some  physical  properties  between 
bulk  and  thin  film  materials  of  the  same  composition  are  of¬ 
ten  the  cause  of  the  lack  of  agreement. 

While  there  exist  many  important  magnetic  properties  to 
be  characterized,  the  net  (i.e.,  remnant)  magnetization  is  by 
far  the  most  important  property  for  screening  of  novel  mag¬ 
netic  materials.  Various  techniques  can  be  used  to  obtain 
magnetic  field  distributions  emanating  from  thin  film  combi¬ 
natorial  samples.  They  include  scanning  SQUID  microscopy 
(SSM)104,107,108  and  the  scanning  Hall  probe  technique.106 
Even  though  magnetic  held  distribution  itself  does  not  yield 
the  value  of  the  magnetization,  in  principle,  one  can  extract 
it  using  calibration  samples.  In  the  absence  of  an  external 
applied  held,  remnant  magnetization  may  be  measured. 

The  magnetic-optical  Kerr  effect  (MOKE)  tech¬ 
nique103,105,106  is  based  on  the  fact  that  the  polarization  of 
light  is  rotated  upon  reflection  from  the  surface  of  a  magnetic 
material,  and  further  that  the  rotation  angle  is  directly  pro¬ 
portional  to  the  magnetization.  MOKE  is  particularly  useful 
since  one  can  map  the  magnetic  hysteresis  loop  from  differ¬ 
ent  positions  on  combinatorial  samples  by  sweeping  the 
applied  magnetic  held,  from  which  one  can  obtain  the  satura¬ 
tion  magnetization  (through  calibration)  as  well  as  the  coer¬ 
cive  held  values.103 

The  scanning  Hall  probe  can  also  be  used  to  map  satura¬ 
tion  magnetization  values.106  In  order  to  obtain  absolute  val¬ 
ues  of  magnetization  from  the  measured  magnetic  held 
without  calibration,  an  inversion  technique  has  been  devel¬ 
oped,  which  can  be  applied  to  any  spatially  resolved  held 
distribution  data.109  This  technique  uses  a  computational 
algorithm,  which  performs  an  inverse  Fourier  transform  to 
calculate  magnetic  pole  densities,  which  in  turn  is  integrated 
to  obtain  magnetization  values  at  any  location  on  the  library. 
This  algorithm  has  been  applied  to  analysis  of  RT  SSM  data 
from  a  variety  of  magnetic  material  libraries.107,108,110  In 
particular,  it  has  been  applied  to  thin  him  composition 
spreads  of  ternary  metallic  alloy  systems  containing  ferro¬ 
magnetic  shape  memory  alloy  (SMA)  phases.  The  remnant 
magnetization  values  determined  for  the  CCS  library  were 
demonstrated  to  be  in  good  agreement  with  values  measured 
by  SQUID  and  vibrating  sample  magnetometers  for  single 
composition,  large  scale  thin  films.108 

B.  Magnetic  oxides 

Since  the  advent  of  high  temperature  superconductivity 
in  layered  perovskite  cuprates,  there  has  been  an  explosion  of 
activity  in  all  TM  oxides.  The  rich  variety  of  properties 
observed  in  metal  oxides,  such  as  colossal  magnetoresistance, 
spin  polarization,  and  long  spin  relaxation  times,  provides  a 
versatile  platform  for  exploring  novel  device  applications. 

Various  sets  of  shadow  mask  designs  have  been 
designed  for  synthesis  of  pseudobinary  and  pseudotemary 
epitaxial  CCS  libraries,111,112  Fig.  27.  This  method  was 
used104  to  fabricate  epitaxial  CCS  libraries  for  metal  oxide 
thin  hlms  on  lattice  matched  substrates  and  enabled  interrog¬ 
ation  of  the  entire  La!_xSrxMn03  composition  space  on  a 


FIG.  27.  Schematic  for  fabrication  of  ternary  composition  spreads  using  a 
set  of  shadow  masks.  The  mechanism  consists  of  three  changeable  targets,  a 
fixed  substrate,  and  a  rotating  mask  system  placed  between  them.  The  mask 
has  three  patterns,  A,  B,  and  C.  The  patterns  A  and  B  are  used  for  ternary 
composition  spread  deposition,  and  the  pattern  C  is  for  binary  composition 
spread  deposition.  The  y  axis  is  along  the  radial  direction  of  the  mask  rota¬ 
tion.  Reproduced  with  permission  from  Y.  Yamamoto  et  al.,  Appl.  Surf.  Sci. 
223,  9  (2004).  Copyright  2004  Elsevier  BV. 

single  substrate,  Fig.  28.  The  Ending  of  infrared  (IR)  reflec¬ 
tivity  at  about  0. 1  eV  is  accompanied  by  a  substantial 
increase  in  magnetization  as  x  increases  from  0.3  to  a  critical 
value  of  0.5,  as  can  be  seen  in  Figs.  28(c)-28(e).  This  verifies 
the  idea  that  a  double  exchange  ferromagnetic  interaction 
develops  in  a  cooperative  manner  at  a  critical  doping  con¬ 
centration  range.  From  Figs.  28(a)  and  28(b),  the  known 
magnetic-electronic  phases  can  be  identified  in  a  straightfor¬ 
ward  manner.  The  phase  diagram  of  Fig.  28(b)  is  in  good 
agreement  with  the  reported  bulk  phase  diagram. 

In  a  similar  experiment,  a  precursor  film  synthesis  tech¬ 
nique  for  mapping  the  phase  diagram  of  the  La!_xCaxMn03 
system  on  a  single  chip  was  employed.113  The  CCS  library 
chip  was  made  by  first  depositing  gradient  thicknesses  of  the 
precursor  layers  La203,  Mn304,  and  CaMn03  at  RT  on  a 
(100)  FaA103  substrate,  followed  by  annealing  at  200  °C  for 
several  days,  400  °C  for  30  h,  and  sintering  at  1000  °C  for 
2h  in  flowing  02;  a  predominantly  epitaxial  library  film 
resulted.  To  map  the  magnetic  phases,  low  temperature  SSM 
was  carried  out  at  3K  to7  K,  and  to  track  the  qualitative 
change  in  electrical  properties  of  this  CCS  film,  SMM  was 
performed  at  RT.81  The  various  property  maps  of  the 
La!_xCaAMn03  library  are  shown  in  Fig.  29.  The  correlation 
between  microwave,  optical,  and  magnetic  properties  at  low 
temperatures  reveals  “boundaries”  between  composition 
regions  where  properties  change  abruptly.  The  phase  boun¬ 
daries  observed  here  were  found  to  coincide  with  predicted 
boundaries  that  describe  the  spin  order  in  the  system  as  a 
function  of  composition  using  a  generalized  Hubbard 
model.114  At  x  =  0.5,  there  is  an  abrupt  change  in  the  optical 
image,  where  A Bz  exhibits  a  minimum.  This  composition  is 
also  believed  to  display  half  charge  ordering,  i.e.,  the  Mn 
orbitals  are  ordered  in  the  structure  because  they  are  evenly 
split  between  +3  and  +4  valences.113  At  x  =  0.93,  where 
A Bz  reaches  another  minimum,  the  optical  image  and  the 
SMM  data  abruptly  change  in  a  concomitant  manner.  This  is 
believed  to  be  a  transition  to  an  electrically  insulating 
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FIG.  28.  Structural-magnetic-electronic  properties  of  a  Lai_xSrxMn03  CCS 
library  as  a  function  of  x.  (a)  A  concurrent  XRD  spectra  ( 9/26  scan)  at  RT. 
The  perpendicular  axis  represents  the  26  angle.  The  colored  contour  lines 
denote  logarithmic  intensity  of  the  XRD  peak.  The  scale  bar  denotes  length 
on  the  him  along  the  x  axis,  (b)  The  out  of  plane  lattice  constants  calculated 
from  (a)  (blue  curve),  (c)  The  strength  of  the  out  of  plane  magnetic  held 
|ABz|  measured  with  scanning  SQUID  microscope  at  3K.  (d)  The  contour 
map  of  infrared  reflectivity  at  293  K  and  (e)  at  5  K.  Color  bars  denote  magni¬ 
tude  of  the  reflectivity.  Reprinted  with  permission  from  T.  Fukumura  et  al., 
Appl.  Phys.  Lett.  77,  3426  (2000).  Copyright  2000  American  Institute  of 
Physics. 

antiferromagnetic  state.  Magnetic  domain  distributions 
measured  at  several  spots  on  the  library  by  low  temperature 
scanning  SQUID  are  also  in  agreement  with  the  changing 
magnetic  states  across  the  phase  diagram. 

Tracking  magnetic  transitions  as  a  function  of  composi¬ 
tion  is  an  extremely  important  aspect  of  mapping  magnetic 
phase  diagrams.  Direct  simultaneous  visualization  mapping 
of  the  magnetic  transition  across  an  entire  composition  spread 
sample  was  demonstrated  using  magneto-optic  imaging  with 
an  indicator  film.115  By  taking  an  image  of  the  library  spread 
film  with  an  indicator  film  over  it,  Fig.  30,  one  can  readily 
identify  the  ferromagnetic  region,  which  exhibits  a  dark  edge 
due  to  stray  field  (in  the  case  of  an  in-plane  aligned  film).  By 
tracking  the  changing  width  and  position  of  the  dark  edge  as 
a  function  of  temperature,  one  can  easily  “see”  the  Curie  tem¬ 
perature  of  the  entire  phase  diagram,  schematically  shown  in 
Fig.  31  for  the  case  of  a  La!_xCaxMn03  library  film.115  The 
Curie  temperature  mapping  across  this  sample  (diamonds)  is 
superimposed  over  a  previously  published  bulk  phase 
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FIG.  29.  Electronic  and  magnetic  transitions  in  the  Lai_xCaxMn03  pseudo¬ 
binary  system  studied  using  SEMP  and  scanning  SQUID  microscope:  (a)  a 
RT  charge- coupled  device  (CCD)  photograph  of  the  continuous 
Lai_xCaxMn03  thin  him  taken  under  white  light  and  a  magnetic  phase  dia¬ 
gram  of  the  system  determined  from  bulk  single  crystals  (the  various  states 
are:  paramagnetic  insulator  (PI),  ferromagnetic  insulator  (FI),  ferromagnetic 
metal  (FM),  charge-ordered  insulator  (COI),  and  antiferromagnetic  insulator 
(AFI));  (b)  SEMP  line-scan  profiles  of  microwave  loss  Al/Q  (related  to  con¬ 
ductivity)  and  relative  frequency  shift  Af/f  (measurement  performed  at  RT); 
(c)  scanning  SQUID  line  profile  obtained  at  3K  (the  oscillation  indicates 
magnetic  domains);  and  (d)  scanning  SQUID  images  taken  at  7  K  at  four  dif¬ 
ferent  composition  regions  showing  the  magnetic  domains  and  the  transition 
from  strong  to  weak  magnetization.  Reproduced  with  permission  from  Y.  K. 
Yoo  et  al.,  Phys.  Rev.  B  63(22),  224421  (2001).  Copyright  2001  American 
Physical  Society. 

diagram,  where  various  acronyms  indicate  different  magnetic 
states.  Although  shifted  compared  to  the  bulk  curve,  one  can 
clearly  discern  the  transition  curve  from  the  paramagnetic  to 
the  ferromagnetic  state,  as  well  as  the  sharp  composition 
boundary  at  x  =  0.5  between  the  ferromagnetic  and  antiferro¬ 
magnetic  states.  Because  thin  film  samples  are  known  to  ex¬ 
hibit  property  deviations  compared  to  bulk  samples  due  to 
strains  and  presence  of  defects,  etc.,  the  observed  shift  is  not 
unreasonable.  This  experiment  clearly  illustrates  that  one  can 
use  a  simple  technique  to  map  Curie  temperature  across  com¬ 
position  spread  libraries.  In  another  example,  the  above  men¬ 
tioned  magneto-optical  technique  has  been  used  to 
successfully  map  ultraviolet  (UV)Msible  light  properties  in 
the  spinel  system  (Zn,Co)Fe204.116 

One  of  the  most  powerful  examples  of  the  efficacy  of 
combinatorial  methodologies  is  the  discovery  of  magnetic 
semiconducting  oxides  in  Co-TiC^.1 17  Ferromagnetic  semi¬ 
conductors,  formed  by  doping  magnetic  impurities  into  host 
semiconductors,  are  key  materials  for  spintronics,  where  the 
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FIG.  30.  Schematic  showing  how  a  magnetic  phase  diagram  is  constructed 
using  magneto-optic  imaging  using  an  indicator  film  technique).  The  length 
and  location  of  dark  lines  at  different  temperatures  indicate  a  magnetic  phase 
region  on  a  composition-temperature  diagram.  Reproduced  with  permission 
from  M.  J.  Turchinskaya  et  al.,  J.  Mater.  Res.  19,  2546  (2004).  Copyright 
2004  Materials  Research  Society. 

correlation  between  the  spin  and  the  charge  of  electrons 
leads  to  spin  dependent  functionalities.  Previously,  most 
known  magnetic  semiconductors  had  been  non-oxides. 
Matsumoto  et  a/.118  fabricated  thin  film  libraries  of  3d  TM 
doped  Ti02  thin  films.  When  they  were  screened  with  an 
SSM,  Co  doped  Ti02  was  found  to  exhibit  ferromagnetism 
at  RT  while  maintaining  optical  transparency,  in  a  single 
phase  microstructure,  up  to  8%  Co.  This  high  throughput  dis¬ 
covery  gave  birth  to  the  field  of  magnetic  semiconducting 
oxides;  since  this  discovery,  exploration  for  novel  magnetic 
semiconducting  oxides  in  a  variety  of  materials  systems  has 
become  an  international  research  endeavor.  High-throughput 
experimentation  continues  to  play  an  important  role  in  deci¬ 
phering  the  exact  mechanism  by  which  doped  magnetic  spe¬ 
cies  brings  about  ferromagnetism.119,120 


Lai  _xCaxMnC>3 


Ca,x 

FIG.  31.  Constructed  magnetic  phase  diagram  of  the  Lai_xCaxMn03  com¬ 
position  spread  on  a  SrTi03  substrate).  The  experimental  points  (diamonds) 
are  superimposed  on  the  bulk  phase  diagram.  PM,  FM,  NM,  AF,  CAF,  CO, 
and  FI  refer  to  paramagnetic,  ferromagnetic,  non-magnetic,  antiferromag¬ 
netic,  canted  antiferromagnetic,  charge  ordered,  and  ferromagnetic  insulat¬ 
ing  states,  respectively. 


C.  Magnetic  metallic  alloys  and  magnetic  phase 
diagrams 

The  search  for  novel  magnetic  materials  lead  to  some  of 
the  original  high  throughput  experiments  and  techniques.13,121 
Magnetic  alloys  were  screened  for  several  different  functional 
properties,  targeting  specific  applications  such  as  recording 
heads  and  inductor  cores,  using  co- sputtering  techniques 
to  rapidly  perform  phase  mapping.  In  one  example,121  the 
(Ni^Fe^)]  _x(Si02)x  pseudoternary  system  was  co-sputtered 
using  the  segmented  target  approach.  The  schematic  in 
Fig.  32  shows  the  three  disk  segments  (Ni,  Fe,  and  Si02) 
placed  on  one  radio  frequency  (RF)  sputtering  target  elec¬ 
trode.  Analysis  of  the  deposition  profile  resulting  from  the 
co-sputtering  of  the  three  segments  indicated  that  a  large  frac¬ 
tion  of  the  phase  diagram  of  this  pseudoternary  system  could 
be  achieved  in  the  library  film.  Resistance,  magnetoresistance, 
permeability,  and  magnetization  were  measured  across  the 
CCS  library;  the  latter  was  measured  by  MOKE.  At  least,  one 
other  example  of  magnetic  properties  mapped  as  a  function  of 
composition  for  a  codeposited  CCS  library  exists.122  A  differ¬ 
ent  approach  for  fabricating  CCS  films  has  been  reported  by 
Priyadarshini  et  al .,123  by  using  a  combination  of  evaporation 
filaments  to  co-deposit  alloy  libraries  containing  as  many  as 
four  elements.  Because  the  filaments  do  not  take  up  much 
space,  the  library  can  be  deposited  in  an  area  less  than  0.5  in. 
in  diameter,  as  compared  to  typical  CCS  libraries  fabricated 
by  co-sputtering,  which  are  3  in.  in  diameter  or  larger. 


Y-  AXIS 


TARGET 


FIG.  32.  Schematic  illustration  of  two  methods  for  fabricating  combinatorial 
thin  films:  (a)  co-sputtering  from  a  three  component  sputtering  target  to  map 
phase  diagrams.  Reproduced  with  permission  from  J.  J.  Hanak,  J.  Mater. 
Sci.  5,  964  (1970).  Copyright  1970  Springer;  and  (b)  co-evaporation  system 
using  electron  beam  vapor  deposition  for  mapping  a  ternary  phase  diagram. 
Reprinted  with  permission  from  K.  Kennedy  et  al.,  J.  Appl.  Phys.  36,  3808 
(1965).  Copyright  1965  American  Institute  of  Physics. 
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To  date,  less  than  5%  of  all  possible  ternary  metal  phase 
diagrams  have  been  investigated;  even  fewer  systems  exist  for 
which  magnetic  properties  have  been  mapped.  Phase  diagrams 
mapped  using  thin  CCS  film  libraries  must  be  regarded  as 
“thin  film  phase  diagrams.”  Composition-stmcture-property 
relationships  deduced  from  thin  film  libraries  can  vary  from 
those  obtained  using  bulk  samples.  Importantly,  most  high 
technology  devices  comprise  highly  integrated  multilayer  thin 
film  structures,  and  thin  film  phase  diagrams  are  often  more 
relevant  than  bulk  phase  diagrams  in  understanding  and  mod¬ 
eling  these  devices.  While  mapping  of  magnetic  properties 
alone  can  provide  great  insight,  correlating  these  properties 
with  phase  structure  across  composition  is  much  more  valua¬ 
ble.  Such  correlations  have  only  been  recently  carried  out;  the 
key  facilitating  factors  for  such  an  integrated  approach  are  the 
advent  of  computer  technology  and  the  use  of  sophisticated 
measurement  instrumentation.  A  good  example  of  compre¬ 
hensive  mapping  of  the  structural  and  magnetic  properties  of 
a  phase  diagram  using  the  CCS  spread  method  is  the  experi¬ 
ment  of  Yoo  et  al .124  The  Fe-Ni-Co  system  was  deposited 
onto  a  triangular  sapphire  substrate  by  sequential  ion  beam 
deposition  of  the  elemental  metals  at  RT,  followed  by  exten¬ 
sive  annealing  in  vacuum  (600  °C  at  10-8  Torr)  to  allow 
complete  interdiffusion  of  the  layers  and  appropriate  phase 
crystallization.  XRD  was  used  to  obtain  the  structural  phase 
distribution.  A  4°  26  range  (42°-46°)  contained  all  of  the 
major  diffraction  features  in  the  phase  distribution.  For  a 
given  structural  phase,  the  X-ray  peak  intensity  allows  one  to 
track  the  compositional  range.  Regions  in  which  phases  coex¬ 
ist  could  be  identified  by  looking  for  diffraction  peaks  arising 
from  multiple  phases.  It  was  observed  that  near  the  phase 
boundaries,  the  full  width  half  maximum  of  the  peaks  tended 


to  increase.  Figure  33(a)  shows  the  distribution  of  phases 
deduced  from  the  library  sample;  fee  Ni,  bcc  Fe,  and  fee  and 
bcc  Ni-Fe  alloys  are  present.  This  result  is  consistent  with  the 
known  bulk  structural  phases  in  this  alloy  system.  The  mag¬ 
netic  properties  of  the  library  film  were  investigated  by 
MOKE,  Fig.  33(b).  A  qualitative  saturation  magnetic  moment 
of  each  point  of  the  ternary  phase  diagram  was  mapped  at 
applied  magnetic  fields  of  +/— 50  Oe.  Comparison  of  the  two 
figures  shows  a  clear  correlation  between  the  phase  distribu¬ 
tion  and  the  magnetization;  the  bcc  Fe  has  the  highest  magnet¬ 
ization.  Further,  the  observed  magnetization  variation  is 
consistent  with  the  known  distribution  of  saturation  induction 
of  the  Fe-Ni-Co  system  deduced  using  the  traditional  bulk 
method,  Fig.  33(c).125  Although  MOKE  analysis  does  not 
allow  direct  quantitative  measurement  of  magnetization,  by 
comparing  the  measured  Kerr  rotation  angle  with  values 
obtained  on  calibration  samples  using  standard  methods  such 
as  SQUID  magnetometry,  one  can  convert  the  Kerr  angle  map 
to  a  magnetization  map.  In  this  particular  ternary  experiment, 
a  mechanical  hardness  map  was  also  obtained  using  nanoin¬ 
dentation,  which  revealed  that  the  Fe-Co  region,  consisting  of 
fee  and  bcc  (two  phase)  or  pure  bcc  (single  phase)  areas,  is 
harder  than  the  higher  Ni  concentration  fee  region. 

Recently,  ternary  systems  containing  Heusler  alloys  such 
as  Ni2MnGa  and  Co2MnGe  have  attracted  much  attention 
due  to  their  unique  functional  properties,  such  as  the 
ferromagnetic  shape  memory  effect  and  high  spin  polariza¬ 
tion.108,126  In  such  alloy  systems,  it  is  of  interest  to  under¬ 
stand  how  the  compositional  variation  and  structural 
properties  affect  the  magnetic  properties.  The  ferromagnetic 
shape  memory  effect  can  occur  in  ferromagnetic  alloys  that 
also  display  a  reversible  martensitic  transition.  In  one  study, 


FIG.  33.  Mapping  of  properties  of  a  Fe-Ni-Co  combinatorial  thin  film  library:  (a)  the  Fe-Ni-Co  phase  diagram  obtained  from  distribution  of  structural  phases 
deduced  from  X-ray  microdiffraction  peaks  of  fee,  bcc,  and  “Co  cubic”  phases,  and  (b)  mapping  of  magnetic  property  using  scanning  magneto-optical  Ken- 
effect  measurement.  The  Kerr  rotation  angle  measured  at  +/—  50  Oe  applied  field  is  plotted.  The  figure  at  left  shows  a  typical  hysteresis  loop.  Reproduced 
with  permission  from  Y.  K.  Yoo  et  al.,  Intermetallics  14,  241  (2006).  Copyright  2006  Elsevier  BV.  (c)  Saturation  induction  mapped  for  the  Fe-Ni-Co  ternary 
phase  diagram  using  the  traditional  individual  alloys  for  comparison.  Reproduced  with  permission  from  R.  M.  Bozorth,  Ferromagnetism  (Wiley-IEEE  Press, 
New  Jersey,  1993).  Copyright  1993  Wiley-IEEE  Press. 
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structural  properties  of  the  Ni-Mn-Ga  system  were  measured 
in  order  to  map  the  ferromagnetic  SMA  composition  regions, 
to  understand  the  origin  of  the  coexistence  of  these  two  prop¬ 
erties,108  Fig.  34.  CCS  libraries  were  deposited  in  an  ultrahigh 
vacuum  (base  pressure  ~  10-9  Torr)  magnetron  cosputtering 
system.  Ni,  Mn,  and  Ni2Ga3  were  used  as  the  three  targets. 
XRD  of  the  fabricated  films  revealed  that  they  were  textured, 
with  the  (110)  orientation  normal  to  the  substrate.  WDS  was 
used  to  accurately  map  the  composition  variation  of  every 
wafer.  It  was  confirmed  that  by  adjusting  the  power  applied 
to  each  cathode,  and  the  gun-substrate  distance  (typically 
12  cm),  different  regions  of  the  ternary  phase  diagram  (as 
well  as  the  entire  ternary  diagram)  could  be  mapped  reprodu- 
cibly.127  A  RT  SSM  was  used  for  rapid  identification  of  ferro¬ 
magnetic  regions,  with  conversion  of  the  field  distribution 


information  into  quantitative  magnetization  values.109  The 
result  of  this  calculation  was  plotted  to  obtain  an  RT  map  of 
remnant  magnetization  across  the  phase  diagram,  Fig.  34(a). 
In  the  figure,  the  solid  triangle- shaped  curve  marks  the  com¬ 
positional  region  mapped  on  this  particular  library  film.  The 
most  strongly  magnetic  region  stretches  from  near  the  middle 
of  the  phase  diagram  toward  the  slightly  Ni-rich  region.  The 
circle  indicates  the  region  surrounding  the  Heusler  composi¬ 
tion,  Ni2MnGa,  which  had  previously  been  extensively 
studied.128-130  Thus,  it  was  shown  that  this  composition  is 
only  one  part  of  a  large  region  exhibiting  a  high  magnetiza¬ 
tion.  As  one  moves  away  from  this  general  region,  the  rem¬ 
nant  magnetization  decreases.  The  magnetization  peak  is 
observed  near  the  center  of  the  phase  diagram,  around  the 
half  Heusler  composition  NiMnGa.  The  RT  magnetization 


Mn 


(d )  Magnetization  at  room  temperature  (e.m.u.  cm-3) 


FIG.  34.  Exploration  of  magnetic  SMAs  in  the  Ni-Mn-Ga  ternary  system,  (a)  RT  mapping  of  remnant  magnetization  of  the  Ni-Mn-Ga  system.  (The  region 
inside  the  triangular  shaped  curve  is  the  compositional  region  mapped  on  a  CCS  library  film;  the  circle  marks  the  compositions  near  the  Ni2MnGa  Heusler 
composition);  (b)  superimposed  functional  phase  diagram  deduced  from  magnetic  mapping  and  SMA  mapping  (the  shaded  region  has  the  compositions  with 
average  electron/atom  ratio  7. 3-7. 8;  the  dotted  line  surrounds  the  region  of  SMAs;  in  the  ferromagnetic  region,  the  red  area  has  the  highest  magnetization);  (c) 
800  °C  isothermal  section  of  the  Ni-Mn-Ga  ternary  system;  (d)  Martensitic  transformation  temperature  plotted  against  RT  saturation  magnetization  (Each  data 
point  correspond  to  a  composition  on  the  spread  wafer.  The  line  is  a  linear  fit  to  the  data).  Reprinted  (Figs.  34(a),  34(b),  and  34(d))  with  permission  from  I. 
Takeuchi  et  al..  Nature  Mater.  2,  180  (2003).  Copyright  2003  Macmillan  Publishers  Ltd. 
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values  extracted  here,  ~ 200-3 00  emu/cm3,  are  consistent 
with  saturation  magnetizations  obtained  using  a  vibrating 
sample  magnetometer  on  single  composition  bulk  samples.  In 
the  same  experiment,  CCS  libraries  were  deposited  on  wafers 
containing  micromachined  arrays  of  cantilevers  to  map  the 
SMA  transition  temperature  regions.  Figure  34(b)  summa¬ 
rizes  the  magnetic  and  the  martensite  phase  diagrams 
deduced  from  the  combined  data.  It  is  clear  that  there  is  a 
large  region  outside  of  the  near  Heusler  composition  that 
exhibits  both  ferromagnetic  and  reversible  martensite  behav¬ 
iors.  It  has  been  shown  that  this  class  of  material  exhibits  a 
martensitic  instability  for  stoichiometries  where  the  (s  +  d 
electron)/atom  ratio  is  ~7.4.  ’  The  shaded  area  in  Fig. 

34(b)  covers  the  composition  region  where  the  (s  +  d  elec- 
tron)/atom  ratio  is  7. 3-7. 8,  and  a  large  overlap  between  the 
shaded  area  and  the  observed  reversible  martensite  region  can 
be  observed.  The  RT  magnetization  is  plotted  against  the 
martensitic  start  temperature  in  Fig.  34(d)  for  each  composi¬ 
tion  in  the  large  composition  range  studied.  A  clear  trend 
between  the  two  parameters  emerges:  the  higher  the  magnet¬ 
ization,  the  lower  the  transformation  temperature.  From  sepa¬ 
rate  measurements,  it  has  been  confirmed  that  magnetization 
at  RT  is  roughly  proportional  to  the  Curie  temperature.  The 
plot  indicates  that  there  is  a  strong  thermodynamic  magneto- 
structural  coupling  in  this  system.  Such  a  coupling  had  previ¬ 
ously  been  observed  in  only  a  limited  compositional  range 
near  Ni2MnGa  (inside  the  circle  in  Fig.  34(a).130  This  experi¬ 
ment  revealed  that  the  relationship  holds  in  general  for  mag¬ 
netic  alloys  with  a  martensitic  instability.  A  comparison  of 
Fig.  34(b)  with  the  800  °C  isotherm  of  the  Ni-Mn-Ga  system 
Fig.  34(c)  (Ref.  131))  shows  that  it  is  the  /?  phase,  which 
largely  transforms  martensitically  at  low  temperatures.  Micro 
XRD  of  the  spread  indicated  that  much  of  the  studied  regions 
have  the  L2X  Heusler  structure,  or  a  distorted  version  of  it. 

Another  Heusler  alloy,  Co2MnGe,  was  predicted  to  be 
fully  spin  polarized  at  the  Fermi  level.  Although  to  date  the 
spin  polarized  behavior  has  not  been  confirmed,  the  possibil¬ 
ity  that  deviations  from  stoichiometry  could  cause  a  reduc¬ 
tion  in  spin  polarization  has  motivated  ternary  mapping  of 
the  Co-Mn-Ge  system.105,132,133  The  MOKE  imaging  tech¬ 
nique  was  used  to  rapidly  characterize  the  magnetic  proper¬ 
ties  of  the  entire  ternary  system.  Figure  35  shows  the 
transverse  differential  MOKE  intensity  image  of  a  25  nm 
thick  continuous  ternary  phase  diagram  chip  measured  using 
a  charge  couple  display  (CCD)  camera  while  applying  a 
magnetic  field  of  +/— 5  KOe.  In  this  figure,  areas  of  high  in¬ 
tensity  correspond  to  magnetic  regions  at  RT  with  strong 
magneto-optic  coupling.  Triangularly  cut  Ge  substrates  of 
different  orientations,  Fig.  36,  were  used  to  map  the  ternary 
phase  diagram  using  molecular  beam  epitaxy  to  deposit 
library  films.  Sub-monolayers  of  Co,  Mn,  and  Ge  were  de¬ 
posited  sequentially  using  a  combination  of  automated  mov¬ 
ing  shadow  masks,  sample  rotation,  and  source  shutters.134 
Three  regions  in  the  phase  diagram  in  Fig.  35  were  deter¬ 
mined  to  be  strongly  ferromagnetic:  (1)  near  pure  Co,  (2) 
near  the  Co-rich  region  centered  around  Co0.7Mn0.iGe0.2, 
and  (3)  in  the  middle  of  the  phase  diagram,  approximately  at 
(Co0.5Mn0.5)i_xGex,  for  x  between  0.2  and  0.6.  The  area 
showing  strong  signals  along  the  MnxGei_x  binary  edge,  for 


FIG.  35.  RT  magnetic  phase  diagram  of  a  25  nm  thick  Co-Mn-Ge  phase  dia¬ 
gram  sample  grown  by  molecular  beam  epitaxy  (MBE).  Image  of  differential 
MOKE  intensity  measured  at  +/—  5  Oe,  which  corresponds  to  the  saturated 
states  of  the  system.  Reproduced  with  permission  from  F.  Tsui  and  P.  A. 
Ryan,  Appl.  Surf.  Sci.  189,  333  (2002).  Copyright  2002  Elsevier  BY. 


x  between  0.6  and  0.9,  corresponds  to  a  known  antiferromag¬ 
netic  region.  Results  are  in  good  agreement  with  earlier 
reports  on  limited  composition  ranges  inside  the  ternary  sys¬ 
tem.  In  order  to  zoom  in  to  study  the  details  of  magnetism  in 
the  (Co!_xMnx)0.8Ge0.2  region,  a  separate  CCS  library  with  a 
smaller  composition  range  was  synthesized.  Reflection  high 
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FIG.  36.  Synthesis  schemes  of  CCS  libraries.  Cross-sectional  schematic  views 
of  (a)  full  binary  sample  (CoxMn!_x)aSii_a  with  0<x<  1  and  (b)  partial 
binary  sample  with  0.3  <  x  <  1 .  (c)  Partial  ternary  sample  and  (d)  composition 
range  of  the  partial  ternary  sample  (shaded  region)  within  the  full  ternary 
phase  diagram.  The  trilayers  of  Si-Co-Mn  submonolayer  wedges  and  “partial” 
wedges  are  typically  repeated  200  times.  Reprinted  with  permission  from  L. 
He  et  al.,  J.  Vac.  Sci.  Technol.  B  29,  03C124  (2011).  Copyright  2011 
American  Institute  of  Physics. 
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energy  electron  diffraction  (RHEED)  and  MOKE  measure¬ 
ments  were  taken  across  the  library,  and  it  was  found  that  the 
robust  ferromagnetism  observed  with  high  saturation  mag¬ 
netization  (and  high  Curie  temperature)  at  Co0.4Mn0.4Ge0.2 
has  strong  correlation  with  the  onset  of  a  structural  order- 
disorder  transition,  revealed  in  an  abrupt  change  in  the 
RHEED  intensity.132  Thus,  as  has  been  demonstrated  in  this 
paper  through  several  examples,  high  throughput  mapping  of 
compositional  phase  diagrams  using  various  characterization 
techniques  can  go  beyond  simply  mapping  phase  boundaries 
and  allows  one  to  develop  important  insights  concerning  the 
physics  of  the  underlying  magnetic  phenomena. 

Bulk  diffusion  couples  (“multiples”)135  can  also  be  used 
to  study  magnetic  properties  as  a  function  of  composition. 
Curie  temperature — composition  relationships  can  be 
obtained  from  bulk  diffusion  multiples  using  magnetic  force 
microscopy  (MFM),  as  shown  in  Fig.  37  for  the  case  of  the 
Co-Cr-Mo  ternary  system.136  Fig.  37(b)  shows  a  SEM  back- 
scattered  electron  image  of  the  Co-Cr-Mo  tri-junction  area. 


FIG.  37.  Application  of  MFM  to  map  magnetic  domains  to  obtain  Curie 
temperature — composition  relationship  for  the  Co-Cr-Mo  ternary  system: 
(a)  schematic  of  a  Co-Cr-Mo-Nb-Ni  diffusion  multiple  with  the  Co-Cr-Mo 
tri-junction  region  highlighted;  (b)  MFM  images  superimposed  on  a  back- 
scattered  electron  image  of  the  Co-Cr-Mo  tri-junction;  and  (c)  the  24  °C 
iso-Curie  temperature  compositions  of  the  fee  phase  in  the  Co-Cr-Mo 
ternary  system  obtained  by  tracing  the  boundary  (dotted  line  in  (b))  between 
the  domain/non-domain  region  with  EPMA.  Reproduced  with  permission 
from  J.  C.  Zhao,  Prog.  Mater.  Sci.  51,  557  (2006).  Copyright  2006 
Pergamon. 


Superimposed  on  the  image  are  MFM  maps  (taken  at 
^24  °C)  displaying  magnetic  domains  in  Co  rich  areas  of  the 
fee  phase;  topographic  features  such  as  polishing  scratches 
and  voids  allowed  the  exact  positions  of  the  MFM  maps  to 
be  superimposed  on  the  SEM  image.  The  ferromagnetic 
region  of  the  fee  phase  could  be  separated  from  the  paramag¬ 
netic  region,  as  shown  by  the  dotted  line  in  Fig.  37(b). 
EPMA  along  this  line  revealed  the  compositions  at  which 
the  magnetic  transition  takes  place  at  24  °C,  i.e.,  an  iso-Curie 
temperature  line  at  24  °C,  Fig.  37(c).  By  performing  MFM 
mapping  at  different  temperatures,  it  would  have  been  possi¬ 
ble  to  obtain  the  constant  Curie  temperature  surface  for  the 
fee  phase. 

While  this  review  article  is  mostly  focused  on  thin  film 
based  combinatorial  techniques,  it  is  important  to  point  out 
the  significance  of  the  diffusion  couple  (or  “multiple”)  sam¬ 
ple  techniques.  In  Sec.  Ill,  the  utility  of  a  bulk  based 
approach  to  mapping  the  properties  of  piezoelectric  ceramics 
was  discussed.  The  necessity  to  map  the  properties  of  materi¬ 
als  in  bulk  form  is  particularly  acute  in  a  variety  of  metallic 
alloy  systems,  where  subtle  change  in  composition  at  differ¬ 
ent  length  scales  can  dramatically  influence  the  mechanical 
properties.  In  addition  to  magnetic  materials,  the  bulk  map¬ 
ping  technique  is  effective  in  investigation  of  structural 
materials  such  as  superalloys.  Diffusion  multiple  samples 
can  be  used  to  determine  phase  diagrams  and  design  func¬ 
tional  alloys.137-141  However,  one  disadvantage  of  the  bulk 
diffusion  multiple  techniques  is  that  property  measurements 
must  be  carried  out  with  very  high  spatial  resolution  of  the 
order  of  1  to  5  micrometers,  due  to  the  limited  scale  of  diffu¬ 
sion  in  bulk  samples. 

The  magnetic  properties  of  thin  films  and  multilayers 
can  change  dramatically  as  a  function  of  thickness,  because 
the  length  scale  that  determines  the  exchange  interaction  of 
magnetic  layers  of  the  order  of  several  nanometers  is  small. 
In  fact,  in  the  magnetic  recording  head  industry,  subnanome¬ 
ter  dimensional  control  of  the  multilayer  thicknesses  is  nec¬ 
essary  for  state-of-the-art  technology.  Because  it  is  relatively 
straightforward  to  prepare  a  thickness  gradient  film,  the  so- 
called  “wedge”  structure  has  long  been  used  to  investigate 
the  effect  of  thickness  (typically  using  a  single  component 
film)  on  magnetic  properties.  This  approach  was  successfully 
implemented  in  a  Fe  (lnm)/Cr  (0-15nm)/Fe  sandwich/ 
wedge  structure  to  systematically  study  the  interlayer 
thickness  dependence  on  exchange  coupling  between  thin 
films.142,143  Figure  38  is  a  scanning  electron  microscopy 
polarization  analysis  (SEMPA)  image  of  the  Fe-Cr-Fe 
trilayer  with  varying  Cr  thickness,  showing  the  periodic 
reversal  of  the  sign  of  the  exchange  coupling  due  to  the 
Ruderman-Kittel-Kasuya-Yosida  (RKKY)  interaction.144,145 

Given  the  urgent  need  for  new  permanent  magnets  with 
higher  energy  products  for  next  generation  motor  vehicle 
and  wind  power  generator  applications,  high  throughput 
approaches  are  expected  to  play  an  important  role  in  the  dis¬ 
covery  of  new  magnetic  materials;  the  need  is  particularly 
acute  for  RE-free  permanent  magnets  with  high  coercive 
fields.146,147  One  proposed  approach  to  synthesizing  high 
energy  product  permanent  magnets  is  to  combine  a  high  co¬ 
ercive  field  hard  magnet  with  a  high  magnetization  soft 
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FIG.  38.  SEMPA  (scanning  electron  microscopy  with  polarization  analysis) 
image  of  a  thickness  wedge  sample,  showing  the  quasiperiodic  oscillations 
of  the  exchange  coupling  between  two  Fe  films  separated  by  a  Cr  film  of 
0-15  nm.  The  light  and  dark  contrast  represents  regions  in  which  the  mag¬ 
netization  of  the  top  Fe  layer  is  aligned  or  antialigned  with  that  of  the  bottom 
layer.  Reproduced  with  permission  from  J.  Unguris  et  al.,  Phys.  Rev.  Lett. 
67,  140  (1991).  Copyright  1991  American  Physical  Society. 

magnet  in  a  tightly  coupled  nanocomposite  configuration,  so 
that  the  overall  material  displays  the  combination  of  high  co¬ 
ercive  field  and  high  magnetization.148  Since  exchange  pin¬ 
ning  of  the  soft  layer  by  the  hard  magnet  is  the  key  to 
achieving  this  effect,  it  is  critical  to  understand  the  mecha¬ 
nisms  and  parameters  that  enable  optimum  coupling.  To  this 
end,  thin  film  multilayers  have  been  used  as  a  model  system 
to  understand  the  micromagnetic  parameters  that  determine 
the  coupling  behavior.149,150  The  gradient  thickness 
approach  has  also  been  used  to  investigate  the  thickness  of  a 
high  magnetization  Fe  layer  coupled  to  a  high  coercivity 


SmCo  layer  using  X-ray  magnetic  circular  dichroism 
(XMCD)  to  determine  the  critical  soft  layer  thickness  pinned 
by  SmCo.103  Such  research  has  further  been  extended  to 
investigate  the  relationship  between  soft  phase  micromag¬ 
netic  constants  and  thickness  in  an  Fe-Co/CoPt  bilayer 
square  library  deposited  by  e-beam  evaporation.  Using  auto¬ 
mated  shadow  masks,  the  thickness  of  the  soft  phase  was 
varied  from  1  to  lOnm  in  one  direction,  and  the  composition 
of  the  soft  phase,  FexCo!_x,  in  the  other,151  Fig.  39(a). 
Figure  39(b)  shows  the  MOKE  hysteresis  loops  taken  at  each 
position  on  the  library.  This  plot  can  be  used  to  visually 
determine  the  effect  of  the  soft  phase  composition  on  the 
exchange  coupling  length,  Fig.  39(c),  i.e.,  the  length  beyond 
which  a  region  of  the  soft  phase  is  decoupled  from  the  hard 
layer.  This  study  revealed  that  the  most  important  parameters 
determining  the  exchange  coupling  length  are  the  hard  phase 
coercive  field  and  the  soft  phase  magnetization. 

Both  ferromagnetic  SMA  and  magneto strictive  materials 
possess  magnetoelastic  properties,  i.e.,  in  response  to  an 
applied  magnetic  field,  a  strain  is  induced  in  the  material, 
and  vice  versa.  Magnetoelastic  materials  are  a  technologi¬ 
cally  important  class  of  materials  with  applications  in  sen¬ 
sors  and  actuators.  Among  magnetoelastic  materials  which 
display  Joule  magnetostriction,  Tb-Dy-Fe  (terfenol)  exhibits 
the  largest  field  induced  strain,  with  saturation  magnetostric¬ 
tion  of  more  than  1000 ppm  (0.1%  strain).  With  the  advent 
of  novel  remote  sensing  technologies  and  the  ubiquity  of 
smart  microactuators  in  microsystems,  there  is  a  need  to  find 
inexpensive,  non-brittle,  RE-free  substitutes  for  terfenol. 
Recently,  Fe!_xGax  was  found  to  display  a  significant  mag¬ 
netostriction  of  3/2(A10o)  >  400  ppm.  Subsequently,  it  was 
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FIG.  39.  (a)  Schematic  of  the  library  for 
investigating  the  coupling  behavior  of 
hard/soft  bilayers  (left),  (b)  magnetic 
hysteresis  loops  (a.u.)  for  different  soft 
layer  thicknesses  (ts)  and  compositions 
on  a  CoPt  (Hc  =  0.64  T)  hard  layer  meas¬ 
ured  on  a  single  library,  and  (c)  the 
coupling  length  zx,  experimentally 
determined  from  the  soft  layer  thickness 
at  which  the  transition  between  one 
phase-like  behavior  and  two-phase-like 
behavior  takes  place.  Reproduced  with 
permission  from  A.  J.  Zambano  et  al., 
Phys.  Rev.  B  75,  144429  (2007). 
Copyright  2007  American  Physical 
Society. 
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found  that  structural  heterogeneity  at  a  compositional  phase 
boundary  can  lead  to  formation  of  a  phase,  which  exhibits 
substantial  magnetostriction.152  A  thin  film  cantilever 
method  to  map  the  magnetos  trie  tive  properties  of  co¬ 
sputtered  Fe-Ga  and  Fe-Ga-Al  CCS  libraries  was  devised153 
in  which  a  laser  beam  reflected  from  the  tip  of  the  cantilever 
was  used  to  determine  its  deflection  as  a  function  of  applied 
magnetic  field.  The  measured  deflection  was  then  used  to 
determine  the  magnetostrictive  coefficient.  This  work 
showed  that  the  observed  bulk  magnetostriction  trend  with 
composition  could  be  reproduced  in  thin  films,  Fig.  40.  In 
another  Fe-Ga-Al  CCS  study,117  Mossbauer  spectroscopy 
was  used  to  understand  the  effects  of  Ga  and  A1  additions  to 
the  magnetic  environment  of  bcc  Fe;  Ga  atom  pairing  is 
believed  to  be  the  reason  for  magnetostriction  enhancement, 
and  this  is  directly  affected  by  A1  substitution.  In  a  related 
magnetostriction  experiment,154  it  was  found  that  upon 
annealing  and  quenching  Co-Fe  films,  a  large  magnetostric¬ 
tion  was  observed  at  compositions  at  the  bcc/(fcc  +  bcc) 
phase  boundary.  This  study  underscores  the  role  structural 
heterogeneity  plays  in  magnetostrictive  properties  and  the 
importance  of  properly  controlling  the  thermal  treatment 
process  of  thin  film  libraries  to  obtain  accurate  property 
trends  as  a  function  of  composition.155 

D.  Superconductors 

Superconductivity  remains  one  of  the  most  intensively 
studied  fields,  and  superconducting  materials  are  a  constant 
target  for  materials  exploration.  Superconducting  materials 
are  used  in  a  wide  variety  of  advanced  technology  applica¬ 
tions  ranging  from  ultrasensitive  magnetometers  and  ultra¬ 
fast  logic  circuits  to  superconducting  power  cables  and 
magnetic  levitation.  In  the  mid-1980s,  the  discovery156  of 
high  temperature  superconductivity  in  cuprates  led  to  a  revo¬ 
lution  in  the  field. 

A  chemical  precursor  based  thin  film  synthesis  method 
to  fabricate  libraries  of  cuprate  superconductors  has  been 
developed.8  This  work  was  the  first  demonstration  of  thin 


FIG.  40.  Magnetostrictive  coefficient  measured  from  Fe-Ga  thin-film  com¬ 
position  spread  plotted  against  values  from  bulk  studies  as  a  function  of  Ga 
content.  The  compositional  trend  agrees  well  with  bulk  values.  Reprinted 
with  permission  from  J.  R.  Hattrick-Simpers  et  al.,  Appl.  Phys.  Lett.  93(10), 
102507  (2008).  Copyright  2008  American  Institute  of  Physics. 


FIG.  41.  A  128-member  matrix  array  cuprate  superconductor  library,  prior 
to  sintering.  Each  site  is  1  mm  by  2  mm;  the  color  of  each  is  the  natural  color 
of  reflected  light  from  a  white  light  source.  The  amorphous  precursor  layers 
are  deposited  using  a  series  of  shadow  masks.  Reproduced  with  permission 
from  X.  D.  Xiang  et  al.,  Science  268,  1738  (1995).  Copyright  1995 
American  Association  for  the  Advancement  of  Science. 


film  libraries  of  inorganic  functional  materials  in  a  discrete 
composition  matrix  array,  Fig.  41.  Different  combinations  of 
amorphous  precursor  layers  were  deposited  using  a  series  of 
shadow  masks.  Screening  of  this  array  was  performed  using 
a  multichannel  four  terminal  resistance  measurement  system, 
and  many  of  the  matrix  elements  (consisting  of  up  to  7  ele¬ 
ments)  were  correctly  identified  as  exhibiting  superconduct¬ 
ing  transitions  from  resistance-temperature  measurements. 
This  work  not  only  demonstrated  the  utility  of  the  combina¬ 
torial  synthesis  of  libraries  of  functional  inorganic  materials 
but  also  showed  that  the  thin  film  amorphous  precursor  tech¬ 
nique  could  be  used  to  synthesize  multi-element  composi¬ 
tions  with  complex  crystal  structures. 

Several  groups  have  adopted  the  parallel  four-terminal 
resistance-temperature  measurement  method  for  rapid 
screening  of  new  superconductors.157,158  Figure  42  shows  a 
196  pin  probe  array,  which  allows  simultaneous  resistance 


FIG.  42.  The  196  pin  layout  that  makes  electrical  contact  to  a  49  sample 
thin  film  library  for  simultaneous  four  terminal  resistivity  measurements. 
The  shortest  distance  between  pins  is  4.64  mm.  Reprinted  with  permission 
from  K.  C.  Hewitt  et  al.,  Rev.  Sci.  Instrum.  76,  093906  (2005).  Copyright 
2005  American  Institute  of  Physics. 
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measurements  of  up  to  49  samples  in  a  liquid  He  cryostat. 
Another  detection  technique  that  has  been  applied  for  map¬ 
ping  superconducting  materials  is  magnetic  field  modulated 
microwave  absorption,159,160  which  takes  advantage  of  the 
fact  that  in  the  presence  of  a  modulating  magnetic  field, 
superconducting  regions  of  a  material  absorb  microwave 
power  due  to  flux  penetration.  A  common  screening 
approach  is  to  apply  this  probe  to  a  thin  film  library.159-163 

Despite  decades  of  international  research,  the  mecha¬ 
nisms  of  superconductivity  in  many  compounds  are  not  fully 
understood.  The  recent  discovery  of  high  temperature  super¬ 
conductivity  in  Fe  pnictides164  has  reignited  the  search  for 
new  superconductors,  and  high  throughput  methodologies 
will  continue  to  play  an  important  role.  Mapping  of  phase 
diagrams  can  be  an  important  part  of  deciphering  the  physics 
of  superconducting  mechanisms.  As  an  example,  the  use  of 
combinatorial  techniques  to  map  phase  space  in  the  Ba-Y- 
Cu-0  system  can  be  found  in  the  work  of  Wong-Ng  et  al.165 

V.  OPTICAL  MATERIALS 
A.  Thermochromics  and  electrochromics 

It  has  been  estimated  that  20%-30%  of  the  energy  used 
to  heat  and  cool  buildings  in  the  US  is  lost  through  windows, 
due  to  undesired  solar  heating  and  thermal  conduction.166 
“Smart”  windows  that  have  either  thermochromic  coatings 
(i.e.,  films  whose  solar  acceptance  varies  as  a  function  of 
temperature),  or  electrochromic  coatings  (i.e.,  films  whose 
solar  acceptance  varies  as  a  function  of  applied  bias)  have 
great  potential  to  minimize  building  energy  losses.  These 
coatings  are  competing  technologies  for  the  “smart”  compo¬ 
nent  of  windows;  electrochromic  windows  require  a  power 
supply,  but  are  easier  to  adjust  than  thermochromic  win¬ 
dows,  which  are  simply  sensitive  to  temperature. 

A  novel  screening  methodology  for  high  throughput 
studies  of  electro-optically  active  materials  has  been 
designed  using  an  optical  screen  containing  an  electrochro¬ 
mic  film  of  W03.167  However,  electrochromic  materials  for 
window  applications  have  not  been  studied  combinatorially. 
Thermochromic  materials  that  undergo  a  reversible  phase 
transformation  close  to  RT  are  attractive  options  for  smart 
window  coatings;  the  low  temperature  phase  of  lightly  doped 
V02  is  transparent  to  IR  radiation,  thus  allowing  natural 
heating,  while  the  higher  temperature  phase  reflects  the  IR, 
thus  limiting  excess  heating  of  the  building  when  the  exter¬ 
nal  temperature  is  high.  Both  phases  are  transparent  to  the 
visible  spectrum;  thus,  the  window’s  clarity  and  esthetic 
appeal  are  maintained.  Ideally,  phase  transition  temperatures 
(Tc)  for  thermochromic  window  coatings  would  be  in  the 
ambient  range,  i.e.,  —  20  °C  to  35  °C.  At  68  °C,  V02  has  the 
closest  Tc  of  the  unary  metal  oxides.  W,  Mo,  Ti,  Nb,  and  Fe 
dopants  can  depress  Tc  at  rates  of  10°C-25  °C/at.  %?168-170 
and  should  be  studied  combinatorially. 

Oxygen  partial  pressure  has  been  varied  combinatorially 
to  study  its  effect  on  VOx  stoichiometry  in  PLD  films.171 
The  processing  window  for  preparing  a  V02  film  from  either 
a  V203  or  V205  target  is  shown  in  Fig.  43.  Thin  V02  ther¬ 
mochromic  films  have  also  been  deposited  by  an  electrostatic 
atomization  technique.172  Although  the  paper  reports  the 
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FIG.  43.  (a)  Experimental  and  simulated  RBS  spectra  for  a  V20x  film  (from 
V203  target  in  vacuum):  x  =  2.99,  thickness  =  83  nm;  (inset)  RBS  sche¬ 
matic,  (b)  photograph  of  V20x  library  with  indicated  values  of  oxygen  par¬ 
tial  pressure  during  deposition,  and  (c)  oxygen  stoichiometry  in  the  library 
as  a  function  of  background  pressure  of  oxygen.  Reproduced  with  permis¬ 
sion  from  N.  D.  Bassim  et  al.,  Appl.  Surf.  Sci.  254,  785  (2007).  Copyright 
2007  Elsevier  B  V. 

deposition  of  film  “libraries,”  it  is  not  clear  what  parameter 
was  varied  in  the  preparation  of  these  libraries. 

B.  Luminescents  and  phosphorescents 

Luminescent  materials  emit  visible  radiation  by  a  vari¬ 
ety  of  mechanisms,  such  as  electrical  stimulation,  absorption 
of  photons,  or  a  chemical  reaction,  but  not  by  heat  (emission 
due  to  heat  is  called  incandescence).  A  subset  of  luminescent 
materials,  phosphorescents,  slowly  (scale  of  minutes  to 
hours)  emits  visible  radiation  after  absorption  of  photons. 
Many  early  combinatorial  thin  film  experiments  were 

designed  to  explore  novel  luminescent  materials,  because 
such  libraries  served  as  excellent  demonstrations  of  the  effi¬ 
cacy  of  the  combinatorial  approach;  initial  high-throughput 
screening  is  very  simple,  as  one  only  had  to  view  optically 
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FIG.  44.  Photoluminescence  images  of  a 
series  of  phosphor  libraries  processed 
under  different  conditions  with  nominal 
compositions  indicated,  (a)  La(or  GdF3)m 
(Sr)nA10x:Euy2+/3+  where  0.375  <  m  <  1, 
0.25  <  n  <  0.4,  1.88  <  y  <  12%  in  atomic 
ratio,  annealed  at  1150°C;  (b)  same  as  (a) 
but  annealed  at  1400  °C;  (c)  La(or 
GdF3)mA10x:Tby3+(Cez3+);  Euh2+/3+ 
where  0.32<m<l,  1.29  <y<  6, 

0.65  <  z  <  4%,  1 .29  <  h  <  8%  annealed  at 
1150°C;  (d)  same  as  (c),  but  annealed  at 
1400  °C;  (e)  La(or  GdF3)m(Sr)n 

A10x:Euy2+/3+  where  0.178  <  m  <  0.714, 
0.17  <  n  <  0.4,  0.75  <  y  <  16.7%  annealed 
at  1150°C.  Each  substrate  is  lin.  square. 
The  images  are  taken  with  UV  lamp  exci¬ 
tation.  Reprinted  with  permission  from  X. 
D.  Sun  et  al.,  Appl.  Phys.  Lett.  70,  3353 
(1997).  Copyright  1997  American  Institute 
of  Physics. 


excited  libraries  and  visually  detect,  which  composition  sites 
were  emitting  light. 

Phosphorescent  materials  are  essential  for  lighting, 
imaging,  and  display  applications  including  plasma,  field 
emission,  electroluminescent,  and  cathode  ray  tube  (CRT) 
projection  displays.  Although  the  photophysical  processes 
leading  to  luminescence  are  relatively  well  understood,  their 
specific  spectral  properties,  luminescence  efficiencies,  and 
operational  lifetimes  depend  on  a  complex  interplay  between 
the  excitation  source,  host  crystal  lattice,  sensitizer,  and 
luminescent  centers.  Identification  of  phosphors  that  satisfy 
the  set  of  requirements  of  a  given  display  application  can  be 
highly  empirical,  and  therefore  combinatorial  methodologies 
can  play  a  significant  role  in  optimizing  their  compositions. 
Beyond  initial  screening,  it  is  important  to  carry  out  quantita¬ 
tive  characterization  of  luminescent  materials  libraries  using 
specific  metrics,  such  as  emission  spectrum,  chromaticity  (a 
measure  of  the  quality  of  a  color  based  on  hue  and  dominant 
wavelength),  and  quantum  efficiency. 

Sun  et  al.119  were  the  first  to  show  the  utility  of  the  com¬ 
binatorial  approach  for  luminescent  materials.  Using  precur¬ 
sor  layers  deposited  by  RF  sputtering  through  a  set  of 
physical  masks,  they  fabricated  thin  film  libraries  to  investi¬ 
gate  RE  activated  refractory  metal  oxides,  Gd(La,Sr)A10x. 
Their  findings  are  summarized  in  Fig.  44;  from  this  study,  a 
green  emitting  composition,  Gdo.46Sro.31Al! 23OxF! 38:Eu0.o6 
was  identified  that  had  better  chromaticity  (as  well  as  a  high 
quantum  efficiency,  >94%),  than  that  of  the  standard 
Y3(Al,Ga)5012:Tb.  The  same  library  has  been  characterized 
by  synchrotron  micro-x-ray  beam  experiments,  where  struc¬ 
tural  properties,  composition  distribution,  and  dopant  ioniza¬ 
tion  states  were  mapped  in  a  single  experiment.180  In  a 
related  study,  a  red  phosphor  with  a  quantum  efficiency  of 
86%  was  found  in  the  (Gd2_xZnx)03_<5:Eu3+  SyStem.173  This 
phosphor  has  a  better  chromaticity  than  the  state  of  the  art 


red  phosphor,  Y203:Eu3+,  and  it  exhibits  enhanced  color 
saturation. 

A  segmented  masking  scheme,  combined  with  the  use 
of  precursor  layers,  has  proven  to  be  particularly  effective 
for  designing  phosphor  composition  libraries.  In  one  study,  a 
quaternary  scheme  was  implemented  using  photolithography 
for  investigation  of  a  large  number  of  different  host  materi¬ 
als,  combined  with  dopant  activation  elements.174  Figure  45 
shows  the  set  of  quaternary  masks  used  to  fabricate  a  1024 
element  quaternary  library.  In  this  scheme,  each  mask  is 
used  four  times,  each  time  rotated  by  90°,  to  segment  partic¬ 
ular  compositional  variations  into  four  distinct  groups.  Thus, 
using  the  five  masks,  compositions  of  the  type  ABCDE  can 
be  mapped,  where  each  letter  A,  B,  C,  D,  or  E  represents  a 
four-fold  doping  element  variation,  leading  to  1024  distinct 
combinations.174  In  this  doping  experiment,  the  mask  pattern 
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FIG.  45.  Masks  for  generating  a  quaternary  library.  Reproduced  with  per¬ 
mission  from  J.  S.  Wang  et  al.,  Science  279,  1712  (1998).  Copyright  1998 
American  Association  for  the  Advancement  of  Science. 
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was  created  using  a  photolithographic  lift-off  technique,  and 
deposition  of  the  precursor  materials  was  by  PLD,  followed 
by  a  series  of  thermal  treatments  to  interdiffuse  the  precursor 
layers  and  crystallize  the  compounds.  A  scanning  spectrome¬ 
ter  was  used  for  quantitative  screening  of  the  library;  a 
well-characterized  phosphor,  Zn2Si04:Mn,  was  used  as  a  ref¬ 
erence  standard.  In  this  study,  a  blue  luminescent  material 
always  sought  by  the  display  manufacturing  community, 
Gd3Ga5012/Si02,  was  discovered  combinatorially. 

A  variety  of  other  techniques  have  been  implemented  for 
synthesizing  libraries  of  luminescent  materials,  including  co¬ 
sputtered  CCS181,182  and  a  solution-based  technique  using  a 
multihead  scanning  ink-jet  delivery  system,175  Fig.  46.  Such 
inkjet  delivery  systems  are  broadly  applicable  for 
combinatorial  synthesis  of  other  materials  systems.  A  key 
advantage  of  the  inkjet  delivery  system  is  that  precursors  can 
be  mixed  at  the  molecular  level  in  the  nozzle  prior  to  con¬ 
trolled  delivery  at  each  site,  thus  avoiding  the  need  for 
thermal  interdiffusion  treatment  of  the  library.183  In  another 
study,176  electron  beam  deposition  in  a  moving  shutter  sys¬ 
tem  was  used  to  create  a  large  library  containing  composition 
gradients  of  multiple  phosphor  hosts  as  well  as  different  RE 
dopants,  Fig.  47.  This  work  led  to  the  discovery  of  a  new 
RE  phosphor,  Sr2Ce04,  containing  one-dimensional  chains  of 
Ce06  octahedra  in  their  crystal  structure.  A  more  recent  com¬ 
binatorial  experiment  on  ZnO-based  phosphors184  involved 
lattice  engineering  via  precise,  unit  cell  deposition  enabled 
by  automated  mask  motion.  The  unit  cell  engineered  deposi¬ 
tion  method  results  in  an  extremely  well  controlled  library, 
and  as  a  result,  luminescent  behavior  can  be  attributed  to 
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FIG.  47.  Image  of  UV  excited  (254 nm)  photoluminescence  from  a  library. 
The  blue-white  emission  in  the  upper  right  corner  is  from  Sr2Ce04. 
Reproduced  with  permission  from  E.  Danielson  et  al.,  Science  279,  837 
(1998).  Copyright  1998  American  Association  for  the  Advancement  of 
Science. 


precise  positions  on  the  library.  A  series  of  new  phosphors 
including  blue  (ZnO:W,  Mg,  and  ZnO:W),  red  (ZnO:Y,  Eu), 
and  yellow  (ZnO:V)  have  been  identified  from  this  study. 
Ink-jet  printing  has  also  been  successfully  used  to  screen  and 
optimize  polymeric  materials  for  organic  light-emitting 
diodes  and  other  optoelectronic  device  applications.185-187 

ZnO  is  a  versatile  wide  bandgap  semiconductor  material; 
to  date,  a  variety  of  devices  including  light  emitting  pn- 
junctions,  lasers,  gas  sensors,  and  photodetectors  have  been 
demonstrated  using  ZnO  and  related  compounds.188  The 
combinatorial  approach  has  proven  to  be  an  effective  way  to 
help  expand  the  utility  of  this  material.  By  controlling  the 
number  of  pulses  used  to  ablate  ZnO  and  Zno.9Mgo.1O 
targets,  bandgap  engineering  has  been  demonstrated  in  an 
epitaxial  library  of  Zn!_xMgxO  thin  films.189  Fig.  48  shows 
the  deposition  cycle  scheme,  which  was  used  to  create  a  dis¬ 
crete  library  consisting  of  nine  different  members  with  differ¬ 
ent  Mg  substitution  concentrations.  In  another  study,  a  CCS 
library  spanning  the  entire  Zn!_xMgxO  binary  system  was 
fabricated  using  a  layer-by-layer  epitaxial  deposition 
approach.190,191  Figure  49  shows  the  continuous  9-29  XRD 
plot  taken  across  the  entire  composition  spread,  showing  the 
evolution  of  structural  phases.  From  this  study,  complete 
bandgap  mapping  of  the  Zn!_xMgxO  was  carried  out,  and  the 
presence  of  a  mixed  phase  region  in  the  middle  of  the  compo¬ 
sition  region  was  revealed.  The  same  CCS  library  was  used 
to  demonstrate  ultraviolet  photodetector  arrays,  where  each 
device  in  the  array  has  a  slightly  different  peak  response 
wavelength  due  to  the  continuously  changing  bandgap. 


FIG.  46.  Schematic  of  a  scanning  multihead  inkjet  delivery  system  for  solu¬ 
tion  deposition  of  optical  material  libraries.  The  system  is  integrated  from 
drop-on-demand  single  nozzle  piezoelectric  inkjets  (see  Ref.  435). 
Reproduced  with  permission  from  T.  X.  Sun,  in  Combinatorial  Materials 
Synthesis ,  edited  by  X.  D.  Xiang  and  I.  Takeuchi  (Marcel  Dekker,  New  York, 
2003),  pp.  141-176. 


C.  Transparent  conducting  oxides  (for  photovoltaics) 

Transparent  conducting  oxides  (TCOs)  enable  a  variety 
of  important  applications  including  photovoltaic  (PV)  cells, 
touch  screens,  flat  panel  displays,  and  light  emitting  diodes. 
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Sample  No. 

FIG.  48.  Schematic  diagram  of  a  combinatorial  laser  molecular  beam  epi¬ 
taxy  deposition  system.  Reproduced  with  permission  from  Y.  Matsumoto 
et  al.,  Jpn.  J.  Appl.  Phys.,  Part  2  38,  L603  (1999).  Copyright  1999  The  Japan 
Society  for  Applied  Physics. 

Their  primary  function  in  photovoltaics  is  as  a  transparent 
electrode,  either  at  the  front  or  back  of  the  cell,  for  collection 
of  photogenerated  charge  carriers.  Therefore,  the  most  im¬ 
portant  properties  of  a  TCO  are  low  resistivity  (<5  x  10-4 
Q-cm,  or  high  conductivity,  2000  (Q-cm)_1)  and  high  trans¬ 
parency  in  the  solar  spectrum  range  (>80%  from  300  nm  to 
1200nm  wavelengths).  Secondary  desired  properties  include 
thermal,  chemical  and  interfacial  stability,  high  carrier 


FIG.  49.  X-ray  diffraction  of  a  MgxZni_xO  composition  spread  in  the  26 
range  of  30°  to  50°.  The  vertical  scale  is  arbitrary,  but  shows  the  relative  in¬ 
tensity  of  the  peak  at  each  26  value  at  a  fixed  composition.  The  upper  left 
shows  the  schematic  of  the  spread  with  the  lattice  constants  of  the  end  com¬ 
positions.  Reprinted  with  permission  from  I.  Takeuchi  et  al.,  J.  Appl.  Phys. 
95,  3840  (2004).  Copyright  2004  American  Institute  of  Physics. 


mobility,  work  function  matching  with  desired  photovoltaic 
materials,  and  the  ability  to  grow  textured  films;192  further 
TCO  selection  criteria  are  based  on  price,  availability,  and 
material  toxicity.  Given  the  large  number  of  elements  that 
may  form  transparent  oxides,  and  the  sensitivity  of  the  final 
material  properties  to  the  presence  of  dopants  as  well  as 
processing,  the  search  for  optimal  TCOs  encompasses  a  large 
multidimensional  parameter  space,  ideal  for  high  throughput 
experimentation.  Since  the  first  combinatorial  experiments 
on  TCOs,193  much  subsequent  research  has  been  devoted  to 
the  search  for  new  or  improved  materials.  The  majority  of 
the  work  to  date  has  been  focused  on  the  ln203  (TIO),  ZnO, 
and  SnO  systems.  Although  this  review  focuses  primarily  on 
TCOs  for  photovoltaic  applications,  there  has  been  recent 
interest  in  searching  for  new  absorber  materials  using  high- 
throughput  density  functional  theory  (DFT)  methods194  and 
combined  high-throughput  DFT  and  experimental  determi¬ 
nations  of  novel  hole  transport  layers.195,196 

1 .  ln203-based  TCOs 

ln203-based  TCOs  are  used  widely  in  display  applica¬ 
tions  and  as  heterojunctions  with  intrinsic  layer  photovoltaic 
cells.  Their  widespread  usage  is  a  result  of  their  high  trans¬ 
parency  and  high  conductivity  (3-10  x  103  (Q-cm)_1)  when 
in  the  crystalline  state.  However,  indium  (In)  is  a  costly  raw 
material,  and  is  difficult  to  dope;  a  systematic  method  for 
doping  ln203  that  preserves  or  enhances  its  properties  has 
been  pursued  for  more  than  10  yr.  Tin  (Sn)  doped  TIO  has 
been  explored  extensively  by  traditional  materials  science 
approaches.  Recently,  using  libraries  prepared  by  liquid  titra¬ 
tion  followed  by  high  temperature  annealing,  Sn  doping  of 
ln203  to  form  highly  crystalline  films  was  reported.197  The 
highest  transmission  reported  in  this  study  was  56.4%  at 
550  nm,  substantially  lower  than  values  reported  elsewhere. 

Extensive  research  has  been  performed  in  alloying  ln203 
with  zinc  (Zn)  and  Sn  to  form  zinc  indium  oxide  (ZIO)  and 
TIO,  which  are  amorphous  TCO  materials.  Although  they 
typically  exhibit  lower  conductivities  than  crystalline  materi¬ 
als,  amorphous  films  are  of  interest  because  they  do  not 
require  annealing  at  high  temperature.  In  a  comprehensive 
body  of  work,193,198-200  CCS  libraries  were  deposited  in  a 
specially  designed  reactive  sputtering  chamber,  using  an  off 
axis  sputtering  configuration.  Multiple  depositions  were  nec¬ 
essary  to  fully  cover  the  compositional  phase  diagram. 
Library  characterization  was  accomplished  with  a  suite  of 
high  throughput  tools  including  EMPA,  UV/visible  (VIS)/ 
near  infrared  (NIR)  and  Fourier  transform  infrared  spectros¬ 
copy  (FTIR)  measurements,  Hall  probe,  and  micro-XRD.  In 
the  ZIO  system,  a  broad  maximum  in  conductivity  was  found 
to  exist  at  Zn  concentrations  between  45  and  80  at.  %.201  The 
conductivity  maximum  was  found  by  Hall  measurements  to 
be  ^3000  (Q-cm)_1  for  In  concentrations  of  ^67%,  Fig.  50. 
The  region  of  enhanced  film  conductivity  roughly  correlates 
with  the  amorphous  region.  All  films  exhibited  optical  trans¬ 
missivities  of  >80%  in  the  visible  region.  It  was  also  found 
that  a  simple  Drude  approximation  could  be  used  to  correlate 
measurements  of  the  plasma  wavelength  to  the  electron  car¬ 
rier  concentration  variations. 
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FIG.  50.  Comparison  of  conductivity  (a),  carrier  concentration  (b),  and  mo¬ 
bility  as  measured  by  combinatorial  tools  (closed  symbols)  and  discrete  Hall 
measurements  (open  symbols).  A  broad  maximum  in  conductivity  and  mo¬ 
bility  is  observed  in  the  region  determined  to  be  amorphous  by  XRD. 
Reproduced  with  permission  from  M.  P.  Taylor  et  al.,  Meas.  Sci.  Technol. 
16,  90  (2005).  Copyright  2005  Institute  of  Physics. 


There  has  also  been  interest  in  doping  ln203  with  TM 
such  as  Ti  and  Mo  to  further  increase  conductivity.  In 
In-Mo-0  libraries,  it  was  found  that  for  films  deposited  at 
350  °C,  the  maximum  conductivity  was  only  1000  (Q-cm)-1 
for  low  concentrations  (~6  at%)  of  Mo;198  this  maximum  is 
attributable  to  the  opposing  compositional  trends  of  carrier 
concentration  and  electron  mobility.  Samples  annealed  at 
550  °C  showed  a  six-fold  increase  in  conductivity  to  6000  (Q 
cm)-1,  although  such  thermal  treatment  makes  the  material 
incompatible  with  organic  photovoltaics  and  raises  some 
questions  about  stability  when  deposited  on  Si-  or  CdTe- 
based  photovoltaics.  In-Ti-0  films  were  also  reported  in  the 
same  study;  in  this  case,  samples  annealed  at  500  °C  exhib¬ 
ited  a  maximum  conductivity  of  6000  (Q  cm)-1  for  Ti  con¬ 
centrations  between  3  and  4  at.  %.200  Between  1.5  and  3  at. 
%  Ti  doping,  a  linear  relationship  between  Ti  concentration 
and  carrier  density  was  observed,  expressed  as  1.06  elec- 
trons/Ti  atom,  Fig.  51.  This  is  about  a  factor  of  five  larger 
than  the  trend  seen  for  Mo,  and  suggests  that  Ti  is  a  promis¬ 
ing  dopant  for  ln203. 

Cosputtering  through  a  mask  was  used  to  create  a  library 
of  discrete  zinc  indium  tin  oxide  (ZITO)  CCS  films  on  poly¬ 
meric  substrates.202  The  samples  were  deposited  from  a  pure 
ZnO  target  and  an  ln203  target  doped  with  10  wt.  %  Sn.  The 
Van  der  Pauw  method  was  used  to  measure  the  resistivity, 
carrier  concentration,  and  mobility  of  each  discrete  film. 
XRD  revealed  samples  that  were  amorphous  across  the 
entire  composition  region  investigated.  A  maximum  conduc¬ 
tivity  of  2.4  x  103  (Q-cm)-1  was  observed  at  a  Zn  content  of 
8.9  at.  %;  this  value  is  comparable  to  values  for  crystalline 
ln203.  A  negligible  compositional  dependence  on  carrier 
concentration  was  observed,  decreasing  only  slightly  for 
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FIG.  51.  Conductivity,  carrier  concentration,  and  mobility  measured  as  a 
function  of  Ti  content  in  ln203  films.  The  carrier  concentration  varies  line¬ 
arly  over  a  range  of  2  to  4  at.  %  Ti,  with  the  conductivity  peaking  in  this 
region.  Reprinted  with  permission  from  M.  van  Hest  et  al.,  Appl.  Phys.  Lett. 
87,  032111  (2005).  Copyright  2005  American  Institute  of  Physics. 

increased  Zn  content.  In  contrast,  the  carrier  mobility  was 
observed  to  show  a  strong  dependence  on  the  overall  Zn  con¬ 
tent.  For  all  compositions  investigated,  the  optical  transmis¬ 
sion  of  ZIO  in  the  wavelength  range  300  nm  to  2400  nm  was 
observed  to  be  larger  than  the  for  conventionally  grown  ITO 
films.  Finally,  PLD  was  used  to  synthesize  libraries  in  the 
In2-2XMe2x03  system,  where  Me  =  Ti,  Zr,  and  Sn.203  It  was 
found  that  Ti  and  Zr  additions  to  ln203  yielded  superior  opti¬ 
cal  properties  compared  to  those  incorporating  Sn,  without 
compromising  resistivity. 

2.  Sn02-based  TCOs 

Sn02  is  the  most  commonly  utilized  TCO  film  material, 
due  to  its  large  band  gap  (4.1  eV)  and  the  high  natural  abun¬ 
dance  of  Sn.  Sn02  is  generally  doped  with  either  fluorine  (F) 
or  antimony  (Sb)  to  increase  its  conductivity  to  about 
1.7  x  103  (Q-cm),-1  although  this  is  still  more  than  a  factor 
of  two  smaller  than  for  ITO.  Despite  its  greater  resistance, 
the  low  cost  of  Sn  has  led  to  a  large  effort  in  identifying 
optimal  doping  and  processing  conditions  for  Sn02. 
Optimization  of  F-doped  Sn02  has  been  reported  using 
atmospheric  pressure  chemical  vapor  deposition  (APCVD), 
with  SnCl4  and  aqueous  hydrogen  fluoride  as  precursors.204 
The  goal  of  this  study  was  to  optimize  the  growth  of  the 
Sn02:F  using  high  throughput  techniques,  followed  by 
implementation  of  the  optimized  films  in  photovoltaic  cells. 
The  combinatorial  APCVD  experiment  was  realized  through 
the  use  of  a  movable  susceptor,  which  heated  local  areas  of  a 
large  wafer;  this  enabled  serial  growth  of  films  at  different 
temperatures  and  precursor  mixtures.  Two  different  types  of 
samples  were  grown  in  this  study.  In  the  first  set,  the  TCOs 
were  grown  on  glass  substrates  for  measurement  of  resistiv¬ 
ity,  transparency,  surface  roughness,  and  optical  absorption 
of  the  Sn02  processed  under  different  conditions.  It  was 
found  that  under  appropriate  conditions,  the  sheet  resistance 
could  be  lowered  to  acceptable  levels,  with  optical 
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absorption  losses  of  1%  at  500  nm  wavelength.205 
Additionally,  it  was  found  that  optimal  film  morphologies 
could  be  obtained  through  the  controlled  addition  of  metha¬ 
nol  and  t-butanol  206  In  the  second  set  of  samples,  the  deposi¬ 
tion  of  the  TCO  layer  was  followed  by  the  growth  of  a 
standard  a-Si:H  single  junction  p-i-n  solar  cell,207  Fig.  52. 
Then,  the  full  quantum  efficiency  of  a  library  of  solar  cells 
incorporating  Sn02  grown  at  different  conditions  was  meas¬ 
ured  using  a  calibrated  light  source  (i.e.,  the  light  source  tests 
the  solar  cell  response  vs.  the  solar  spectrum  when  the  direct 
optical  path  length  through  the  Earth’s  atmosphere  is  1.5 
times  as  long  as  the  path  length  at  the  zenith).  Energy  con¬ 
version  efficiencies  of  up  to  9.3%  with  high  short-circuit  cur¬ 
rents  (14.9mA/cm2)  and  external  quantum  efficiency  (i.e., 
the  proportion  of  photons  that  produce  collected  charge  car¬ 
riers)  values  of  86%  at  540  nm  wavelength  were  achieved. 

3.  ZnO-based  TCOs 

Other  commonly  used  TCOs  are  based  on  the  ZnO  sys¬ 
tem.  With  appropriate  doping  ZnO  has  been  reported  to  have 
conductivities  as  high  as  5  x  103  (Q-cm)_1,  comparable  to 
ITO,  and  Zn,  like  Sn,  is  much  less  expensive  than  In.  An  addi¬ 
tional  benefit  is  that  ZnO  exhibits  high  thermal  and  chemical 
stability.  To  date,  there  have  been  combinatorial  studies  for 
various  dopants  (gallium  (Ga),  aluminum  (Al),  nitrogen  (N), 
and  cobalt  (Co))  in  ZnO,  and  the  alloying  of  ZnO  with  SnO. 
Experiments  were  carried  out  on  ZnO  doped  with  N  and  Al  to 
demonstrate  in  situ  materials  processing  monitoring  of  TCOs 
during  annealing.193  It  was  found  that  high  spatial  resolution 
(10  fim  spot  size)  fiber  optic  reflectance  measurements  of  indi¬ 
vidual  spots  on  a  library  during  annealing  could  be  obtained 
in  less  than  15  ms,  Fig.  53.  This  measurement  system  was 
then  applied208  to  investigate  the  properties  of  ZnO  alloyed 
with  Sn02,  in  compositions  along  the  ZnO-Sn02  tie  line.  Two 
local  maxima  in  conductivity  were  reported,  for  Zn/Sn  ratios 
of  2:1  and  1:1,  corresponding  to  the  formation  of  the  amor¬ 
phous  compositions  Zn2Sn04  and  ZnSn03;  both  phases 
showed  good  transmission  (^80%).  Other  researchers197  used 
sol-gel  methods  to  dope  ZnO  with  Al.  After  annealing  at 
1073  K,  XRD  peaks  attributed  to  ZnAl204  were  identified  for 
Al  contents  as  low  as  2  at.  %,  and  UV/VIS  measurements 
demonstrated  transmissions  as  high  as  80%  at  550  nm.  Using 


(a)  (b) 


FIG.  52.  Images  of  F  doped  Sn02  samples  deposited  by  APCVD  on  (a)  a 
glass  substrate  and  (b)  as  a  series  of  solar  cells.  Reproduced  with  permission 
from  U.  Dagkaldiran  et  al..  Mater.  Sci.  Eng.,  B.  159-160,  6  (2009). 
Copyright  2009  Elsevier  S.A. 
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FIG.  53.  Schematic  of  high  spatial  resolution  in  situ  reflectance  measure¬ 
ment  system.  Reproduced  with  permission  from  Q.  Wang  et  al.,  Appl.  Surf. 
Sci.  189,  271  (2002).  Copyright  2002  Elsevier  BY. 


robot  assisted  sol-gel  methods,  a  series  of  mixtures  of  Zn-Sn- 
O  along  the  same  tie  line  were  prepared  in  further  work.209 
Their  results  showed  that  in  contrast  to  PVD  and  CVD  meth¬ 
ods,  sol-gel  samples  showed  unacceptably  low  conductivities. 

Substitutional  doping  of  ZnO  with  Ga  via  combinatorial 
sputtering  was  also  studied,210  with  maximum  Ga  contents 
up  to  8  at.  %.  The  position  of  the  conductivity  peak  was 
found  to  shift  to  lower  Ga  concentrations  with  increasing 
annealing  temperature,  with  a  maximum  value  of  1.4  x  10-3 
(Q-cm)_1  for  samples  containing  7.5  at.  %  Ga  annealed  at 
250  °C,  Fig.  54.  Further,  transmissivities  were  greater  than 
80%  for  the  entire  library.  Jung  et  al211  reported  doping  of 
highly  oriented  (002)  ZnO  films  doped  with  Ga203.  The 
libraries  were  grown  on  glass  substrates  at  RT,  and  focused 
on  compositions  with  less  than  13.6  wt.  %  Ga203.  They 


FIG.  54.  Compositional  variation  of  conductivity  in  ZnO  samples  doped 
with  Ga  at  different  deposition  temperatures.  Reproduced  with  permission 
from  C.  W.  Gorrie  et  al.,  33rd  IEEE  Photovoltaic  Specialists  Conference 
(2008),  Vol.  1-4,  p.  635.  Copyright  2008  IEEE. 
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report  an  order  of  magnitude  decrease  in  the  conductivity, 
from  ^103  (Q  cm)-1  for  the  lowest  Ga203  concentrations,  to 
^102  (Q  cm)-1  in  samples  with  7.0  wt.  %,  Fig.  55.  The  best 
sample  contained  0.8  wt.  %  Ga203  and  had  a  conductivity 
6.8  x  102  (Q  cm)-1  with  an  optical  transmittance  at  550  nm 
of  90%.  A  direct  comparison  of  the  maximum  conductiv¬ 
ities210,211  of  materials  with  similar  compositions  shows  a 
factor  of  two  differences  in  the  values  reported,  while  the 
values  and  trends  for  optical  transmission  compare  well  at 
RT. 

In  an  effort  to  explore  p-type  doping  of  ZnO,  combina¬ 
torial  studies  of  Co  doping  of  ZnO  have  been  conducted. 
Co  has  a  similar  ionic  radius  to  Zn  and  thus  has  a  large  sol¬ 
ubility.  Further,  there  is  the  potential  that  the  d6  electron 
configuration  of  Co3+  in  an  octahedral  crystal  environment 
could  aid  in  providing  a  conduction  path  for  hole  carriers 
at  the  top  of  the  valence  band.  Suh  et  al.  produced  a  series 
of  combinatorial  libraries  at  different  oxygen  partial  pres¬ 
sures  and  correlated  structural  information  gathered  by 
XRD  with  composition  through  the  use  of  principle  com¬ 
ponents  analysis.212  From  this,  it  was  determined  that  the 
presence  of  oxygen  vacancies  in  the  samples  was  related 
not  only  to  the  oxygen  incorporation  rate  during  deposition 
but  also  to  the  overall  Co  content  in  the  film.  The  highest 
conductivity  samples  were  those  of  the  structure  type 
ZnCo204. 

4.  Other  TCO  systems 

Combinatorial  libraries  of  the  IR  transparent,  TCO 
NiO-Mn203-Co304  ternary  system  were  screened  for  their 
conducting  properties  as  a  function  of  temperature.213 
Conductivity  measurements  were  performed  with  80  /im 
resolution  by  using  a  carefully  selected  sharp  graphite  rod, 
which  was  kept  in  constant  contact  with  the  film  during  scan¬ 
ning.  By  controlling  the  hardness  of  the  graphite,  electrical 
contact  could  be  maintained  without  scratching  the  film 
surface  or  damaging  the  probe  tip.  They  found  that  the  conduc¬ 
tivities  of  a  series  of  ternary  alloys  along  the  line  Mn0.45Ni0.63 
Co1>9204-Mn0.6oNio.72Co1-6804-Mn0.69Nio.8iCo1>504,  while 
comparable  to  NiCo204,  exhibited  much  lower  temperature 
sensitivity  (about  1.5%/K). 


FIG.  55.  Resistivity  and  film  thickness  as  a  function  of  position  on  a 
Zn0-Ga203  CCS  library.  Reproduced  with  permission  from  K.  Jung  et  al., 
Appl.  Surf.  Sci.  256,  6219  (2010).  Copyright  2010  Elsevier  BY. 


Further  research  used  PLD  to  deposit  CCS  libraries  of 
Nb-Sn-Ti-O.214  Their  results  showed  that  films  with  Sn  lev¬ 
els  less  than  20  at%  maintained  low  resistivity  and  high 
transparency,  while  reducing  the  refractivity  of  the  coating 
(a  reduction  in  the  refractivity  reduces  the  effective  path 
length  of  light  through  a  material,  thus  reducing  the  probabil¬ 
ity  of  absorption).  Using  metal-organic  chemical  vapor  depo¬ 
sition  (MOCVD),  linearly  varying  CCS  libraries  of  Co-Sn-0 
were  deposited.215  In  this  work,  transmission  greater  than 
80%  from  500  nm  to  lOOOnm  wavelengths  across  the  entire 
library  was  reported.  Using  standard  Hall  measurement  tech¬ 
niques,  high  carrier  concentrations  (1020cm-3)  and  mobili¬ 
ties  (60  cm2  V-1  s_1)  were  found  for  all  compositions. 


D.  Phase  change  materials 

Phase-change  materials  (PCMs),  most  often  semicon¬ 
ducting  or  semi-metallic  alloys  (also  known  as  chalcogenide 
compounds)  containing  S,  Se,  or  Te,  have  both  optical  and 
electrical  storage  applications.216  Upon  heating,  PCMs 
undergo  a  reversible  phase  transformation  from  a  super¬ 
cooled,  metastable  amorphous  phase  to  a  crystalline  phase. 
The  amorphous  phase  has  significantly  lower  reflectivity 
than  the  crystalline  phase,  making  them  the  material  of 
choice  for  rewriteable  optical  storage  media  (e.g.,  compact 
discs).217,218  PCMs  also  have  a  low  electrical  resistance  in 
the  crystalline  phase,  which  may  increase  by  many  orders  of 
magnitude  in  the  amorphous  phase,  making  them  very  use¬ 
ful  for  storage  applications  as  non-volatile  flash  memories 
in  integrated  circuits.219  In  fact,  the  original  intended 
application  for  PCMs  was  for  integrated  circuit  memory 
storage.220,221 

To  achieve  high  write  and  erase  rates  for  memory  stor¬ 
age,  the  kinetics  of  the  phase  transformation  between  the 
amorphous  and  crystalline  phases  (in  both  directions)  must 
be  very  fast.  Non-combinatorial  crystallization  kinetics  stud¬ 
ies  led  to  the  selection  of  Ge-Sb-Te  ternary  alloys  for 
optical  storage  applications.  ’  Alloys  for  nonvolatile 
memory  applications  also  require  fast  kinetics,  but  other  fac¬ 
tors  such  as  low  power  consumption  and  cycling  endurance 
are  also  important.  Ge2Sb2Te5  is  the  alloy  of  choice  for  this 
application  at  this  time.  However,  since  a  large  range  of 
composition  space  is  being  explored  for  this  application,  it  is 
interesting  to  note  that  there  is  a  dearth  of  combinatorial 
experimentation  in  this  area.  It  is  likely  that  the  data  required 
for  this  application  are  not  readily  obtained  by  combinatorial 
methods,  or  that  much  of  this  work  is  proprietary  and  has  not 
been  published.  One  group  performed  combinatorial  investi¬ 
gations224  on  sputter  deposited  Ge-Sb-Te  films.  The  CCS 
film  library  had  concentration  gradients  centered  about  the 
Ge2Sb2Te5  phase.  Crystallization  times  that  varied  between 
220  and  500  ns,  and  write  times  that  varied  between  20  and 
120  ns,  both  as  a  function  of  composition,  were  measured. 
Figure  56  illustrates  a  map  of  amorphization  (writing)  time 
as  a  function  of  composition.  Further  work  by  this  group225 
expanded  the  materials  optimization  to  CCS  libraries  based 
around  three  other  compounds:  GeSb4Te7,  GeSb2Te4,  and 
Ge4SbTe5.  It  was  determined  that  ternary  compositions  at  or 
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FIG.  56.  Map  of  amorphization  (writing) 
time  on  a  Ge-Sb-Te  CCS  film  library. 
Reproduced  with  permission  from  S. 
Kyrsta  et  al.,  Thin  Solid  Films  398,  379 
(2001).  Copyright  2001  Elsevier  S.A. 


near  the  GeTe-Sb2Te3  pseudobinary  tie  line  exhibit  fast 
(70-150  ns)  crystallization  behavior. 

VI.  ENERGY-RELATED  MATERIALS 
A.  Thermoelectrics 

Thermoelectric  materials,  which  enable  the  direct  solid 
state  conversion  of  thermal  to  electrical  energy,  have  impor¬ 
tant  applications  for  utilization  of  waste  heat  energy  (particu¬ 
larly  vehicular  engine  heat)  via  the  Seebeck  effect,  as  well  as 
solid  state  refrigeration  via  the  Peltier  effect.226,227  However, 
at  present,  the  best  thermoelectric  conversion  materials  and 
devices  are  not  efficient  enough  for  widespread  commercial 
use.  The  availability  of  higher  conversion  efficiency  thermo¬ 
electric  materials  will  play  a  significant  role  in  their 
implementation.  High  conversion  efficiency  thermoelectric 
devices  require  materials  that  possess  large  figures  of 
merit,  ZT,  which  is  equal  to  S2ctT/k ,  where  S  the  Seebeck 
coefficient,  a  the  electrical  conductivity,  k  the  thermal  con¬ 
ductivity,  and  T  the  absolute  temperature.  Combinatorial 
approaches  have  been  used  to  optimize  and  discover  thermo¬ 
electric  materials  through  measurement  of  S,  a,  and  therefore 
S 2o  (another  thermoelectric  figure  of  merit,  known  as  the 
power  factor).  To  obtain  ZT,  the  thermal  conductivity  of  the 
film  must  also  be  measured.  This  is  a  notoriously  difficult 
measurement  to  perform,  especially  on  a  combinatorial  thin 
film  library,  because  the  mass  of  the  substrate  is  far  greater 
than  that  of  the  film;  however,  high  throughput  measure¬ 
ments  of  thermal  conductivity  have  also  been  performed,  as 
will  be  discussed  below. 


Using  PLD,  a  CCS  library  film  of  (Ca1_xBax)3Co409 
was  prepared,  where  x  varied  between  0  and  0.2. 228  The 
Seebeck  coefficient  and  resistivity  of  the  striped  and  pat¬ 
terned  films  in  the  library  were  then  measured  using  a  multi¬ 
channel  four  point  probe  array,  Fig.  57,  where  one  side  of 
the  film  stripe  was  heated  to  create  the  required  AT  for 
measurement,  and  the  other  side  was  maintained  at  RT. 
Barium  substitution  was  found  to  significantly  decrease 


Thermocouple 


FIG.  57.  Schematic  diagram  of  a  multichannel  power  factor  measurement  sys¬ 
tem  for  thermoelectric  materials.  Reproduced  with  permission  from  H.  Minami 
et  al.,  Appl.  Surf.  Sci.  197,  442  (2002).  Copyright  2002  Elsevier  BY. 
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FIG.  58.  Schematic  diagram  of  measurement  probe  to  measure  electrical 
conductivity  and  Seebeck  coefficient.  Reprinted  with  permission  from  M. 
Otani  et  al.,  Appl.  Phys.  Lett.  91,  132102  (2007).  Copyright  2007  American 
Institute  of  Physics. 

the  conductivity,  but  had  no  effect  on  the  Seebeck  coeffi¬ 
cient.  The  Sr  and  La  substituted  Ca3Co409  system, 
(Ca1_x_ySrxLay)3Co409,  where  x  varied  continuously 
between  0  and  approximately  0.35,  has  also  been  studied 
combinatorially. 229,230  PLD  was  used  to  deposit  the  ternary 
film  library,  and  a  four  point  screening  probe,  with  one 
heated  leg,  Fig.  58,  was  used  to  measure  both  the  Seebeck 
coefficient  and  the  electrical  resistivity.  The  probe  spatial  re¬ 
solution  was  about  1  mm;  over  one  thousand  points  in  the 
library  were  measured  in  about  6  h.  Seebeck  coefficient  cali¬ 
bration  films  measured  with  this  probe  were  found  to  be  in 
good  agreement  with  those  measured  using  a  commercial 
thermoelectric  measurement  apparatus.  The  rapid  construc¬ 
tion  of  a  thermoelectric  property  diagram  is  illustrated  in 
Fig.  59  for  the  Ca3_x_ySrxLayCo409  system.  One  can  see 
that  the  power  factor  “sweet  spot”  is  easily  identified  at 
about  x  =  0.25  and  y  =  0.10.  Such  high  throughput  thermo¬ 
electric  materials  analysis  will  likely  facilitate  commerciali¬ 
zation  of  this  class  of  materials. 

Another  high  throughput  technique  for  thermoelectric 
oxides  has  been  reported231,232  in  which  thick  film  discrete 
libraries  were  synthesized  by  depositing  aqueous  solutions  of 
M!-M2-Co-0  (M  =  metal)  on  an  alumina  plate,  followed  by 


FIG.  60.  Graph  of  Seebeck  coefficients  from  the  La-Ni-O  library  heated  at 
1073  K  for  10  h  in  air.  Reproduced  with  permission  from  R.  Funahashi  et  al., 
Appl.  Surf.  Sci.  223,  44  (2004).  Copyright  2004  Elsevier  BV. 

annealing.  The  macroscopic  films,  approximately  one  centi¬ 
meter  in  length,  were  then  probed  with  a  handheld  two  termi¬ 
nal  Seebeck  probe.  This  technique  enabled  the  screening  of 
up  to  about  1000  samples  per  day.  An  effective  example  of 
its  use  for  the  La-Ni-0  system  is  shown  in  Fig.  60. 

Another  class  of  thermoelectric  materials  that  was  stud¬ 
ied  is  the  MgxSiyGei_y  ternary  system,  where  x  was  varied 
between  2.3  and  4.5,  and  y  between  0.45  and  1,  on  a  discrete 
PLD  film  library  sample.233  The  library  was  deposited  onto  a 
prefabricated  alumina  substrate  into  which  the  local  heater 
and  contact  leads  had  been  integrated.  Temperature  gradients 
across  the  sample,  as  well  as  the  base  temperature  of  the 
samples  (which  varied  between  RT  and  673  K),  were  meas¬ 
ured  with  an  IR  camera.  The  measurements  were  first  cali¬ 
brated  using  a  constantan  film,  whose  Seebeck  coefficient 
and  resistivity  were  found  to  be  very  close  to  the  bulk  values. 
Seebeck  coefficient  and  resistivity  results  for  this  ternary 
library  film  are  depicted  in  Fig.  61. 

More  traditional  thermoelectric  metal  alloys  have  also 
been  studied  using  high  throughput  methods.234  Ni!_xCux 
library  samples,  where  x  spanned  the  entire  range  from  0  to 
1,  were  prepared  in  two  different  ways.  The  first  set  was  pre¬ 
pared  by  powder  metallurgy  techniques  involving  ball  mill¬ 
ing  of  eleven  different  compositions,  followed  by  spark 


Ca3Co40g  (Ca2La)Co40g  CagC^Og  (Ca2La)Co40g  Ca3Co40g  (Ca2La)Co40g 


FIG.  59.  (a)  Electrical  conductivity,  (b)  Seebeck  coefficient,  and  (c)  power  factor  of  the  library  film  (Cai-^Sr^La^Co^c),  where  x  and  y  vary  between  0 
and  approximately  0.35.  Reprinted  with  permission  from  M.  Otani  et  al.,  Appl.  Phys.  Lett.  91,  132102  (2007).  Copyright  2007  American  Institute  of  Physics. 
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FIG.  61.  Compositional  dependence  of 
the  Seebeck  coefficient  and  electrical  re¬ 
sistivity  at  different  temperatures  for  the 
Mg-Si-Ge  library,  (a-d)  and  (1-4)  refer 
to  various  compositional  points  on  the 
library  film.  Reproduced  with  permis¬ 
sion  from  M.  Watanabe  et  al.,  J.  Comb. 
Chem.  10,  175  (2008).  Copyright  2008 
American  Chemical  Society. 
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sintering  of  powder  layers  of  different  compositions  to  form 
a  bulk  sample  with  a  compositional  gradient.  This  bulk 
library  sample  was  then  sliced  into  thin  wafers  to  expose  the 
compositional  gradient.  In  the  second  method,  discrete  com¬ 
positions  were  formed  in  microwells  in  a  graphite  sample 
holder  by  melting  and  slicing  the  pellets.  The  thermoelectric 
properties  of  the  Ni!_xCux  alloys  were  measured  from 
323  K  to  950  K,  and  the  results  are  shown  in  Fig.  62.  In 
particular,  Ni0.7Cu0.3  is  found  to  possess  a  large  power  fac¬ 
tor,  0.012  Wm_1K“2  at  around  950  K  (calculable  from 
Figs.  62(a)  and  62(b). 

In  addition  to  the  power  factor  figure  of  merit,  S 2cr, 
whose  components  are  reasonably  easy  to  measure  (even 
combinatorially),  the  more  complete  thermoelectric  power 
factor,  ZT  (equal  to  S2cfT/k)  is  often  desired.  Thus,  the  ther¬ 
mal  conductivity,  /c,  must  be  determined.  Using  Ni-based 
bulk  diffusion  multiple  sample  libraries,  Fig.  63,  thermal 
conductivity  of  Ni  alloys  were  determined  combinatorially 
by  time  domain  thermoreflectance  (TDTR)  measurements, 
with  a  spatial  resolution  of  about  four  micrometers,  and  a 
sampling  rate  of  3  points/s.235  The  authors  concluded  that 
their  measured  thermal  and  electrical  conductivities  for  the 
Ni  alloys  could  not  be  well  described  by  theory,  as  their 
measured  and  calculated  values  deviated  by  as  much  as  1 1  % 
to  29%,  depending  upon  the  alloy.  Thus,  their  combinatorial 
experiments  were  absolutely  necessary  as  a  consistency 
check  on  theoretical  values.  It  is  also  important  to  note  that 
their  data,  Fig.  64,  show  excellent  consistency  with  non- 
combinatorially  determined  values  for  certain  alloys. 

High  throughput  TDTR  measurements  have  also  been 
extended  to  oxide  thin  films.230  In  this  case,  thermal 


effusivity,  b ,  equal  to  ( Kcp)11 2  (where  k  is  the  thermal  con¬ 
ductivity,  c  is  the  specific  heat,  and  p  is  the  density),  was 
measured  on  a  metal  coated  oxide  film,  with  a  ten  microme¬ 
ter  spot  size.  The  authors  were  able  to  obtain  a  thermal  effu¬ 
sivity  calibration  curve  using  various  bulk  samples,  i.e., 
Si02,  SrTi03,  LaA103,  A1203,  and  Si.  Using  this  calibration, 
the  thermal  conductivity  of  an  yttria  (Y203)  containing  ter¬ 
nary  oxide  film  was  then  determined.  The  measured  value, 
11.96  J  (Kms)-1,  was  within  about  10%  of  the  reported  value 
of  12.87  J(Kms)-1.  This  tool  will  clearly  be  a  very  useful 
thermoelectric  (as  well  as  other  materials)  screening  tool. 

B.  Lithium  battery  materials 

Battery  powered  high  technology  products  are  ubiqui¬ 
tous.  Since  the  advent  of  rechargeable  nickel  metal  hydride 
(NiMH)  batteries  enabled  environmentally  friendly  portable 
electronics,  the  market  for  rechargeable  batteries  has  dramati¬ 
cally  expanded  236,237  Li  ion  batteries  have  emerged  as  the 
preferred  power  source  for  portable  electronics,  laptops,  cell 
phones,  and  high  performance  electric  vehicles  due  to  their 
superior  gravimetric  and  volumetric  energy  densities,  reduced 
self-discharge,  improved  cycling  capabilities,  and  thermal  sta¬ 
bility.  State  of  the  art  Li  ion  batteries  exhibit  remarkable  char¬ 
acteristics,  with  the  following  properties  (maxima)  reported: 
energy  storage  capacities  of  2400  mA-h;  specific  energies  of 
250W-h/kg;  energy  densities  of  620W-h/l;  and  specific 
powers  of  340W/kg.  Further,  they  can  be  cycled  between 
four  and  twelve  hundred  times.  Typical  cells  are  composed  of 
a  graphite  anode  that  intercalates  Li  ions  given  off  by  a 
LiCo02  cathode  during  discharge. 


FIG.  62.  The  temperature  dependence  of 
resistivity  (a)  and  Seebeck  coefficient 
(b)  of  Ni-Cu  alloys.  Reproduced  with 
permission  from  A.  Yamamoto  et  al., 
Thermoelectric  Power  Generation, 
273-278  (2008).  Copyright  2008 

Cambridge  University  Press. 
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FIG.  63.  Schematic  diagram  of  the  Ni-Cr-Pd-Pt-Rh-Ru  diffusion  multiple. 
The  five  thick  bars  mark  locations  where  the  composition  and  thermal  con¬ 
ductivity  profiles  were  acquired.  Reproduced  with  permission  from  X. 
Zheng  et  al.,  Acta  Mater.  55,  5177  (2007).  Copyright  2007  Pergamon. 

The  continued  demand  for  enhanced  battery  perform¬ 
ance  has  led  to  the  pursuit  of  improved  anode  and  cathode 
materials.  The  primary  reasons  for  modifying  the  cathode 
are  to  increase  cathode  capacity,  address  the  LiCo02  toxicity 
issue,  and  replace  Co  with  a  less  expensive  metal.  Much  of 
the  combinatorial  work  to  date  has  focused  on  materials  such 
as  LiMn03  and  Li-Ni-Co-Ti  oxides.  For  the  anode,  the  pri¬ 
mary  reason  is  the  desire  to  intercalate  more  Li  atoms  per  C 
(graphite)  atom,  thus  increasing  the  total  capacity  of  the  bat¬ 
tery,  while  reducing  the  degradation  of  battery  performance 
as  a  function  of  charge/discharge  cycles  (“cycle  fade”). 
Whereas  6  C  atoms  are  required  to  intercalate  one  Li  atom, 
the  same  can  be  accomplished  with  only  4  Sn  atoms.238  The 
primary  source  of  cycle  fade  for  Li  ion  battery  anodes  is  the 
large  volume  change  that  must  be  accommodated  during  Li 
intercalation,  which  introduces  large  stresses  and  therefore 
fracture  of  the  anode,  with  a  resultant  capacity  reduction. 
Thus,  two  major  strategies  have  been  pursued  to  preserve  the 
capacity  while  mitigating  cycle  fade,  i.e.,  placing  the  active 
metals  in  a  passive  matrix,  and  the  use  of  amorphous 
materials. 

The  large  compositional  phase  space  represented  by 
anodic  and  cathodic  materials  for  Li  ion  batteries  has  moti¬ 
vated  a  number  of  academic  and  industrial  groups  to  utilize 
high  throughput  methods  to  identify  useful  materials. 
Typical  cathode  libraries  are  prepared  by  UHV  sputter¬ 
ing238,239  and  ink  jet  deposition.240  Other  studies  employed 
more  traditional  bulk  synthesis  techniques  such  as  liquid 
phase  precursor  deposition,  where  libraries  are  created 
utilizing  automated  robots.241  Anode  libraries  tend  to  be  de¬ 
posited  via  sputtering  to  take  advantage  of  the  high  quench¬ 
ing  rates  inherent  to  that  technique.  Screening  of  both  anode 
and  cathode  libraries  is  typically  done  either  through  paral¬ 
lel  techniques  such  as  imaging,240  Fig.  65(a),  the  use  of 
automated  serial  arrays,242-245  Fig.  65(b),  or  the  use  of  serial 
electrochemical  tests  via  automated  scanning  probe 

246  247 

measurements. 
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FIG.  64.  (a)  Composition  dependence  of  the  thermal  conductivity  at  room 
temperature  for  Ni(Cr),  Ni(Pd),  Ni(Pt),  Ni(Rh),  and  Ni(Ru)  solid  solutions, 
(b)  Composition  dependence  of  the  thermal  conductivity  for  solute  concen¬ 
trations  between  0.1  and  10  at.  %  plotted  on  a  log-log  scale.  Open  triangles, 
solid  diamonds,  open  squares,  dotted  triangles,  and  solid  circles  are  data  for 
solutes  of  Pd,  Pt,  Rh,  Ru,  and  Cr,  respectively.  Open  circles  are  for  earlier 
non-combinatorial  data  reported  for  Ni(Cr)  alloys  (see  Ref.  436). 
Reproduced  with  permission  from  X.  Zheng  et  al.,  Acta  Mater.  55,  5177 
(2007).  Copyright  2007  Pergamon. 


1.  Cathode  materials 

LiNi02  is  considered  a  potential  replacement  material 
for  LiCo02  due  to  its  increased  capacity  and  the  lower  cost 
of  Ni.  However,  during  cycling,  LiNi02  undergoes  structural 
rearrangements,  which  eventually  lead  to  cycle  fade.  It 
has  been  proposed  that  the  cycling  capabilities  could  be 
improved  by  the  substitution  of  other  elements  for  Ni. 
Combinatorial  computational  chemistry  calculations  were 
used  to  investigate  the  structural  properties  of  LiCo02, 
LiNi02,  and  Lio.66Nio.5Meo.5O2,  where  Me  =  Cr,  Al,  and 
Co.248  To  investigate  the  de-intercalation  of  Li  from  the 
starting  materials,  structures  of  the  form  LixM02  with  x  =  1, 
0.66,  0.33,  and  0  were  created  by  removing  Li  ions  and 
relaxing  the  structure.  For  the  case  of  LiNi02,  as  the  Li  was 
de-intercalated,  the  Ni  ions  went  from  a  trivalent  to  a  tetrava- 
lent  state,  with  a  substantial  change  in  the  Ni-0  bond  length. 
It  was  found  that  with  a  1 : 1  ratio  of  Ni  and  Co,  the  Co-0 
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FIG.  65.  Example  of  use  of  (a)  fluorescence  libraries.  Reproduced  with  per¬ 
mission  from  E.  Reddington  et  al.,  Science  280,  1735  (1998).  Copyright 
1998  American  Association  for  the  Advancement  of  Science,  (b)  Apparatus 
for  parallel  electrochemical  measurements  to  screen  battery  materials. 
Reproduced  with  permission  from  M.  D.  Fleischauer  et  al.,  J.  Electrochem. 
Soc.  150,  A1465  (2003).  Copyright  2003  The  Electrochemical  Society,  Inc. 


bond  lengths  were  between  those  of  the  two  types  of  Ni 
bonds,  and  it  was  thus  postulated  that  this  could  mitigate 
the  distortion  of  the  LiNi02  structure  during  charging/ 
discharging. 

The  experimental  synthesis  of  Li-0  based  cathode  libra¬ 
ries  was  realized  through  the  use  of  an  automated  solution 
atomization  process,  called  the  “M-ist”  approach.  Using  an 
array  of  syringe  pumps,  liquid  precursors  were  dosed  into 
and  mixed  in  a  manifold,  delivered  to  a  robotically  posi¬ 
tioned  stainless  steel  nozzle,  and  then  atomized  onto  a 
substrate  biased  with  a  high  voltage.  With  this  approach  hun¬ 
dreds  of  samples  with  diameters  on  the  order  of  3  mm  were 
deposited;  the  total  amount  of  material  in  each  sample  was 
about  0.1  mg.241  It  has  been  proposed  that  layered  oxides,  or 
oxides  with  the  ramsdellite  structure  would  have  superior 
electrode  properties;  therefore,  initial  screening  of  the  mixed 
oxide  materials  libraries  focused  on  scanning  XRD  studies. 
In  the  Li-Ni-Co  ternary  oxide  system,  a  large  region  in  the 
center  of  the  phase  diagram  was  found  to  consist  of  single 
phase  layered  oxides.  In  the  pseudo  four  component  system 
Li-Ni-Co-Ti-O,  a  large  region  of  potential  cathode  materials 
with  layered  structures  of  the  form  Lix(NiyCozTiw)02,  cover¬ 
ing  the  range  of  compositions  (x^l,  0  <  y  <  1 ,  0<z<l, 
0<w<0.2)  was  found.249  A  large  region  containing  the 
ramsdellite  structure  was  also  found  by  substituting  CrOL5 
for  Ni  in  the  LiO0.5-X-TiO2  system.250 

More  than  25  000  potential  cathode  material  compounds 
were  explored  by  a  combination  of  high  throughput  model¬ 
ing  and  data  mining.251  As  a  result,  it  was  suggested  that 
oxide  materials  such  as  LiZn02  and  LiA102  (having  the  hy¬ 
pothetical  a-NaFe02  structure)  could  exhibit  higher  interca¬ 
lation  voltages  than  LiCo02.252  It  was  also  found  in  this 
work  that  phosphate  materials  generally  exhibit  larger  Li 
insertion  voltages  and  comparable  capacities  compared  to 
their  oxide  counterparts,  Fig.  66.  Through  further  computa¬ 
tional  and  experimental  methodologies,  it  was  shown  that 
exchanging  Co  with  Al  results  in  a  cathode  material  with  an 
increased  cell  voltage,  while  reducing  its  overall  density  and 


J  Phosphates 
•  Known  oxides 


5 


4.5  V 


0 

CD 

(C 

I 


4 

3 

2 


FIG.  66.  Comparison  of  calculated  range 
•  over  which  different  cations  can  be 

oxidized  by  Li  removal  to  voltages  of 
known  oxides.  Reproduced  with  permis¬ 
sion  from  G.  Ceder,  MRS  Bull.  35, 
693  (2010).  Copyright  2010  Materials 
Research  Society. 
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cost.253  The  key  computational  discovery  that  led  to  this 
breakthrough  was  the  finding  that  the  oxygen  atoms,  not 
the  metals,  were  responsible  for  accepting  the  electrons  dur¬ 
ing  Li  intercalation.  A  similar  approach  was  used  to  theoreti¬ 
cally  predict  and  experimentally  verify  that  Li(Ni0.5Mn0.5)O2 
cathodes  had  higher  rate  capability  than  LiCo02.254 
Further,  using  liquid  phase  reactions  and  post  synthesis  gra¬ 
dient  annealing  of  arrays  of  cathode  materials,  it  was  demon¬ 
strated  that  LiC0.33Ni.33Mn  3302  could  provide  competitive 
electrochemical  performance  with  standard  cathode 
materials.255,256 

2.  Anode  materials 

In  order  to  increase  the  resistance  of  anodes  to  cycle 
fade,  a  great  deal  of  research  has  been  devoted  to  amorphous 
alloys  containing  Si  and  Sn.  The  quest  to  identify  amorphous 
regions  of  candidate  phase  diagrams  has  led  to  the  almost 
exclusive  use  of  sputtering  to  deposit  libraries,  since 
extremely  high  quench  rates  can  be  achieved.  However, 
one  difficulty  of  the  sputtering  process  is  subtle  run-to-run 
compositional  variations  in  the  libraries,  which  can  add  com¬ 
plexity  to  trends  observed  when  correlating  multiple  meas¬ 
urements  taken  on  a  series  of  library  samples.  A  novel 
orbiting  sputtering  system  that  allows  for  the  simultaneous 
deposition  of  up  to  five  different  samples  in  a  single  deposi¬ 
tion,242,244  Fig.  67  has  been  developed  to  circumvent  this 
issue.  Thus,  in  a  single  deposition,  libraries  for  XRD,  SEM, 
gravimetric  analysis  (accomplished  by  depositing  the  materi¬ 
als  directly  onto  an  array  of  25  pre- weighed  Cu  disks),  and 
electrochemical  cell  analysis  could  be  realized.  This  system 
was  used  to  demonstrate  the  presence  of  a  large  range  of 
amorphous  compositions  in  the  Si-Al-Sn  system,  for  compo¬ 
sitions  containing  more  than  ~66  at.  %  Si.257 

A  large  amorphous  phase  region  was  observed  in  Si-Al- 
Mn  libraries,244  as  had  been  observed  before  in  rapidly  cooled 
bulk  alloys.258,259  Within  this  region,  Sio.3gAlo.52Mno.10  was 
identified  as  exhibiting  the  best  resistance  to  cycle  fade,  a 


gravimetric  capacity  superior  to  that  of  graphite,  and  a  reversi¬ 
ble  capacity  (based  on  standard  Li  battery  cell  type  18  650) 
40%  higher  than  for  graphite  based  cells.  The  predicted  resist¬ 
ance  to  cycle  fade  was  confirmed  by  melt  spun  ribbon  (bulk 
material)  experiments260  and  was  attributed  to  the  formation 
of  a  supersaturated  solid  solution  of  Si  and  Mn  in  a-Al. 

A  design  of  experiment  approach  was  developed  to 
facilitate  discovery  and  optimization  of  useful  TM  additives 
to  Sn-C  anodes,  which  led  to  the  selection  of  the  Sn-Co-C 
system,261,262  Fig.  68.  The  experimental  design  first  involved 
only  selecting  TMs  that  were  known  to  form  amorphous 
compounds  with  Sn.  A  second  criterion  was  the  lack  of  for¬ 
mation  of  TM  carbides.  XRD  showed,  for  instance,  that 
although  Sn-Co,  Sn-Ti,  and  Sn-V  are  all  amorphous  for 
Sn/TM  ratios  of  1,  the  addition  of  small  amounts  of  Carbon 
results  in  precipitation  of  the  TM  carbide.  It  was  inferred 
that  the  presence  of  nanoparticle  carbides  in  the  Ti  and  V 
libraries  facilitated  the  precipitation  of  Sn  out  of  the  amor¬ 
phous  matrix;  while  in  the  case  of  Co  additives,  no  Sn  pre¬ 
cipitation  was  observed;  further,  for  the  Co  case,  minimal 
cycle  fade  was  observed  up  to  27  cycles.  At  least  two  groups, 
working  with  bulk  materials,  have  confirmed  the  high 
cycling  stability  of  the  Sn-Co-C  alloys.263,264  This  work  has 
recently  been  extended  to  include  the  Sn-Si-C  system,  where 
the  C  serves  to  destabilize  the  two-phase  regions,  thus  reduc¬ 
ing  the  effect  of  volume  expansion  while  improving 
cycleability.265 

Sn-Cu  compounds,  particularly  Cu6Sn5,  have  been  pro¬ 
posed  as  potential  crystalline  anodes,  but  they  exhibit  signifi¬ 
cant  cycle  fade  within  the  first  50  cycles,266  thought  to  be 
due  to  the  large  stress  associated  with  the  volume  during  Li 
intercalation.  It  has  been  proposed  that  nano-sized  materials 
would  be  better  suited  to  accommodating  large  volume 
changes.267  In  a  separate  experiment,  C  particles  that  could 
mitigate  grain  growth  were  added,  and  a  capacity  retention 
in  excess  of  90%  was  observed.268  The  explanation  for  this 
behavior  is  the  formation  of  the  Cu6Sn5  phase,  which 
becomes  the  principal  Li  intercalation  material  within  a  C 
matrix,  and  prevents  grain  aggregation. 
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FIG.  67.  Schematic  of  orbital  sputtering  system.  Reproduced  with  permis¬ 
sion  from  M.  D.  Fleischauer  and  J.  R.  Dahn,  J.  Electrochem.  Soc.  151, 
A1216  (2004).  Copyright  2004  The  Electrochemical  Society. 


FIG.  68.  Roadmap  to  the  selection  of  the  Sn-Co-C  system  as  the  best  choice 
among  Sn-M-C  systems  for  negative  electrodes  in  Li-ion  batteries. 
Reproduced  with  permission  from  A.  D.  W.  Todd  et  al.,  J.  Electrochem. 
Soc.  154,  A597  (2007).  Copyright  2007  The  Electrochemical  Society. 
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C.  Hydrogen  storage  materials 

The  penetration  of  fuel  cell  powered  vehicles  into  the 
marketplace  is  contingent  upon  the  ability  to  store  substantial 

^fx\  r)ric\ 

amounts  of  hydrogen  within  existing  chassis  designs. 
Although  hydrogen  may  be  stored  either  under  high  pressure 
as  a  gas,  cryogenically  as  a  liquid,  or  in  the  solid  state  as  a 
hydride,  the  latter  option  offers  the  best  theoretical  combina¬ 
tion  of  gravimetric  and  volumetric  capacity,  stability,  and 
safety.  Due  to  its  low  density,  hydrogen  gas  must  be  stored 
at  pressures  as  high  as  70  MPa  for  a  three  hundred  mile  range 
vehicle.  Cryo-compressed  hydrogen  systems  can  provide 
sufficiently  large  gravimetric  and  volumetric  densities  of 
hydrogen,  although  substantial  losses  of  the  liquid  fuel  due 
to  imperfect  insulation  will  occur.  On  the  other  hand,  hydro¬ 
gen  can  be  stored  as  the  hydride  of  light  elements,  such  as 
Li,  B,  Mg,  N,  Al,  and  Na,  and  hydrogen  concentrations  as 
high  as  18  wt.  %  may  be  achieved  in  some  hydrides.271  The 
volumetric  density  of  gaseous  hydrogen  at  RT  and  1  bar  is 
8.2  x  10-5  kg/1,  while  metal  hydrides  such  as  LiBH4  can 
store  as  much  as  1.2  x  10“ 1  kg/1.  Even  at  70  MPa,  the  volu¬ 
metric  density  of  gaseous  hydrogen  is  5.7  x  10-2  kg/1, 
approximately  half  that  of  LiBH4.  The  efficient  renewable 
production  of  sufficiently  pure  hydrogen  represents  a  signifi¬ 
cant  challenge  that  must  be  addressed  in  conjunction  with 
the  search  for  suitable  hydride  materials.  High-throughput 
efforts  in  the  field  of  solar  energy  splitting  of  water  have 
been  developed  by  several  groups,272-274  and  a  recent  critical 
review  on  this  subject  has  been  published.275 

Despite  these  advantages,  the  implementation  of  metal 
and  chemical  hydrides  has  been  hindered  by  the  dehydrogen¬ 
ation  thermodynamics  and  kinetics  of  typical  materials.  A 
large  portion  of  the  initial  effort  of  the  US  Department  of 
Energy  (DOE)  sponsored  metal  hydride  work  focused  on 
studying  known  hydrides,  with  or  without  the  addition  of  TM 
based  oxides  and  chlorides,  to  expedite  kinetics.276  First,  the 
addition  of  TiCl3  to  NaAlH4  had  been  shown  to  substantially 
change  its  thermodynamics  and  kinetics  and  thus  an  entire 
community  formed  around  the  deliberate  doping  of  known 
hydrides.277,278  Second,  the  dehydrogenation  mechanisms  of 
the  existing  materials  had  not  been  well  studied  and  extensive 
work  was  needed  to  gain  a  basic  understanding  of  the  chemi¬ 
cal  mechanisms  involved.  Finally,  it  was  thought  that  mix¬ 
tures  of  thermodynamically  stable  complex  hydrides  would 
not  form  single  phases,  but  rather  multiphase  mixtures. 

The  landscape  of  the  hydrogen  storage  community  was 
significantly  altered  by  the  reassessment  of  some  early 
research279  that  demonstrated  the  enthalpy  of  formation  for 
MgH2  could  be  lowered  from  76kJ/mol  H2  to  63kJ/mol  H2 
by  alloying  with  Ni.  Thus,  a  niche  for  combinatorial  materi¬ 
als  science  was  formed  to  discover  rapid  methods  of  screen¬ 
ing  additives  to  metal  hydrides  based  on  their  efficacy  in 
promoting  hydrogenation/dehydrogenation  kinetics.  More 
recently,  it  was  shown  that  a  1:1  mixture  of  LiBH4  and 
MgH2  resulted  in  a  substantial  reduction  in  the  enthalpy  of 
formation  of  pure  LiBH4,  from  90kJ/mol  H2  to  50kJ/mol 
H2  280,281  EXperiments  involving  mixing  of  other  complex 
hydrides,  such  as  LiNH2  with  MgH2  or  LiH,  or  LiBH4  and 
LiNH2  soon  followed.282-286  Surprisingly,  in  1 : 1  mixtures  of 


LiNH2  and  MgH2,  it  has  since  been  demonstrated  that  after 
only  one  dehydrogenation/hydrogenation  cycle,  no  LiNH2 
remains  in  the  powder;  subsequent  hydrogen  cycling 
involves  only  the  compounds  Mg(NH2)2,  LiH,  and  MgH2.287 
Further,  it  has  been  shown  experimentally  that  the  initial 
ratio  of  LiNH2  to  MgH2  can  affect  the  overall  reaction  mech¬ 
anism.288,289  For  the  case  of  LiBH4  and  LiNH2  mixtures,  it 
was  surprisingly  discovered  that  almost  complete  solubility 
between  the  two  phases  could  be  achieved,  with  free  inter¬ 
change  of  the  NH2  and  BH4  anions.290  Investigations  into 
ternary,  quaternary,  and  higher  order  mixtures  were  found  to 
yield  still  more  interesting  materials  with  previously  unin¬ 
vestigated  chemistries.291-293 

The  multidimensional  compositional  space  of  importance 
to  the  hydrogen  storage  community,  encompassing  cation  (Li, 
Mg,  Na,  Ti,  K,  etc.)  and  anion  substitutions  (A1H4,  BH4,  NH2, 
etc.),  catalyst  composition  (MgH2,  TiCl3,  etc.),  nanostructure, 
thermodynamics,  and  kinetics  is  inaccessible,  from  a  practical 
standpoint,  using  traditional  materials  science  techniques. 
Therefore,  the  hydrogen  storage  community  has  embraced 
high  throughput  methodologies  to  enhance  the  probability  of 
discovering  breakthrough  materials.  Their  efforts  have  spanned 
the  spectrum  from  database  searches  and  first  principles  DFT 
calculations,  to  the  demonstration  of  high  throughput  thin  film 
measurement  techniques;  the  latter  are  especially  useful  for  ex¬ 
perimental  verification  of  predictive  modeling  results. 

What  follows  in  this  section  is  divided  into  three  subsec¬ 
tions:  data  mining  and  modeling,  used  to  rapidly  screen  mil¬ 
lions  of  possible  reactions  for  the  relatively  small  number 
that  are  theoretically  possible;  thin  film  combinatorial 
experiments,  aimed  at  developing  a  comprehensive  high 
throughput  scheme  that  enable  working  with  nanogram  to 
microgram  quantities  of  material;  and  combinatorial  powder 
studies,  which  mimic  the  bulk  materials  that  would  actually 
be  used  in  a  hydrogen  storage  system. 

1.  Data  mining  and  modeling 

The  advent  and  general  availability  of  DFT  software 
packages  such  as  Vienna  Ab  Initio  Simulation  (vasp),294  as 
well  as  large  scale  tabulations  of  experimentally  obtained 
crystal  structure  databases,295-297  have  enabled  computa¬ 
tional  screening  of  millions  of  possible  hydrogenation/dehy¬ 
drogenation  reactions.  These  methods  provide  a  powerful 
means  of  searching  an  enormous  phase  space  in  a  short 
(days)  time  span.  For  example,  more  than  one  hundred  dehy¬ 
drogenation  reactions  involving  49  compounds  were  compu¬ 
tationally  screened,298  Fig.  69.  From  these  reactions,  five 
were  identified  with  reaction  enthalpies  that  would  allow  the 
release  of  hydrogen  at  pressures  and  temperatures  relevant  to 
vehicular  applications.  Through  detailed  comparison  of  their 
calculated  reaction  enthalpies  and  values  available  through 
the  literature,  the  authors  determined  that  their  error  bar  was 
on  the  order  of  ±10kJ/mole  H2;299  however,  depending 
upon  the  system  chosen,  this  error  can  correspond  to  a  varia¬ 
tion  in  predicted  equilibrium  pressure  of  two  orders  of  mag¬ 
nitude.  An  extremely  comprehensive  screening  of  16  x  106 
reactions300,301  was  also  carried  out;  it  is  worth  noting  that  in 
this  case,  only  44  candidate  reactions  were  identified. 
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FIG.  69.  DFT  calculated  enthalpies  of  formation  for  a  large  number  of  possi¬ 
ble  one  step  hydrogen  release  reactions.  The  reactions  that  fall  within  the 
bounds  of  the  dotted  lines  were  identified  as  being  promising  enough  to 
merit  experimental  investigation.  Reproduced  with  permission  from  S.  V. 
Alapati  et  al.,  J.  Phys.  Chem.  C  111,  1584  (2007).  Copyright  2007  American 
Chemical  Society. 

A  similar  methodology  was  used  to  investigate  destabili¬ 
zation  of  LiBH4  and  Ca(BH4)2  by  TMs  and  simple  metal 
hydrides.280  Twenty  candidate  reactions  were  investigated, 
and  it  was  found  that  additives  such  as  Cr  and  ScH2  could 
sufficiently  destabilize  Ca(BH4)2  to  the  temperature  and 
pressure  requirements  for  vehicular  applications.  Further, 
prior  to  this  work,  it  was  assumed  that  borohydrides  such  as 
LiBH4  decomposed  into  a  mixture  of  Li,  B,  and  H2. 
Subsequent  work  demonstrated  that  an  entire  series  of  inter¬ 
mediate  compounds,  such  as  Li2B12H12  and  CaB12H12  are 
formed  during  hydrogen  cycling  of  many  borohy- 

ary)  qnq 

drides.  ’  The  effect  of  destabilizers  on  the  reaction  paths 


of  these  newly  discovered  reaction  products  were  also  stud¬ 
ied,304  where  it  was  found  that  the  intermediate  phases  were 
very  stable,  with  enthalpies  greater  than  70kJ/mol  H2,  but 
could  be  substantially  destabilized  in  the  presence  of  basic 
metal  hydrides  such  as  MgH2.  In  order  to  perform  DFT  cal¬ 
culations  for  M(B12H12)  materials,  their  structure  must  be 
calculated.305 

Further  work  in  the  simulation  and  modeling  area 
was  directed  towards  reaction  pathways,  enthalpies,  and 
hydrogen  release  in  the  MgH2-LiH-NH3  ternary  system, 
Fig.  70. 306  An  outcome  of  another  simulation  was  the  first 
report  of  some  phase  structures  in  the  Li-Mg-Al-H  system, 
e.g.,  the  unstable  CaGaF5-type  structure  of  LiMg(AlH4)3; 
here,  it  was  determined  that  LiMg(AlH4)3  decomposes  to 
LiMgAlH6  at  RT.307  Recent  interest  in  ammonia  borane 
(NH3BH3)  and  its  intermediate  decomposition  products 
spurred  calculations  of  the  structure  and  reaction  enthalpies 
of  (NH4)2B12H12.308 

2.  Thin  film  combinatorial  work 

The  goal  of  any  hydrogen  storage  program  is  to  manu¬ 
facture  materials  on  the  kilogram  scale.  However,  a  number 
of  groups  have  utilized  high  throughput  thin  film  methodolo¬ 
gies  to  screen  hydrogen  storage  materials  on  the  microgram 
scale308-311  and  have  been  able  to  determine  the  phase  dia¬ 
gram  of  a  hydride  system  on  a  single  library  sample.  Also, 
thin  films  can  mitigate  the  role  of  kinetics  in  screening 
hydride  samples;  the  films  can  be  grown  with  arbitrarily 
small  grain  sizes,  allowing  for  rapid  cycling  of  materials,  or 
the  effective  grain  size  can  be  deliberately  modulated  to 
determine  the  presence  of  rate  limiting  steps  during  hydro¬ 
genation.  In  thin  film  experiments,  the  compositional  resolu¬ 
tion  is  generally  limited  by  the  size  of  the  sample  and  the 
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FIG.  70.  Calculated  hydrogen  release  (in  wt.  %)  across  the  Li-Mg-N-H  system  as  a  function  of  temperature.  Each  figure  corresponds  to  a  specific  temperature 
range  (a)  T  <  30  K,  (b)  1 30  K  <  T,  426  K,  (c)  426  K  <  T  <  458  K,  (d)  458  K  <  T  <  606  K,  (e)  606  K  <  T  <  734  K,  and  (f)  743  K  <  T  <  780  K.  Reproduced  with 
permission  from  C.  Wolverton  et  al.,  J.  Phys.  Condens.  Matter  20,  064228  (2008).  Copyright  2008  Institute  of  Physics. 
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spatial  resolution/sampling  volume  requirements  of  the  mea¬ 
surement  system.  The  number  of  samples  to  be  characterized 
from  a  given  library  can  vary  from  100  to  1000. 3 11-3 13 

However,  there  are  significant  challenges  to  the  thin  film 
approach.  For  example,  a  phase  that  can  be  synthesized 
through  thin  film  deposition  techniques  may  not  be  capable 
of  being  produced  in  kilogram  quantities  using  traditional 
bulk  synthesis  techniques.  Also,  in  thin  film  form,  the  influ¬ 
ences  of  crystallographic  texture  and  film  stress  may  alter 
the  thermodynamics  of  the  materials  being  investigated. 
Finally,  in  thin  films,  measurements  of  hydrogen  storage 
properties  cannot  be  obtained  using  traditional  volumetric 
and  gravimetric  measurement  techniques.  To  circumvent  this 
issue,  the  combinatorial  community  has  utilized  a  series  of 
indirect  measurement  techniques,  such  as  film  transparency, 
IR  emissivity,  XRD,  and  cantilever  deflection  to  characterize 
thin  film  libraries.310,314-316  Using  these  tools,  combinatorial 
thin  film  studies  can  provide  an  important  bridge  between 
theoretical  and  bulk  powder  work.  Thus,  thin  film  demon¬ 
stration  of  a  predicted  reaction  pathway,  with  an  enthalpy 
close  to  the  expected  value,  can  provide  evidence  that  such  a 
reaction  is  not  kinetically  limited.  In  the  bulk  powder  case, 
demonstrating  that  a  given  reaction  does  not  occur,  but  rather 
decomposition  proceeds  through  a  different  reaction  path 
can  help  theorists  refine  their  models  and  add  new  com¬ 
pounds  to  their  databases  for  subsequent  studies. 

Many  materials,  such  as  Mg  and  Y,  undergo  a  metallic 
to  insulator  transition  upon  exposure  to  hydrogen.317  The 
consequent  increase  in  resistivity  results  in  an  increase  in  the 
amount  of  IR  radiation  by  the  part  of  the  library  sample  that 
has  absorbed  hydrogen.  IR  emissivity,  Fig.  71,  was  used  to 
screen  for  hydrogen  absorption  temperatures  in  thin  film 
libraries  of  Mg-Ni-Fe,318  Fig.  72.  In  a  subsequent  study,  IR 
emissivity  mapping  was  used  in  tandem  with  TEM  to  study 
hydrogenation  in  Mg-Ni  libraries.319  They  found  that  over  a 
broad  range  of  Ni  concentrations,  a  Ni  stabilized  fee  Mg 


FIG.  71.  Cartoon  of  a  typical  IR  emissivity  measurement  system  for  thin 
film  hydride  characterization.  Reprinted  with  permission  from  C.  H.  Oik 
et  al.,  J.  Appl.  Phys.  94,  720  (2003).  Copyright  2003  American  Institute  of 
Physics. 


phase  coexisted  with  precipitates  of  Mg2Ni  and  MgNi2, 
which  resulted  in  an  increased  hydrogenation  pressure.  It 
was  further  demonstrated  that  this  technique  could  be  used 
to  monitor  metal  to  semiconductor,  as  well  as  metal  to  metal 
transformations;  this  finding  is  important  since  many  hydro¬ 
gen  storage  materials  can  exhibit  multiple  decomposition 
steps,  and  not  all  of  the  products  are  insulating. 

CCS  libraries  of  hydrogen  storage  materials  enable  mea¬ 
surement  of  the  thermodynamic  properties,  but  to  determine 
reaction  kinetics,  it  is  useful  to  employ  thickness  wedge 
libraries,  where  the  importance  of  grain  size  on  the  reaction 
rates  can  be  studied.  Materials  such  as  MgH2  strongly  inhibit 
the  diffusion  of  hydrogen,  and  thus  present  a  substantial  bar¬ 
rier  to  efficient  hydrogenation.  The  wedge  library  technique 
was  used  to  great  advantage321  to  demonstrate  that  for  a 
Mg98.4Ti16  film  the  MgH2  layer  grows  quickly  (0.3nm/s), 
until  it  is  about  150  nm  thick.  At  this  stage,  the  hydride 
growth  rate  drops  by  an  order  of  magnitude;  thus,  the  opti¬ 
mal  particle  size  for  the  bulk  powder  alloy  would  be  in  the 
range  of  300  nm,  since  hydride  formation  in  the  bulk  would 
progress  by  radial  hydrogen  diffusion.  Thus,  the  information 
obtained  from  thickness  gradient  libraries  can  provide  the 
manufacturers  of  the  bulk  powder  alloys  with  important  par¬ 
ticle  size  information. 

The  measurement  of  the  change  in  optical  transmissivity 
in  thin  film  libraries  due  to  hydrogen  reaction  has  been 
coined  “hydrogenography”  and  has  proved  to  be  a  robust 
screening  technique.  ’  In  a  two  phase  mixture  of  a  metal 
and  its  hydride,  the  concentration  of  hydrogen  is  proportional 
to  the  log  of  the  optical  transmission.323  Thus,  by  measuring 
the  change  in  transmission  as  a  function  of  pressure,  the 
enthalpy  of  hydride  formation  for  a  thin  film  can  be  directly 
measured.  Much  of  the  hydrogenography  work  has  been 
reviewed.324  Similar  to  IR  emissivity,  hydrogenography  can 
also  detect  metal  to  metal  transitions  during  hydrogenation 
of  films;325  this  feature  makes  hydrogenography  useful  for 
studying  other  classes  of  materials,  such  as  hydrogen  separa¬ 
tion  membranes.326  Hydrogenography  has  also  been  utilized 
as  a  switchable  mirror  technology,  taking  advantage  of  the 
change  in  transmissivity  that  results  from  hydrogenation,  as 
well  as  other  applications.327-330 

Hydrogenography  has  been  applied  for  compositional 
mapping  of  a  number  of  Mg  based  systems  such  as 
Mg-Ni,309,331  Mg-Ni-Ti,332  Mg-Al,333  and  Mg-Al-Ti.334  A 
notable  result  from  this  body  of  work  was  the  discovery316  of 
the  alloy  Mg69Ni26Ti5,  Fig.  73,  which  has  an  hydrogenation 
enthalpy  of  —  39.2kJ/mol  H2,  quite  close  to  the  ideal  value  of 
— 40kJ/mol  H2  for  on  board  vehicular  hydrogen  storage.  A 
subsequent  DFT  study  validated  the  experimentally  measured 
enthalpy  and  demonstrated  that  in  the  hydrogenated  state  the 
alloy  existed  as  a  Ti-doped  Mg2NiH4  alloy,  which  upon  dehy¬ 
drogenation  decomposes  to  a  mixture  of  Mg2Ni  and  a  Ti-Ni 
intermetallic  phase.  In  a  related  study,  it  was  found  that  using 
a  Ti  buffer  layer  in  Mg-Al  films  facilitated  the  hydrogenation 
of  Mg(AlH4)2  below  1  bar.  Unfortunately,  knowledge  of  such 
catalytic  behavior  has  not  been  able  to  inform  the  practice  of 
using  milling  to  dope  bulk  powder  lots. 

Other  approaches  for  investigating  thin  film  hydrogen 
storage  materials  libraries  include  XRD,  prompt  gamma 
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FIG.  72.  Representative  IR  emissivity 
data  from  a  Mg-Ni  composition  spread. 
Reprinted  with  permission  from  C.  H. 
Oik  et  al.,  J.  Appl.  Phys.  94,  720  (2003). 
Copyright  2003  American  Institute  of 
Physics. 


activation  analysis  (PGAA),  and  RBS. 314,335,336  Although 
these  methods  are  useful,  there  are  disadvantages  to  each.  In 
the  case  of  PGAA  and  RBS  large,  expensive  and  dedicated 
equipment  are  required.  XRD  has  limited  usefulness  because 
many  hydrogen  storage  materials,  e.g.,  borohydrides,  tend  to 
decompose  into  amorphous  phases.  Another  intriguing 
approach  to  combinatorial  screening  of  hydrides  is  the  mea¬ 
surement  of  microelectromechanical  (MEMS)  cantilever 


displacement,  which  results  from  the  expansion/contraction 
of  the  material  as  it  is  cycled.310,337  Thus,  far  this  technique 
has  only  been  applied  to  simple  Mg  hydrides. 

3.  Combinatorial  studies  of  hydride  powders 

The  goal  of  high  throughput  experimentation  on  powder 
samples  is  to  demonstrate  their  properties  on  the  gram  to 


FIG.  73.  Representative  data  taken  via 
hydrogenography  (a)  an  image  of  a 
backlit  combinatorial  sample  is  obtained 
at  a  set  temperature  and  pressure,  (b)  the 
transmission  of  each  point  is  mapped  to 
the  ternary  phase  diagram,  (c)  the  pla¬ 
teau  pressure  of  each  composition  is 
evaluated  as  a  function  of  temperature, 
and  (d)  the  enthalpy  is  plotted  on  the 
phase  diagram.  Reproduced  with  permis¬ 
sion  from  R.  Gremaud  et  al.,  Adv. 
Mater.  19,  2813  (2007).  Copyright  2007 
Wiley— VCH  Verlag  Gmbh  &  Co. 
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kilogram  scale.  In  high  throughput  ball  milling  and  charac¬ 
terization  experiments,  the  number  of  different  material  sam¬ 
ples  is  limited  to  the  number  of  milling  jars.  The  price  of 
such  systems  rapidly  increases  with  the  number  of  samples 
to  be  investigated  per  experiment.  Typical  bulk  powder  syn¬ 
thesis  and  processing  involves  physical  mixing,  via  ball  mill¬ 
ing,  of  base  materials  (e.g.,  LiBH4),  with  a  destabilizing 
material  to  alter  its  thermodynamics  (e.g.,  MgH2),  and  finally 
a  catalyst  to  promote  kinetics  (e.g.,  TiCl3).  The  ball  to  pow¬ 
der  ratio,  milling  speed  (energy),  ball  composition,  and  order 
in  which  the  components  are  milled  can  dramatically  affect 
their  hydrogen  storage  properties.  For  example,  in  the  ter¬ 
nary  LiBH4-LiNH2-MgH2  system,  the  overall  hydrogen  reac¬ 
tion  pathways  can  be  altered  depending  upon  whether  all 
components  are  milled  simultaneously  in  a  single  jar,  or  are 
first  milled  individually  to  achieve  smaller  starting  grain 

9Q1  ?Q^ 

sizes  and  then  milled  together. 

Several  industrial  scale  high  throughput  efforts  have  been 
implemented  to  screen  for  powder  materials,  using  parallel 
ball  milling  apparatus.338,339  These  systems  can  synthesize 
and  measure  the  pressure-composition-isotherm  characteris¬ 
tics  of  up  to  forty-eight  individual  powder  samples  at  a  time. 
In  a  typical  experiment,  a  sample  of  known  mass  is  held  in  a 
vessel  of  known  volume  at  a  constant  temperature.  As  the 
pressure  is  increased  in  discrete  increments,  the  change  in  the 
system  pressure  due  to  hydrogen  absorption/desorption  is 
measured.  Using  this  methodology,  a  self-catalyzing  hydride 
system,  LiNH2-LiBH4-MgH2,338  was  discovered,  Fig.  74. 
The  optimal  combination  of  low  temperature  hydrogen 
release  (<180  °C),  weight  percent  hydrogen  released  (less 
than  4wt.  %),  and  reversibility  were  observed  at  ratios  of 
2:1:1  and  2:0.5: 1  in  this  system.340  An  interesting  aspect 
of  work  on  this  ternary  system  was  the  comprehensive  use  of 
both  experiment  and  theory  to  first  identify  promising  sys¬ 
tems,  followed  by  larger  scale  powder  experiments. 

For  most  hydrogen  containing  compounds,  vibrational 
spectroscopy  can  provide  important  structural  and  local 
binding  site  information  about  the  material;  such  informa¬ 
tion  is  complementary  to  gravimetric  or  volumetric  meas¬ 
urements.  In- situ  high  throughput  spectroscopy  has 
increasingly  become  a  part  of  the  array  of  tools  used  by  the 
hydrogen  storage  community.312,341,342  It  has  been  applied 
to  discrete  powder  arrays  as  well  as  to  diffusion  couple  sam¬ 
ples.  Micro-Raman  spectroscopy  was  used  to  measure  the 
diffusion  of  atomic  and  molecular  hydrogen  in  the  borohy- 
dride  cluster,341,342  and  it  was  demonstrated,  by  radioiso¬ 
topic  substitution,  that  hydrogen  was  transported  almost 
exclusively  by  the  diffusion  of  BH4-,  the  exchange  of  single 


hydrogen  ions  occurring  at  a  rate  ten  orders  of  magnitude 
slower. 

D.  Catalysts  for  fuel  cell  anodes 

Fuel  cells,  which  convert  chemical  energy  into  electrical 
energy,  have  long  been  considered  a  potential  device  for 
both  mobile  and  stationary  energy  generation  applica¬ 
tions.343  In  the  most  basic  case,  hydrogen  and  oxygen  react 
to  form  water,  and  electrical  energy  is  generated.  However, 
fuel  cells  can  run  on  fuels  other  than  hydrogen,  as  is  the  case 
illustrated  in  Fig.  75,  in  which  a  mixture  of  methanol  and 
water  impinge  on  the  anode  of  the  cell,  and  are  electrochemi- 
cally  reformed  and  reduced,  respectively,  into  H2  gas  and 
C02.  The  proton  is  conducted  through  an  electrically  insulat¬ 
ing  layer,  while  the  electron  is  conducted  through  an  electri¬ 
cal  circuit,  creating  a  current.  Finally,  at  the  cathode, 
molecular  oxygen  is  reduced  and  combined  with  the  proton 
and  electron  to  form  water.  Carbonaceous  fuels  such  as 
methanol  (CH3OH)  can  serve  as  alternatives  to  molecular 
hydrogen;  these  fuels  will  result  in  the  formation  of  H2  and 
C02.  The  advantage  of  alternate  fuels  is  that  the  transporta¬ 
tion  infrastructure  has  a  century  of  experience  in  handling 
liquid  fuels,  and  virtually  none  with  hydrogen.  Additionally, 
many  sources  of  hydrogen  are  available,  from  methane 
(CH4)  to  more  complex  hydrocarbon  fuels.  Much  effort  has 
been  directed  to  finding  appropriate  catalysts  for  solid 
oxide  fuel  cells  (SOFCs)  and  direct  methanol  fuel  cells 
(DMFCs).24, 344-347  However,  oxidation  of  carbonaceous 
fuels  takes  place  in  multiple  steps,  generally  including  C-H 
bond  activation,  C-0  oxidation,  and  water  adsorption,  and 
no  single  element  catalyst  has  thus  far  shown  ideal  activity 
and  selectivity.  The  focus  of  a  large  portion  of  the  combina¬ 
torial  work  done  in  this  field  has  been  on  the  identification  of 
novel  cathodes  and  anodes  for  SOFC  and  DMFC  applica¬ 
tions;  however,  concerns  over  the  oxidation  of  carbon  sup¬ 
ports  in  polymer  electrolyte  membrane  (PEM)  fuel  cells  has 
led  to  some  combinatorial  work  in  that  field.348 

Several  high  throughput  approaches  to  fuel  cell  anode 
catalyst  discovery  and  optimization  have  been  pursued.  One 
method  uses  sputtering  to  create  CCS  libraries.349-352  Inkjet 
based  solution  synthesis  techniques,  in  which  metal  salts  are 
solution  deposited  onto  an  electrode,  have  also  been 
employed.353-355  The  arrays  are  then  reduced  to  the  metallic 
state,  frequently  by  borohydrides  such  as  NaBH4.  One  con¬ 
cern  with  the  uniform  usage  of  NaBH4  as  a  reducing  agent  is 
its  inability  to  provide  the  high  surface  area  catalysts  eventu¬ 
ally  required  for  implementation.355  Another  technique  is 


FIG.  74.  Weight%  hydrogen  desorbed  from  differ¬ 
ent  mixtures  of  LiNH2-LiBH4-MgH2  as  a  function 
of  cycle  and  temperature,  (a)  First  cycle  at  220  °C, 
(b)  second  cycle  at  285  °C,  and  (c)  fourth  cycle 
350  °C.  Reproduced  with  permission  from  G.  J. 
Lewis  et  al.,  J.  Alloys  Compd.  446,  355  (2007). 
Copyright  2007  Elsevier  BY. 
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FIG.  75.  Schematic  of  a  single  fuel  cell  (see  Ref.  437).  Reproduced  with 
permission  from  S.  Eccarius  et  al.,  J.  Power  Sources  179,  723  (2008). 
Copyright  2008  Elsevier  S.A. 


wet  impregnation,  followed  by  calcination  at  high  tempera¬ 
tures  to  promote  oxidation  and  mixing.356,357  Other  groups 
employ  electrodeposition,  in  which  a  fluid  containing  the 
catalyst  precursors  is  pipetted  onto  an  array  of  electrodes, 
followed  by  the  electrodeposion  of  the  catalyst  on  the  elec¬ 
trode  surface  using  a  potentionstat,358  Fig.  76.  In  subsequent 
work,  gel  transfer  technology  was  used  in  combination  with 
electrodeposition  to  design  high  density  CCS  libraries.359  In 
this  method,  a  gel  containing  an  uniform  concentration  of  a 
single  metal  salt  is  impregnated  at  a  single  point  with  a  sec¬ 
ond  metal  salt;  the  diffusion  of  the  point  source  into  the 
semi-infinite  gel  then  creates  the  desired  composition  gradi¬ 
ent.  Samples  are  then  electrodeposited  using  normal  methods 
and  the  gel  is  removed,  yielding  samples  with  continuously 
varying  composition  gradients. 

Electrochemical  interrogation  of  the  catalyst  libraries  is 
accomplished  by  either  qualitative  screening  for  the  C-0  oxi¬ 
dation  potential  or  quantitative  measurements  of  catalytic 
activity  by  monitoring  the  current  as  a  function  of  voltage.  A 
popular  characterization  method  for  DMFC  anodes  is  fluores¬ 
cence.240  In  this  technique,  the  library  is  deposited  onto  a  con¬ 
ductive  substrate,  which  acts  as  a  common  electrode  for  each 
sample.  The  entire  sample  is  then  immersed  in  a  solution  of 
quinine  and  methanol,  and  the  voltage  is  slowly  raised.  When 
the  methanol  is  oxidized,  it  releases  up  to  4  protons,  which 
are  absorbed  by  the  quinine,  causing  it  to  fluoresce.  This  tech¬ 
nique  has  been  widely  used  as  a  primary  screen  for  anodes, 
although  it  suffers  from  some  disadvantages;  it  is  qualitative, 
and  offers  poor  spatial  resolution  due  to  diffusion  of  the  acti¬ 
vated  quinine  through  the  liquid. 167,353, 354, 360-362  mitigate 
diffusion  issues,  some  groups  have  created  semi-discrete 
arrays  by  micromachining.349  Another  method  of  circumvent¬ 
ing  the  limited  spatial  resolution  of  fluorescence  is  to  raster  a 
small  counter  electrode,  which  acts  as  a  scanning  electro- 
chemical  microscope,  ’  because  areas  of  the  library  with 
increased  activity  will  generate  more  hydrogen  and  thus  a 
larger  current.  Similar  electrochemical  screening  techniques 
have  been  employed  to  characterize  the  effect  of  corrosion  in¬ 
hibiting  additives  to  the  electrolyte  on  different  metallic 
wires,  to  great  effect.364"366 


FIG.  76.  Automated  system  for  electrosynthesis,  (a)  Rapid  serial  deposition 
system;  here,  a  moving  arm  is  used  to  sequentially  address  each  well,  (b) 
parallel  deposition  system;  here  an  array  of  counter  electrodes  are  used  to 
simultaneously  electrodeposit  each  well.  Reproduced  with  permission  from 
B.  Sung-Hyeon  et  al.,  Meas.  Sci.  Technol.  16,  54  (2005).  Copyright  2005 
Institute  of  Physics. 

A  more  quantitative  approach  employs  multielectrode 
fuel  cell  stacks  to  enable  the  direct  measurement  of  C-V  char- 

onri 

acteristics  of  the  catalyst.  “  Statistical  analysis  of  electro¬ 
chemical  experiments  run  on  a  40  cell  stack  showed  results 
that  were  consistent  with  single  cell  experiments.370  Further, 
in  this  research,  the  cell  stack  was  optimized  for  catalyst  load- 

QZQ  QHfA 

ing  and  the  concentration  of  methanol  in  the  feed. 
Another  study  used  a  well  characterized  array  membrane  as¬ 
sembly  fuel  cell,  capable  of  characterizing  up  to  twenty-five 
separate  electrodes,  in  combination  with  a  statistical  experi¬ 
mental  design  to  benchmark  several  potential  synthesis 
techniques  for  Pt-Ru  DMFC  anode  catalysts  against  commer¬ 
cially  available  catalysts,367,368’372,373  Fig.  77;  the  commer¬ 
cially  available  Pt-Ru  catalysts  prepared  by  Johnson  Matthey 
outperforming  the  other  synthesis  techniques 

Of  all  the  single  element  catalysts,  Pt  shows  the  highest 
activity  for  methanol  oxidation  at  temperatures  compatible 
with  fuel  cells  with  polymer  membranes.  However,  CO, 
which  forms  during  the  oxidation  of  methanol,  strongly 
adsorbs  to  Pt  surfaces,  inhibiting  usage  of  catalytically  active 
sites.  The  combination  of  Pt  with  an  oxophilic  material,  such 
as  Ru,  can  alleviate  this  problem.355  When  Ru  is  substituted 
into  Pt,  it  has  a  tendency  to  be  leached  from  the  surface  of 
the  catalyst  during  electrochemical  cycling.  To  circumvent 
this  issue,  the  ordered  intermetallic  phases  PtBi  and  PtPb 
have  recently  come  to  attention.  Using  fluorescence  imaging 
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FIG.  77.  Schematic  of  an  array  fuel  cell  with  common  counter  electrode  and 
array  working  electrode.  Reproduced  with  permission  from  E.  S.  Smotkin 
and  R.  R.  Diaz-Morales,  Ann.  Rev.  Mater.  Res.  33,  557  (2003).  Copyright 
2003  Annual  Reviews. 

and  SEM  to  screen  Pt-Bi-Pb  ternary  and  Pt-Ti-Co-Pb  quater¬ 
nary  combinatorial  thin  film  libraries  for  their  activity  to 
methanol  decomposition,349,360,374  it  was  found  that  by 
lightly  doping  a  PtPbo.53  alloy  with  Bi  the  onset  potential  of 
methanol  reduction  could  be  decreased  by  50  mV.  XRD  and 
XPS  analyses  showed  the  material  to  have  a  structure  similar 
to  PtPb.  Although  the  full  oxidation  of  methanol  produces 
only  C02  and  H2,  several  other  partial  oxidation  products 
can  be  generated,  such  as  formic  acid  and  formaldehyde.  An 
ideal  catalyst  would  have  maximum  selectivity  towards  full 
oxidation,  but  the  characterization  techniques  discussed 
above  can  only  determine  overall  activity.  A  scanning  differ¬ 
ential  electrochemical  mass  spectrometer,  which  enables 
spatially  resolved  (on  the  order  of  the  diameter  of  the  capil¬ 
lary)  investigations  of  the  product  phases  during  electro¬ 
chemical  catalysis,  was  devised.375  In  this  technique,  the 
rastering  electrode  was  a  hollow  capillary,  which  fed  directly 
into  a  mass  spectrometer,  and  a  hydrophobic  membrane  was 
used  to  isolate  the  electrolyte  from  the  high  vacuum  of  the 
MS.  High  sensitivity  measurements  were  demonstrated  by 
detection  of  the  oxidation  of  a  monolayer  of  CO  on  a  Pt 
surface.351 

The  effects  of  physically  mixing  different  compositions 
of  Pt-Ru  catalysts  with  various  proportions  of  Nation  (a  sul- 
fonated  tetrafluoroethylene  based  fluoropolymer  with  high 
protonic  conductivity),  through  a  multi-tiered  screen  were 
investigated.376  The  best  catalysts,  as  identified  by  cyclic 
voltammetry  and  impedance  spectroscopy  were  evaluated  in 
a  DMFC  application.  It  was  found  that  a  Pt:Ru  ratio  of  1:1 
mixed  with  36.4  wt.  %  Nation  exhibited  the  highest  activity 
and  the  lowest  initial  activation  and  CO  desorption  poten¬ 
tials;  a  synergistic  effect  between  the  Nation  and  the  active 
metal  was  proposed  that  involves  Nation  providing  protonic 
sites  to  promote  CO  oxidation. 

Attempts  to  further  increase  current  densities,  lower  the 
onset  potential,  and  stabilize  the  anodes  to  the  acidic  envi¬ 
ronment  of  the  DMFC  have  resulted  in  several  investigations 
of  ternary  and  quaternary  additions  to  Pt-Ru.  An  array  of  64 
individually  addressable  electrodes  was  implemented  to 
measure  structure-composition-activity-stability  maps  of  the 
Pt-Ru-Co  ternary  system.352  The  results  indicated  that  thin 
him  catalysts  with  the  nominal  composition  of  Pt18Co62Ru2o 


exhibited  a  10  to  30  times  improvement  in  activity  over  pure 
Pt,  with  negligible  evidence  of  chemical  attack.  A  thin  him 
catalyst  with  similar  composition  (Pti5Co5oRu38)  was  also 
reported.350  The  proposed  catalyst  exhibited  current  densities 
2.5  times  those  of  conventional  Pt90Ru10,  in  addition  to  a 
favorable  onset  potential  of  0.22  V,  compared  to  a  saturated 
calomel  electrode.  The  discovery  of  a  Pt44Rui2W44  catalyst 
formulation350  which  exhibited  current  densities  1.5  times 
those  of  the  best  binary  Pt-Ru  catalysts  was  also  reported; 
however,  the  high  W  content  in  the  alloy  tended  to  make  it 
unstable  in  an  acidic  environment.  Further,  a  scanning  elec¬ 
trochemical  microscope  was  used  to  analyze  activity  of 
Pt-Ru-Mo  catalysts  for  electro-oxidation  of  hydrogen  and 
CO.377  Fluorescence  imaging  was  used  to  determine  the 
activity  of  Pt77Rui7Mo4W2,  a  promising  quaternary  catalyst 
alloy  for  DMFC  applications.  This  catalyst  was  then  bench- 
marked  against  a  standard  Pt50Ru5o  catalyst  formulation,353 
and  showed  roughly  twice  the  current  density  after  one  hour 
of  exposure  to  methanol.  The  Pt-Ru-Ni-Ti  quaternary  system 
was  also  investigated,  with  particular  attention  to  mitigating 
the  effects  of  surface  morphology  and  segregation  effects.378 
A  low  Pt  content  catalyst  (Ni3iZr13Pt33Ru23)  was  discov¬ 
ered378  which  exhibited  similar  current  densities  to  Pt84Ru16 
catalysts,  with  an  onset  voltage  that  was  150  mV  lower, 
Fig.  78.  XPS  analysis  revealed  a  substantially  shifted  Pt 
binding  energy,  possibly  reflecting  a  change  in  the  electron 
orbitals  directly  responsible  for  the  catalytic  behavior  typi¬ 
cally  observed  in  Pt. 

Arrays  of  different  noble  and  TMs  on  a  variety  of  sup¬ 
ports379  were  studied  to  investigate  the  effect  of  metal/oxide 
combinations  on  methanol  decomposition  and  the  water  gas 
shift  reaction  (i.e.,  the  reaction  of  CO  and  H20  to  form  C02 
and  H2).  It  was  found  that  Pd/Mn02  performed  the  best  for 
both  reactions.  Further  ranking  of  the  catalyst/support  com¬ 
binations  revealed  discrepancies  between  reactions,  suggest¬ 
ing  incomplete  understanding  of  the  underlying  mechanism. 
A  64  element  half-cell  array  was  used  to  investigate  the 
effect  of  Pt  loading  and  particle  size  on  carbon  supports.380 
It  was  found  that  smaller  particle  sizes  resulted  in  fine  struc¬ 
ture  in  the  voltammetry  data  that  were  not  present  in  samples 
with  larger  particle  sizes,  most  likely  associated  with  an 
increase  in  edge  sites  on  the  small  particles.  Physical  mix¬ 
tures  of  Au/Ti02,  a  high  activity  CO  oxidation  catalyst,  with 
PtRu/C  were  investigated  to  increase  methanol  reduction 
activity.381  A  reduced  oxidation  potential  as  compared  to 
pure  PtRu/C  was  observed  for  proportions  of  PtRu:Au/Ti02 
of  9:1,  8:2,  and  7:3.  Further,  electrochemically  synthesized 
mixtures  of  Pt-W03  were  studied  in  a  combined  synthesis 
and  characterization  high  throughput  system.382  Pt-Zn  alloys 
were  characterized,  and  it  was  found  that  Zn  in  Zn-rich 
alloys  leached  until  a  PtZn  phase  with  a  presumed  Pt  rich 
skin  was  formed.383  There  has  only  been  sparse  work 
devoted  to  the  cathode  (oxygen  reducing)  side  of  the  poly¬ 
mer  fuel  cell  membrane.384  One  study  was  devoted  to 
screening  the  Pt-Ni-Zr  ternary  system  for  oxygen  reduction 
reaction  properties.385  In  a  related  work,  5  Pt  group  metals 
(Pt,  Ru,  Os,  Ir,  and  Rh)  systems  were  screened  for  their 
bifunctional  oxygen  reduction  and  water  oxidation  proper¬ 
ties;  it  was  demonstrated  that  replacing  most  of  the  Ir  in  Ptlr 
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FIG.  78.  (a)  Plot  comparing  the  current  density  of  a  Ni-Zr-Pt-Ru  catalyst  as 
a  function  of  temperature  for  a  catalyst  composition  discovered  by 
Whitacre.  (b)  Same  data  normalized  to  mole  fraction  Pt  per  unit  area. 
Reproduced  with  permission  from  J.  F.  Whitacre  et  al.,  J.  Electrochem.  Soc. 
152,  A 1780  (2005).  Copyright  2005  The  Electrochemical  Society,  Inc. 

alloys  with  Ru  improved  catalyst  activity.386  A  coupled 
theoretical-experimental  approach  was  recently  demon¬ 
strated  by  Lee  et  al.  in  discovering  novel  binary  oxygen 
reduction  reaction  catalysts.387  By  calculating  the  oxygen 
adsorption  energy  for  a  series  of  compounds  and  comparing 
them  to  the  value  for  Pt,  a  new  Pd-Co  alloy  was  predicted 
and  subsequently  experimentally  verified. 

Finally,  ceria  (Ce203)  based  mixed  oxide  conductors 
have  been  under  extensive  investigation  over  the  past  10  yr  as 
potential  combined  catalysts/anodes  for  SOFC’ s. 388  The  goal 
is  to  replace  Ni/Zr02  anodes,  which  suffer  from  coke  forma¬ 
tion  during  hydrocarbon  reforming,  which  plugs  the  SOFC 
and  causes  power  loss.  Metal  doping  of  Ce203  based  catalysts 
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and  its  effect  on  methane  reforming  has  been  studied.  ’ 
Ink  jet  printing  of  oxide/carbonate  suspensions  was  used  to 
create  libraries  in  the  La-Sr-Fe-Mn-0  system,  with  particular 
emphasis  on  identifying  new  conducting  perovskite  phases.390 
Subsequent  ionic  conduction  measurements  were  attempted, 


but  significant  contamination  of  the  perovskite  surfaces 
prevented  definitive  correlations  between  composition  and 

orn  o Q 1 

ionic  conductivity  from  being  established.  ’  Although  it 
is  well-known  that  several  companies  developed  high- 
throughput  techniques  to  screen  potential  fuel  cell  materials, 
there  have  been  few  reports  in  the  literature  regarding  their 
efforts.  One  published  example  involved  the  development  of 
a  high-throughput  method  for  synthesizing,  evaluating,  and 
characterizing  new  low  temperature  SOFC  materials.392 

VII.  HIGH-THROUGHPUT  MAPPING  OF  PHASE 
DIAGRAMS 

As  discussed  in  Sec.  IV  C,  mapping  of  compositional 
phase  diagrams  is  an  integral  part  of  any  materials  explora¬ 
tion  effort.  Not  surprisingly,  combining  information  regard¬ 
ing  distribution  of  structural  phases  across  phase  diagrams 
with  mapping  of  functional  properties  represents  a  particu¬ 
larly  powerful  way  to  rapidly  identify  new  compositions  and 
compounds  with  enhanced  physical  properties.  In  the  fields 
of  metallurgy  and  ceramics,  mapping  binary  and  higher- 
order  phase  diagrams  through  bulk  compound  synthesis  is  a 
traditional  and  critical  task.  However,  despite  concerted 
global  efforts  by  materials  scientists,  a  large  fraction  of  all 
possible  ternary  and  multinary  alloys  remain  experimentally 
unexplored.  In  fact,  the  history  of  materials  science  is  paved 
with  efforts  to  understand  and  map  phase  diagrams.  The 
high-throughput  thin  film  library  approach,  combined  with 
scanning  XRD  experiments  represents  a  tremendous  oppor¬ 
tunity,  as  had  been  recognized  as  long  as  about  50  years 
ago.12,393  Bulk  high-throughput  efforts  such  as  diffusion 
multiples136  and  directed  laser  deposition394  have  also  been 
successfully  implemented  to  map  structural  phases  in  com¬ 
position  gradient  samples. 

There  have  been  many  demonstrations  of  rapid  mapping 
of  structural  phases  using  synchrotron-based  XRD.  Chu 
et  al.395  studied  epitaxially  grown  CCS  library  films  in  the 
Co-Mn-Ge  system.  It  has  also  been  used  to  map  structural 
phases  across  a  Ni-Fe-Cr  composition  spread;396  comparison 
of  the  measured  phase  distribution  from  a  CCS  library  sam¬ 
ple  annealed  at  800  °C  corresponds  well  to  the  equilibrium 
diagram  calculated  for  the  same  temperature,  as  shown  in 
Fig.  79.  The  synchrotron  XRD  mapping  exercise  is  particu¬ 
larly  useful  when  combined  with  x-ray  fluorescence  for  com- 
position  determination.  ’  ’  The  same  combination  of 

techniques  (although  not  using  a  synchrotron)  has  been 
effectively  used  to  map  composition  regions  of  metallic 
glasses  and  amorphous  metal  systems.399 

Given  the  speed  with  which  structural  data  can  be 
acquired  in,  for  example,  the  aforementioned  synchrotron 
experiments,  the  rapid  analysis  of  large  amounts  of  raw 
XRD  data  becomes  a  challenge.  There  have  been  a  number 
of  efforts  directed  at  developing  algorithms  that  allow  rapid 
identification  of  key  trends  and  features  in  diffraction  data 
taken  across  multiphase  CCS  libraries,  and  subsequent  attri¬ 
bution  of  such  data  to  individual  phases.400^07  By  perform¬ 
ing  cluster  analysis,  similar  diffraction  patterns,  even  among 
hundreds  of  spectra,  can  be  rapidly  identified.  Similar  XRD 
patterns  usually  indicate  the  presence  of  similar  phases  in 
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FIG.  79.  (a)  Calculated  ternary  diagram 
of  Ni-Fe-Cr  system  at  875  °C  and  (b) 
measured  ternary  diagram  from  a  CCS 
library  film.  Reproduced  with  permis¬ 
sion  from  A.  Rar  et  al.,  Meas.  Sci. 
Technol.  16,  46  (2005).  Copyright  2005 
Institute  of  Physics. 


a  particular  composition  region.  Thus,  one  goal  of  these 
algorithms  is  to  automatically  identify  mixtures  of  crystallo¬ 
graphic  phases,  even  when  significant  diffraction  peak  shift¬ 
ing,  for  example  due  solid  solution  substitutional  effects,  is 
present.  These  analysis  techniques  are  extremely  useful  in 
directly  and  visually  delineating  the  composition-structure- 
property  relationships  of  materials  systems.  Many  of  the 
algorithms  developed  for  handling  large  amount  of  XRD 
data  are  also  applicable  to  rapid  analysis  of  other  spectra 
such  as  Raman,  XPS,  and  FTIR  data.408 

As  a  general  comment,  computational  combinatorial 
approaches  have  recently  received  much  attention  with  the 
arrival  of  the  Materials  Genome  Initiative  (MGI).409  To  take 
advantage  of  the  power  of  high-throughput  methodologies, 
an  integrated  strategy  combining  computational  methods 
with  rapid  experimental  verification  work  is  required.  Thus, 
going  forward,  the  challenge  for  combinatorial  methodology 
is  the  effective  coupling  of  synthesis,  characterization,  and 
theory.410-414  This  will  almost  certainly  require  an  entirely 
new  type  of  flexible  management  system  for  “big  data,” 
which  can  rapidly  handle  large  amount  of  data  in  a  variety  of 
formats. 

Accurate  stmctural  phase  mapping  can  directly  facilitate 
rapid  discovery  of  functional  materials,  whose  interesting  prop¬ 
erties  are  derived  from  phase  transformations.  Morphotropic 
phase  boundary  piezoelectrics,  discussed  in  Sec.  Ill,  are  a  good 
example.  Further,  SMAs  undergo  stmctural  transitions  as  a 
result  of  martensitic  phase  transformations,  going  from  high 
temperature,  higher  symmetry  parent  phases  to  low  tempera¬ 
ture,  lower  symmetry  martensites.  The  technological  applica¬ 
tions  of  many  SMAs  depends  on  operating  the  alloys  in  the 
temperature  range  of  their  phase  transformations;  for  example, 
superelasticity,  a  property  that  derives  from  stress-induced 


transformations,  can  be  used  to  accommodate  large  stains  in 
mechanical  framework  for  a  variety  of  applications  including 
dental  braces,  peripheral  vascular  stents,  and  eyeglass 
frames  415  In  one  study,  a  combinatorial  investigation  was  car¬ 
ried  out  to  identify  those  SMA  compositions  with  minimal 
martensitic  transformation  thermal  hysteresis  416  It  is  known 
that  minimizing  the  hysteresis  would  increase  the  fatigue  life  of 
SMAs,  and  thus,  significantly  improve  the  long-term  perform¬ 
ance  of  many  SMA  applications.  Guided  by  the  non-linear 
theory  of  martensite  developed  by  Ball  and  James  417  combina¬ 
torial  mapping  of  the  lattice  constants  of  Ni-Ti-Cu  SMAs  was 


FIG.  80.  Thermal  hysteresis  (measured  by  resistive  transitions)  vs.  middle 
eigenvalues  (determined  from  lattice  constants)  for  different  Ni-Ti-X-Y  sys¬ 
tems  shows  that  making  the  middle  eigenvalue  go  to  1  (by  satisfying  lattice 
constant  conditions  set  by  the  non-linear  theory  of  martensite)  reduces  the 
thermal  hysteresis  of  martensitic  transformation.  Each  data  point  corre¬ 
sponds  to  a  different  composition.  Reproduced  with  permission  from  R. 
Zarnetta  et  al.,  Adv.  Funct.  Mater.  20,  1917  (2010).  Copyright  2010 
Wiley— VCH  Verlag  Gmbh  &  Co. 
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(a) 


Cluster  analysis 


FIG.  81.  (a)  Rapid  cluster  analysis  of  diffraction  spectra  taken  from  a  Ni-Ti-Cu  composition  spread.  Different  colors  indicate  spectra  with  similar  patterns 
reflecting  the  different  structural  regions.  Reproduced  with  permission  from  R.  Zarnetta  et  al.,  Intermetallics  26,  98  (2012).  Copyright  2012  Elsevier  BV.  (b) 
High-temperature  phase  diagram  of  Ni-Ti-Cu  obtained  using  bulk  samples.  The  clear  resemblance  between  (a)  and  (b)  attest  to  the  validity  of  the  composition 
spread  and  the  cluster  analysis  technique  for  quick  delineation  of  structural  phase  distributions  across  large  compositional  ranges  (see  Ref.  438).  Reproduced 
with  permission  from  F.  J.  J.  Vanloo  et  al.,  J.  Less-Common  Met.  57,  111  (1978).  Copyright  1978  Elsevier  Publishing  Co. 


carried  out  through  synchrotron  XRD.  This  study  revealed 
that  satisfying  the  lattice  constant  conditions  set  forth  by  the 
theory  indeed  leads  to  SMAs  with  minimal  hysteresis.  In  a 
follow-on  study,  the  quaternary  system  Ni-Ti-Cu-Pd  was  inves¬ 
tigated  to  further  improve  the  composition,  and  focus  on 
“zero-hysteresis”  compounds,  as  determined  by  thin  film  meas¬ 
urements.418  Figure  80  summarizes  this  research;  it  can  be 
seen  that  minimization  of  the  thermal  hysteresis  was  achieved 
through  tailoring  the  lattice  constants  of  the  alloys.  As  part 
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FIG.  82.  Schematic  of  the  parallel  nano-scanning  calorimeter  (see  Ref.  439). 
Reproduced  with  permission  from  P.  J.  McCluskey  et  al.,  Acta  Mater.  59, 
5116  (2011).  Copyright  2011  Pergamon. 


of  the  same  study,  it  was  shown  that  the  “zero-hysteresis” 
compositions  identified  from  the  CCS  experiments,  when  proc¬ 
essed  as  bulk  samples,  indeed  showed  substantially  enhanced 
functional  stability  due  to  minimized  fatigue  behavior.  Zametta 
et  al 419,420  have  also  investigated  the  details  of  phase  transfor¬ 
mation  behavior  across  a  variety  of  related  SMA  systems. 
Figure  81  shows  one  such  study  where  data  cluster  analysis 
was  performed  on  XRD  spectra  from  a  Ni-Ti-Cu  CCS  library. 
The  clear  resemblance  to  the  known  bulk  high-temperature 
phase  diagram  (Fig.  81(b))  speaks  to  the  effectiveness  of  the 
combinatorial  approach.421 

Other  methods  of  mapping  phase  diagrams  of 
SMAs  involve  their  mechanical  properties.  Since  in  their  mar¬ 
tensitic  state,  they  are  mechanically  soft,  nanoindentation 
mapping  of  elastic  modulus  has  been  effectively  used  to  iden¬ 
tify  martensite  composition  regions  422  Another  method  for 
detection  of  the  martensitic  transition  is  through  implementa¬ 
tion  of  arrays  of  cantilevers,  since  the  change  in  the  elastic 
modulus  of  the  film  deposited  atop  the  cantilever  as  a  result  of 
the  martensitic  transition  results  in  an  abmpt  deflection.108 
Indeed,  cantilever  arrays  have  served  as  a  versatile  combinato¬ 
rial  platform  for  mapping  a  variety  of  mechanical  properties 
such  as  film  stresses,  glass  transitions,  and  fatigue.311,423^126 
Finally,  other  dense  MEMS  sensor  arrays,  such  as  microhot¬ 
plates  for  nanocalorimetry  and  gas  sensors,  have  been  very 
useful  in  high  throughput  measurements.427^129  Figure  82 
illustrates  a  nanocalorimeter  array  device. 

Another  class  of  materials  that  has  advanced  as  a  result 
of  high  throughput  mapping  of  properties  are  ternary  com¬ 
pounds  with  a  composition  of  Mn+1AXn  (“MAX”  phases) 
where  n  =  1,  2,  or  3,  M  is  a  TM,  A  is  a  III- A  or  IV- A  Group 
element,  and  X  is  carbon  or  nitrogen.  Currently,  more  than  50 
distinct  MAX  phases  have  been  synthesized  in  bulk  430  There 
is  much  interest  in  this  class  of  refractory  materials  due  to 
their  machinability  and  extreme  hardness.431^133  Further, 
since  they  exhibit  low  friction  and  oxidation  resistance,430 
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these  materials  are  good  candidates  for  protective  coatings. 
However,  much  remains  unknown  about  MAX  phases,  espe¬ 
cially  with  regard  to  the  properties  of  solid  solutions.  In 
recent  research,  nanoindentation  mapping  was  carried  out  on 
CCS  films  of  (Ti!_xNbx)2AlC  and  (Cr|=xVx)2GeC,  and  it  was 
discovered  that  some  solid  solutions  have  larger  elastic  mod¬ 
ulus  compared  to  their  end  members.434,435 
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NOMENCLATURE 

ALD  =  atomic  layer  deposition 
APCVD  =  atmospheric  pressure  chemical  vapor  deposition 
BHmax  =  maximum  energy  product 
CCD  =  charge  coupled  display 
CCS  =  continuous  composition  spread 
CMOS  =  complementary  metal  oxide  semiconductor 
CRT  =  cathode  ray  tube 
C-V  =  capacitance-voltage 
CVD  =  chemical  vapor  deposition 
DFT  =  density  function  theory 
DMFC  =  direct  methanol  fuel  cell 
EPMA  =  electron  probe  microanalysis 
FTIR  =  Fourier  transform  infrared  spectroscopy 
FYNES  =  fluorescence  yield  near  edge  spectroscopy 
ITO  =  indium  tin  oxide 

ITRS  =  International  Technology  Roadmap  for 
Semiconductors 
I- V  =  current-voltage 
IR  =  infrared 

MEMS  =  microelectromechanical  systems 
MFM  =  magnetic  force  microscopy 
MOCVD  =  metal-organic  chemical  vapor  deposition 
MOKE  =  magnetic-optical  Kerr  effect 
MPB  =  morphotropic  phase  boundary 
NIR  =  near  infrared 

NMOS  =  n-channel  metal  oxide  semiconductor 
PCM  =  phase  change  material 
PEM  =  polymer  electrolyte  membrane 
PFM  =  piezoresponse  force  microscopy 
PGAA  =  prompt  gamma  activation  analysis 
PLD  =  pulsed  laser  deposition 
PMOS  =  p-channel  metal  oxide  semiconductor 
PV  =  photovoltaic 
PZT  =  Lead  zirconium  titanate 
RBS  =  Rutherford  backscattering 
RE  =  rare  earth 
RF  =  radio  frequency 

RHEED  =  reflection  high  energy  electron  diffraction 
RT  =  Room  temperature 
SEM  =  scanning  electron  microscopy 
SEMPA  =  scanning  electron  microscopy  polarization 
analysis 


SMA  =  shape  memory  alloy 
SMM  =  scanning  microwave  microscopy 
SOFC  =  solid  oxide  fuel  cell 
SQUID  =  superconducting  quantum  interference  device 
SSM  =  scanning  SQUID  microscopy 
STO  =  strontium  titanate 
TCO  =  transparent  conducting  oxide 
TDTR  =  time  domain  thermoreflectance 
TEM  =  transmission  electron  microscopy 
TIO  =  tin  indium  oxide 
TM  =  transition  metal 
UV  =  ultraviolet  light 

VASP  =  Vienna  Ab  Initio  Simulation  Program 
VIS  =  visible  light 

WDS  =  wavelength  dispersive  spectroscopy 
XMCD  =  X-ray  magnetic  circular  dichroism 
XPS  =  x-ray  photoemission  spectroscopy 
XRD  =  x-ray  diffraction 
ZITO  =  zinc  indium  tin  oxide 
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